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The  phenomena  of  physical  sputtering  have  been 
extensively  studied  since  it  was  first  mentioned  in 
Philosophical  Transactions  by  Grove  in  1352.  Sputtering  is 
the  ejection  of  atomic  material  from  a target  resulting  from 
bombardment  by  an  energetic  particle.  Most  scientific 
investigation  has  been  restricted  to  ionic  bombardment 
because  the  ion  can  easily  be  accelerated  to  any  desired 
velocity.  Sputtering  has  defied  complete  physical 
description  because  of  complicated  mechanisms,  difficulty  in 
treating  the  behavior  of  a many-body  system,  and  an 
incomplete  understanding  of  the  interacting  forces  and 
potentials  in  crystal  structures. 

Early  experiments  typically  determined  the  "sputtering 
ratio"  of  single  and  polycrystalline  materials.  The 
"sputtering  ratio"  (interchangeable  with  sputtering  yield 
Y)  is  defined  as  the  number  of  target  atoms  ejected  from  the 
surface  per  incident  ion. 

In  1923,  Kingdon  and  Langmuir  [1]  bombarded  thoriaced 
tungsten  with  various  ions  in  a glow  discharge  tube.  This 
was  a special  case  of  sputtering  since  the  thin  surface  film 
of  thorium  on  a tungsten  substrate  was  sputtered  rather  than 
the  tungsten  itself.  The  results  of  this  experiment 
qualitatively  suggested  an  ejection  mechanism. 

A few  years  later,  Blechschmidt  and  von  Hippie  [2,3] 
proposed  a theory  which  described  sputtering  as  an 
evaporation  of  surface  atoms  (an  emission  mechanism). 
Earlier,  von  Hippie  [4]  had  found  by  spectroscopic  analysis 
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that  at  least  soma  sputtered  atoms  were  in  an  excited  state. 
The  sputtering  theory  showed  that  atoms  in  the  region  of 
impact  could  rapidly  acquire  thermal  energy  if  the  kinetic 
energy  of  the  incident  ion  was  converted  to  thermal  energy 
at  the  target  surface.  If  an  atom  acquired  a sufficient 
amount  of  thermal  energy,  it  would  then  evaporate  from  the 
surface,  some  atoms  evaporating  while  in  an  excited  state. 

Direct  application  of  statistical  methods  to  sputtering 
was  made  by  Harrison  [5]  who  envisioned  the  interaction  of 
two  distribution  functions  (the  crystal  lattice  and  the  ion 
beam) . These  models,  based  on  statistical  methods, 
implicitly  accept  ejection-type  mechanisms. 

One  of  the  most  important  contributions  to  the  study  of 
sputtering  was  made  by  Wehner  [6]  in  1953.  He  discounted 
the  evaporation  theory  and  presented  strong  evidence  for  a 
momentum  transfer  process.  Shortly  after  Wehner's  findings 
were  reported,  Henschke  [7]  proposed  a theory  of  sputtering 
based  solely  on  classical  theory,  treating  normal  and 
oblique  incidence  sputtering  separately.  His  concept  of 
normal  incidence  sputtering  required  many-body  collisions  in 
which  the  ion  was  eventually  reflected  outward  to  sputter 
surface  atoms.  This  theory  was  plausible  for  oblique 
incidence  sputtering;  however,  the  case  for  normal  incidence 
sputtering  required  that  the  ion  be  rerlected  inside  the 
crystal . 

Wehner  and  Rosenberg  [8]  reported  that  in  the  sputtering 
of  polycrystalline  metals  at  low  ion  energy  (normal 
incidence)  more  material  is  ejected  obliquely  than  in  a 
direction  normal  to  the  target  surface.  This  is  due  to  the 
fact  that  less  than  a full  (180°)  overall  direction  change 
of  momentum  is  more  favorable  for  sputtering  because  this 
requires  fewer  number  of  successive  collisions  for  atom 
ejection.  In  1960,  Wehner  and  Anderson  [9]  reported 
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sputtering  patterns  from  single  crystal  targets  which  showed 
an  anisotropic  spread  of  energy  from  a collision  center  and 
atoms  which  are  preferentially  ejected  in  close- packed 
directions. 

When  an  atom  collides  with  another  which  has  a higher 
mass,  the  scattering  angle  can  vary  between  0°  and  180°.  In 
a collision  with  a lower  mass  atom  the  scattering  is 
confined  to  between  0°  and  90°  (thus  backward  reflection 
cannot  occur) . In  low-energy  sputtering  of  multi-element 
targets.  Winters  and  Sigmund  [10]  showed  theoretically  that 
a significant  contribution  to  sputtering,  called  "single 
reflective  collisions",  arises  from  the  reflection  of  lower 
mass  atoms  from  underlying  heavier  atoms.  For  a detailed 
discussion  of  momentum  reversal,  see  Harrison  and  Magnuson 
[ 11  ]. 

In  1962,  Wehner  et  al  [12]  published  a compilation  of 
sputtering  yields  for  metals  and  semiconductors  in  the 
100-600  eV  range.  In  1966,  Wehner,  Anderson,  and  KenKnight 
[13]  published  the  results  of  a study  that  covered  almost  a 
decade  of  research  in  low-energy  sputtering.  A few  of  the 
numerous  topics  discussed  included  more  accurate  and 
improved  techniques  for  analyzing  the  sputtered  mass,  the 
temperature  dependence  of  sputtering  yields  of  variolous 
metals  bombarded  by  atomic  and  polyatomic  ions  in  the  2.8  to 
10  keV  range,  the  effect  of  oxygen  on  ion  ejection  patterns, 
and  the  spatial  distribution  of  ejected  atoms  in  metals  and 
semiconductors. 

A few  low-energy  sputtering  experiments  [14,  15,  16,  17, 
13],  which  lend  themselves  to  direct  comparison  with  these 
simulations,  investigate  the  behavior  of  the  sputtering 
ratio  as  a function  of  ion  energy,  lattice  orientation, 
target  material,  and  ion  angle  of  incidence. 
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Paralleling  the  experimental  investigation  of  the 
sputtering  yield  and  the  development  of  mass  spect rogra phic 
analysis  was  the  attempt  to  model  the  sputtering  mechanism 
through  the  use  of  computer  simulation.  In  1950,  Gibson, 
Goland,  Milgram,  and  Vineyard  [19],  the  first  investigators 
to  use  computer  simulation  as  a means  of  studying  radiation 
damage,  set  forth  the  criteria  that  must  be  satisfied  before 
the  simulation  method  could  be  applied  to-  crystals. 
Reference  [19]  gives  a very  detailed  description  of  the 
simulation  method.  Basically,  they  built  a computer  model 
to  represent  metallic  copper  and  studied  radiation  damage 
events  at  low  and  moderate  ranges  of  energy  up  to  400  eV. 
In  their  model,  one  atom  in  the  stationary  lattice  was  given 
an  arbitrary  kinetic  energy  and  direction  of  motion,  as 
though  it  had  been  struck  by  a bombarding  particle,  and  was 
allowed  to  interact  with  the  remaining  atoms  in  the  lattice. 
This  resulted  in  cascades  of  independent  binary  collisions. 
It  is  well  established  that  sputtering  of  solids  by  ion 
bombardment  is  the  result  of  successive  independent 
approximately  binary  collisions  somewhat  as  if  these 
collisions  take  place  with  or  between  separate  gas  atoms. 
Low-energy  ion  scattering  studies  [20]  have  shown  that  the 
approximately  binary  scattering  behavior  still  holds  at 
bombarding  energies  as  low  as  those  approaching  the 
sputtering  thresholds. 


In  1967,  Harrison,  Levy,  Johnson,  and  Sffron  [21]  used  a 
computer  to  simulate  the  bombardment  of  f ace-centered-cubic 
(fee)  copper  single  crystals  by  an  argon  ion.  From  their 
studies,  the  sputtering  mechanisms  in  the  fee  crystal  were 
isolated  and  identified.  They  reported,  for  ion  energies 
less  than  10  keV,  the  sputtering  process  occurred 
predominately  within  three  atomic  layers  of  the  surface. 
The  computer  simulation  included  only  a repulsive  force 
between  the  atoms.  The  equations  of  motion  for  the 
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many-body  system  were  computed  using  an  "average- force" 
algorithm  rather  than  a "central-difference"  algorithm,  see 
Harrison,  Gay,  and  Effron  [22].  In  1972,  Harrison,  loore, 
and  Holcombe  [23]  in  continuing  efforts  to  develop  a more 
realistic  computer  .model  incorporated  a Cu-Cu  potential 
function  which  included  an  attractive  potential,  surface 
layer  relaxation,  and  surface  binding  energy  for  the  (100), 
(110)  , and  (111)  copper  planes.  They  reported  improved 
agreement  between  experimental  and  simulation  results,  but 
the  original  interpretation  of  the  low-energy  sputtering 
mechanisms  in  fee  copper  crystals  remained  unchanged. 

Computer  simulation  of  physical  processes  is  still  in  a 
very  infant  state  because  the  available  processing  speed  of 
existing  computers  severely  limits  the  sophistication  of  the 
models.  However,  past  investigations  were  quite  successful 
in  identifying  the  mechanisms  which  occur  when  energetic 
particles  interact  with  crystal  atoms.  The  ability  to 
observe  these  mechanisms  is  an  advantage  peculiar  to  tie 
computer  simulation  because  each  crystal  atom  is  identified 
by  a number  for  the  mathematical  calculation,  and  its 
complete  track  can  be  labeled  and  recorded.  The  tracks  of 
selected  atoms  can  be  examined  to  determine  the  dynamics  of 
the  atoms  which  exit  the  front  surface. 


II . 


RESEARCH  OBJECTIVES 


A.  ANGULAR  DISPERSION  OF  ISOTOPES 

Recently,  Wehner  [24]  reported  that  the  isotope  ratio  in 
metal  deposits  obtained  by  spattering  a flat  target  at  low 
ion  energy  (100e7)  under  normal  incidence  is  a function  of 
the  ejection  angle.  He  investigated  the  bombardment  of  a 
Cu  (63-65)  target  with  100  eV  Hg  ions  . He  reported  that  the 
lighter  isotopes  are  enriched  in  a direction  normal  to  the 
target  surface.  This  research  was  a computer  simulation 
investigation  of  Wehner' s studies  of  angular  dispersion  of 
isotopes  in  a monocrystalline  copper  lattice. 


B.  LIQUID  SPLASH  EFFECT 

In  preliminary  NPS  simulation  studies,  after  analysis  of 
data,  a predominance  of  sputtered  atoms  occurred  at  certain 
distances  from  the  center  of  the  lattice  representative  area 
in  which  the  ion  impacts  occur.  Figures  1 and  2 f 25  ] depict 
the  top  layer  of  a copper  fee  crystal  respectively  bombarded 
by  neon  and  gold  ions  at  2 keV . The  numbers  represent  the 
number  of  times  that  the  atom  was  sputtered  in  36  events 
(shots) . The  dark  triangle  in  the  center  represents  the 
target  impact  area.  All  36  events  are  directed  into  this 
representative  area.  The  darkened  elipses  highlight  the 
atoms  most  often  sputtered  and  suggest  a ring  outward  from 
the  center  (especially  for  the  Cu-Au  interaction). 


21 


Macklin  and  Metaxas  [ 26  ] experimentally  investigated  the 
splashing  of  liquid  drops  on  deep  and  shallow  liquid  layers. 
They  published  high  speed  photographs  which  showed  liquid 
splashes  that  might  be  analogous  to  surface  atoms  being 
ejected  from  the  surface  and  a few  eventually  sputtering. 

This  research  was  an  investigation  of  "ring  formation" 
or  "liquid  splash"  effect  of  sputtered  atoms.  Determination 
was  made  of  the  probability  that  a congregate  of  atoms  at  a 
given  distance  from  the  impact  area  reference  will  sputter. 


C.  ENERGY  CONTAINMENT 


As  with  any  computer  simulation,  a trade-off  must  be 
made  between  the  accuracy  and  detail  desired  and  the  size 
and  time  requirements  imposed  by  computer  availability. 
Because  of  this  constraint  a limit  exists  on  the  size  of  the 
monocrvstalline  lattice  used  in  the  simulation.  This 
research  was  an  investigation  to  determine  the  size  lattice 
required  to  contain  a low  energy  event. 


D.  ION/ATOM  MASS  RATIO 

In  July  of  1976,  Eer  Nisse  [27],  reported  experimental 
confirmation  of  a monotonic  dependence  of  Y with  the  atomic 
number.  He  investigated  the  yield  of  a gold  target  impacted 
by  hO  ke7  ions  up  to  the  atomic  number  of  gold.  This 
monotonic  dependence  is  predicted  by  a sputtering  theory  of 
Sigmund  [28].  However,  computer  simulation  data  of  previous 
N PS  studies  [29]  indicate  that  the  monotonic  dependence  of 
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the  sputtering  ratio  on  the  ion  mass  ceases  to  be  valid  for 
ion/atom  mass  ratios  greater  than  one.  This  research  was  an 
investigation  of  the  sputtering  ratio  versus  ion  atomic 
number  for  a wide  range  of  ion/atom  ratios. 


III.  SIMULATION  MODEL  MODIFICATIONS 


The  basic  simulation  model  used  in  this  research  has 
been  developed  by  the  NPS  group  over  the  last  decade  [21, 
23,  31,  32  ].  A brief  description  of  this  model  and  the 
alterations  made  to  it  for  this  research  follow.  For 
detailed  discussions  of  the  simulation  method,  the  reader  is 
directed  to  references  19,  20,  21,  and  23.  A detailed 
description  of  the  mathematical  model  has  been  published 
[22]. 

In  the  development  of  the  computer  model,  the  lattice 
unit  (LU)  proves  to  be  a convenient  unit  of  length,  and  is 
used  extensively.  One  LU=a„/2,  where  a0  is  the  cubic 
lattice  constant,  the  atomic  spacing  as  determined  by  x-ray 
crystallographic  studies.  For  copper,  a0  =3.615  angstroms. 

Certain  properties  of  the  program  and  results  are  best 
considered  in  the  framework  of  a specific  run.  Each  run  is 
begun  by  constructing  a target  microcrystallite  whose  sites 
represent  the  equilibrium  positions  of  copper  atoms  near  the 
surface  of  a fee  crystal.  The  target  atom  - target  atom 
interaction  is  represented  by  the  3orn-Mayer  Gibson  type  II 
(Gibson-II)  repulsive  potential  function. 

Each  ion  of  the  beam  is  represented  by  a single  neutral 
atom  of  the  desired  element  [30].  The  ion  is  placed  at  "-he 
point  on  its  trajectory  such  that  its  velocity  vector 
determines  the  point  of  impact.  The  term  bullet  and  ion  are 
used  interchangeably  throughout  this  report  for  the  incident 
particle  to  avoid  confusing  it  with  target  atoms. 
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The  bullet-target  potential  function 
the  Gibson-II  or  that  determined  by  the 
Electron  (KSE-B)  Potential  [30]-  The 
functions  seem  to  represent  a reasonable  r 
"soft"  potential  functions. 

As  reported  by  Harrison,  Moore,  and  Ho 
attractive  forces  between  the  atoms  in 
small,  when  compared  to  the  forces  whic 
moderately  energetic  collisions,  that 
neglected  in  sputtering  simulations.  Also 
the  smaller  reduced  target  areas,  the  ion 
incident.  Even  though  the  computer  progra 
the  capability  to  utilize  attractiv 
obliquely  incident  ions,  they  were  no 
investigation . 
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A.  MICROCRYSTALLITE  SIZE 


Initial  discussio 
revealed  the  advanta 
lattice  that  could  be 
sputtered  by  ions  i 
would  become  available 
layer  a square  and  ha 
boundary  was  needed, 
previously  used  co 
alterations. 


n to  establish  research  objectives 
ges  gained  from  using  a symmetrical 
folded  and  rotated  so  that  atoms 
ncident  at  all  possible  impact  points 
for  analysis.  A lattice  with  each 
ving  an  odd  number  of  atoms  alon j each 
The  lattice  generation  program  as 
uld  handle  this  change  with  no 


Figures  3,  4,  and  5 show  a roughly  isometric  projection 
of  the  labeled  atoms  in  a 9x4x9  (100)  orientation,  a 9x5x9 
(110)  orientation,  and  a 9x4x9  (111)  orientation  of  a fee 
microcrystallite  respectively.  This  type  of  projection, 
first  used  by  Harrison  and  Delaplain  [31],  with  the  distance 
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between  layers  in  the  Y direction  exaggerated,  and  the  atoms 
represented  by  20°  ellipses  has  proved  to  be  a valuable  tool 
for  analyzing  the  computer  programs  and  for  displaying  data. 
The  first  number  in  each  lattice  is  atom  number  two,  with 
atom  number  one  being  the  bullet. 

From  previous  studies,  it  was  clear  that  at  the  low 
energy  levels  considered  the  sputtering  mechanisms  were 
confined  to  the  surface  layers  [21].  For  this  reason,  the 
(100)  and  (111)  orientation  microcrystallites  were  only  four 
layers  deep.  The  (110)  orientation  lattice  generation 
program  constructs  a lattice  which  has  significantly  fewer 
atoms  and  its  layers  are  closer  together.  In  an  attempt  to 
increase  the  number  of  atoms  in  the  lattice  while  keeping 
the  required  computer  time  below  that  of  the  (100)  and  (111) 
orientation  programs,  the  (110)  orientation  microcr ystallite 
was  set  at  six  layers  deep. 

3.  CRYSTAL  BOUNDARIES 


The  previous  simulation  program  contained  a 
microcrystallite  with  established  boundaries  that  were  not 
consistent  with  the  development  of  the  potential  functions. 
Since  the  potential  function  is  truncated  at  the 
nearest- neighbor  separation  (ROE) , all  boundaries  (except 
the  top  surface)  were  set  at  ROE.  This  is  justified  by  the 
fact  that  an  atom  outside  these  boundaries  is  farther  from 
the  crystal  than  the  potential  truncation  radius.  The  upper 
surface  boundary  was  set  by  starting  a small  distance  from 
the  lattice  surface,  and  increasing  the  distance  in 
succeeding  runs  until  the  total  yield  was  independent  of  the 
distance.  For  computational  efficiency  this  distance  should 
be  set  as  small  as  possible. 
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Figure  6 shows  a drawing  of  the  microcrystallite  in 
solid  lines  and  the  plane  boundaries  in  dashed  lines.  The 
boundaries  which  are  defined  in  the  Glossary  (Appendix  C) 
are  shown  as  planes  at  the  distance  from  the 
microcrystallite  surface  cited  above. 


C.  ISOTOPE  GENERATION 

In  order  to  study  the  angular  dispersion  of  sputtered 
isotopes,  a method  for  generating  a microcrystallite  with 
the  correct  isotope  ratio  and  masses  was  needed.  The  method 
used  assigned  an  additional  variable  to  each  atom  that 
identified  the  atom  "type".  The  atoms  were  classified 
according  tc  two  types;  type  1 - least  abundant,  heavier 
isotope  and  type  2 - most  abundant,  lighter  isotope. 
Appendix  A is  a detailed  discussion  of  how  the  correct 
isotope  ratio  was  obtained  and  each  atom  properly  identified 
with  the  correct  mass. 


D.  TARGET  IMPACT  AREAS 


A representative  area  must  contain  all  possible  impact 
points  upon  a particular  surface,  that  is  any  point  in  the 
surface  must  be  projectable  into  a point  of  the 

representative  area.  Figure  7 shows  the  intrinsic  geometry 
of  the  three  fee  orientations  and  the  complete  sets  of 
impact  points  for  each  representative  area.  Figure  8 shows 
the  reduced  sets  of  impact  points  from  which  the  complete 
sets  are  generated  by  reflection  and  rotation.  Figures  9, 
10,  and  11  show  an  expanded  view  of  the  targe4-  impact 
and  the  actual  impact  points  used.  The  individual 
points  are  those  used  in  previous  simulations  for 
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impact 

normal 


incidence.  This  simplification,  which  greatly  reduces  the 
computer  running  time,  is  only  possible  when  the  ion  is 
normally  incident. 

E.  LATTICE  FOLDING 


In  order  to  investigate  the  ring  formation  or  "liquid 
splash"  effect,  it  was  necessary  to  develop  a means  of 
folding  the  results  obtained  from  the  reduced  target  area. 
This  requirement  was  satisfied  by  writing  three  subroutines; 
KFOLDA,  KFOLDB,  and  KFOLDC.  Appendix  B gives  a detailed 
description  of  the  "folding"  logic  and  the  geometry  of  each 
orientation.  Figure  12  shows  how  the  (100)  and  (110) 
orientation  reduced  representative  areas  were  folded  to 
obtain  the  complete  representative  areas.  Figure  13  shows 
how  the  (111)  orientation  reduced  representative  area  was 
folded  and  summed  to  obtain  the  complete  representative 
area. 


F.  ATOM  £ POTTERING  PROBABILITY 

A means  of  talcing  the  folded  data  and  calculating  a 
conditional  probabilty  that  an  atom  would  sputter  from  a 
particular  ring  was  needed.  This  was  accomplished  by 
writing  the  subroutine  KPROB.  In  this  subroutine,  the 
distance  from  each  atom  to  the  reference  (0,0)  point  in  the 
reduced  representative  impact  area  was  calculated.  then  all 
atoms  at  each  distance  were  checked  and  the  atom  sputtering 
count  was  summed.  This  sum  was  divided  by  the  total  number 
of  sputs  from  all  the  shots.  This  quotient  was  the 
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conditional  probabilty  that  a given  ring  at  a known  distance 
would  sputter  an  atom  given  the  atoms  that  had  already 
sputtered . 

.*r>R03  performs  the  above  calculations.  However,  a very 
interesting  result  is  obtained  by  dividing  the  conditional 
probabilty  obtained  above  by  the  number  of  atoms  in  a 
particular  ring.  This  result  was  a conditional  probability 
that  an  atom  at  a given  distance  from  the  impact  point  would 
sputter. 


G.  EJECTION  ANGLES 


The  ejection  angle  is  defined  as  the  angle  between  the 
atom's  velocity  vector  and  the  surface  normal.  A tally  was 
kept  of  each  atom's  ejection  angle,  grouped  in  five  degree 
intervals.  This  tally  was  recorded  separately  for  each 
isotope.  iO  provide  a quick  analysis  of  the  data,  a program 
(HISTO)  was  written  to  print  out  a histogram  of  the  number 
of  sputtered  atoms  ejected  in  each  interval  for  each 
isotope . 


H.  VARIABLE-SHOT  SELECTION 


The  previous  program  was  only  flexible  enough  to  allow 
either  a single-shot  run  or  all  36  shots.  As  the  size  of 
the  lattice  and  the  ion's  energy  were  increased,  the 
computer  run  time  for  a set  of  36  shots  became  unmanageable 
(12  - 18  hours) . At  this  point,  a variable  size  shot 
"package"  was  developed.  The  features  included  the  ability 
to  specify  the  number  of  shots  desired  and  the  starting 
impact  point  of  the  "package".  The  final  summary  for  the 
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particular  package  run  was  printed  out.  The  packages  had  to 
be  combined  by  hand  calculations  in  order  to  get  the  total 
information  for  a particular  run.  Even  so,  this  procedure 
allowed  runs  to  be  made  with  larger  lattices  and  higher 
energies  (i.e.  the  containment  studies)  that  previously 
would  have  been  prohibitively  long. 


I.  RANDOMLY  SELECTED  IMPACT  POINTS 


To  obtain  more  angular  dispersion  data  for  the 
low-energy,  (100)  orientation  events,  a modification  was 

si  made  to  the  simulation  model.  This  modification  replaced 

tr 

the  systematic  method  of  incrementing  the  impact  point  over 
a grid  of  points  in  the  representative  target  area  with  a 
random-shot  capability.  Again,  the  IMSL  GGUB  routine  was 
used  to  generate  uniformly  distributed,  psuedo-random 
numners.  The  first  random  number  was  labeled  the  x 

coordinate  in  the  impact  area.  The  second  random  number  was 
labeled  the  z coordinate  and  tested  to  determine  if  it  was 
inside  the  impact  area  z boundaries  determined  by  the  value 
of  the  x coordinate.  If  this  test  was  passed  sucessfuliy, 
these  values  for  the  x and  z coordinates  were  used. 
However,  if  the  test  did  not  pass  (the  x and  z coordinates 
were  outside  the  impact  area),  the  two  original  numbers  were 
discarded,  two  more  were  selected,  and  the  process  was 
repeated  again. 

This  random  impact  point  did  not  affect  the  variable  - 
shot  selection  feature.  3ut,  enough  random  numbers  had  to 
be  generated  to  allow  for  those  points  that  did  not  fall 
into  the  impact  area. 


The  results  of  randomly  bombarding  the  impact  area  will 
be  discussed  later.  This  feature  was  incorporated  in 
response  to  other  investigators  suggestion  that  the  grid  of 
impact  points  be  compared  with  a randomly  bombarded  impact 
area . 
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IV  . 


RESULTS 


The  three  lattice  orientation  models  for  the  copper 
target  were  utilized  in  computer  simulation  runs  in  which 
100-eV,  200-eV,  500-eV,  and  1-keV  normally  incident  ions 
bombarded  the  target.  The  ion  masses  varied  from  4 amu 
(He+)  to  over  200  amu  (Hg+)  . The  lattice  was  constructed 
with  copper  atoms  of  either  the  natural  mass  (63.54  amu)  or 
the  proper  ratio  of  isotopes  (63.93  amu  and  64.93  amu)  . The 
lattice  size  varied  from  a 9x4x9,  (100)  orientation  of  162 
atoms  to  a 17x4x17,  (100)  orientation  of  576  atoms. 


Computer  runs  were  made  for  various  permutations  of  ion 
energy,  ion  mass,  lattice  orientation,  lattice  size,  and 
potential  function.  The  summary  of  each  run  contained  the 
number  of  atoms  sputtered  from  each  impact  point,  the  number 
of  times  each  atom  sputtered,  a histogram  showing  the 
angular  dispersion  of  the  sputtered  atoms  for  each  isotope, 
the  number  of  sputs  for  each  atom  after  the  folding  process, 
the  distance  from  each  ring  to  the  impact  area  reference 
point,  the  number  of  atoms  and  the  sum  of  their  sputs  in 
each  ring,  and  the  probability  that  each  ring  contained  a 
sputtered  atom. 

In  the  sputtering  summaries  the  numbers  in  the  ellipses 
indicate  the  number  of  times  that  particular  atom  sputtered 
in  the  simulation.  The  dark  areas,  a triangle  for  the  (100) 
and  (111)  orientations  and  a rectangle  for  the  (110) 
orientation,  indicate  the  location  of  the  target  impact 
area.  The  atom  that  is  struck  by  the  ion  is  called  the 
"primary  knock-on  atom"  (hereafter  PKA)  . Since  different 
seeds  were  used  to  generate  the  random  isotope  sites,  a 
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letter  L or  H denotes  whether  the  sputtered  atom  was  a 
light  or  heavy  isotope.  Wherever  possible,  the  results  for 
the  Gibson-II  and  KSE-B  potential  functions  were  combined. 
In  the  figures  displaying  the  sputtering  summary  and  the 
mechanism  tracing,  some  of  the  lower  layers  were  omitted. 
However,  all  atoms  were  included  in  all  calculations.  In 
the  100-eV  results  all  the  sputtered  atoms  are  shown. 
Wherever  a direction  in  the  lattice  is  specified,  it  is 
always  with  respect  to  that  particular  lattice  axes.  The 
top  surface  is  always  an  XZ  plane. 

A.  9X9  H IC R OCR YSTAL LITE 

Initially,  all  simulations  used  a 9x9  lattice  to 
determine  the  sputtering  ratio  as  a function  of  the  ion  mass 
at  100  eV  and  1 keV.  The  model  was  modified  to  include  the 
isotope  generator,  the  symmetrical  folding  of  the  summary, 
the  angular  dispersion  of  sputtered  isotopes,  and  the  "ring 
sputtering"  probabilities. 


1 • 1J00L  Orientation 

The  (100)  orientation  micr ocrystallite  was 
constructed  in  a lattice  consisting  of  162  atoms.  The 
lattice  was  9 atoms  in  the  X and  Z directions  and  4 atoms 
deep  in  the  Y direction,  hence  9x4x9.  Figure  3 shows  an 
isometric  projection  of  the  lattice  and  its  atom  numbering 
scheme.  Figure  9 shows  an  expanded  view  of  the  grid  of  36 
impact  points  used. 


33 


100-eV  Ions 


a . 


Of  primary  importance  was  the  investigation  of 
the  angular  dispersion  of  isotopes  reported  by  Wehner  [24]. 
His  experiments  were  conducted  with  100-eV  mercury  ions 
bombarding  a copper  target.  These  conditions  were  ideal  to 
simulate  with  modifications  to  the  existing  model. 

The  sputtering  simulations  for  a normally 
incident,  100-eV  ion  were  run  using  He+,  Ne+,  Ar+,  Cu+, 
Kr+,  Xe  + , Au+,  and  Hg+  ions.  Figure  14  shows  the 
sputtering  ratio  as  a function  of  the  ion  mass  for  both 
potential  functions.  As  mentioned  in  the  objectives,  a 
decrease  in  the  yield  was  evident  for  the  heavier  ions.  The 
KSE-B  curve  shows  a more  natural,  smooth  behavior  than  that 
of  Gibson-II. 


(1)  Helium  Ions . The  results  when  100-eV 
helium  icns  impact  the  target  area  are  shown  in  figures  15, 
16,  and  17.  Figure  15  shows  the  sputtering  summaries  using 

(a)  Gibson-II,  resulting  in  a sputtering  ratio  of  0.444,  and 

(b)  KSE-3,  resulting  in  a sputtering  ratio  of  0.111.  The 
predominance  of  the  sputtered  atoms  in  the  <100>  and  < 1 0 1 > 
directions  is  quite  evident.  Figure  16  shows  the 
conditional  probability  that  a ring  of  atoms,  as  a function 
of  the  ring  distance  from  the  impact  area  reference  point, 
will  contain  an  atom  that  sputtered  for  (a)  Gibson-II  and 
(b)  KSE-B.  The  term  "ring  sputtering"  probability, 
henceforth  ring  probability,  was  used  as  a label  to  mean  the 
above  conditional  probability.  The  relative  amplitudes  of 
the  ring  probabilties  could  be  misleading  because  the 
calulaticn  does  not  take  into  effect  the  number  of  atoms  in 
each  ring.  Again,  the  predominance  of  rings  that  contain 
atoms  in  the  <100>  and  < 1 0 1 > directions  is  evident.  The 
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peak  of  the  KSE-3  curve  was  closer  to  the  center  than  the 
Gibson-II  peaks,  which  were  spread  out  all  the  wa y to  the 
lattice  edges.  Figure  17  shows  the  angular  dispersion 
results  for  (a)  Gibson-II  and  (b)  KSE-3.  The  particular 
seed  used  in  the  psuedo-r andom  number  generator  to  locate 
the  isotope  random  sites  was  such  that  only  light  isotopes 
sputtered  in  this  case. 

(2)  Neon  Ions^  The  results  when  100-eV  neon 
ions  impact  the  target  area  are  shown  in  figures  18,  1 Q,  and 
20.  Figure  19  shows  that  for  the  light  mass  of  the  neon 
ion,  atoms  on  the  lattice  boundary  are  sputtered  using  (a) 
Gibson-II  and  (b)  KSE-3.  Gibson-II  resulted  in  a sputtering 
ratio  of  0.889,  while  the  sputtering  ratio  for  KSE-3  was 
0.806.  Figure  19  shows  that  the  Ne+  sputtered  atoms  were 
further  away  from  the  lattice  canter  than  the  He+,  for  both 
(a)  Gibson-II  and  (b)  KSE-B.  This  was  confirmed  by  the 
increased  number  of  next-nearest  neighbor  atoms  (in  the 
< 1 0 1 > direction)  sputtering.  Figure  20  indicates  that  some 
of  the  heavier  isotopes  sputtered  for  (a)  Gibson-II  and  (b) 
KSE-B.  There  were  not  enough  sputtered  atoms  to  notice  any 
trend  in  the  isotope  angular  dispersion. 

(3)  A rgon  lons^  The  argon  simulations  w°re 
used  to  investigate  the  "reflective  collision"  theory. 
Argon  ions  of  100-eV  energy  bombarded  the  target  area  using 
both  of  the  potential  functions.  In  the  random  isotope  site 
locations,  the  PKA  (atom  22)  was  a heavy  isotope,  wnile 
three  (52,  58,  and  67)  of  the  four  atoms  directly  underneath 
were  light  isotopes.  Four  simulations  were  run  using  each 
of  the  potential  functions;  (1)  the  random  isotope  sites, 
(2)  forcing  the  PKA  to  a light  isotope,  (3)  forcing  the  PKA 
to  a heavy  isotope  and  all  four  atoms  underneath  to  a light 
isotope,  and  (4)  forcing  the  PKA  to  a light  isotope  and  all 
four  atoms  underneath  to  a heavy  isotope.  The  remainder  of 
the  lattice  isotopes  were  unchanged. 
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Figure  21  shows  the  sputtering  summaries 
for  all  four  combinations  using  (a)  Gibson-II  and  (b)  KSE-B. 
In  all  four  cases,  identical  atoms  sputtered  independent  of 
isotope  locations.  Gibson-II  resulted  in  sputtering  ratios 
of  1.056,  while  KSE-B  resulted  in  sputtering  ratios  of 
0.944.  Figure  22  illustrates  the  ring  probabilty  in  all 
four  cases  for  (a)  Gibson-II  and  (b)  KSE-B.  Figure  23 
illustrates  the  angular  dispersion  of  ejected  atoms,  using 
Gibson-II,  for  (a)  case  (1),  (b)  case  (2),  (c)  case  (3),  and 
(d)  case  (4)  . No  light  isotopes  were  sputtered.  Figure  24 
illustrates  the  same  results  for  KSE-B. 

Since  at  least  four  atoms  were  ejected  into 
a different  five  degree  window  for  Gibson-II,  the  sputtering 
mechanisms  were  checked.  The  most  active  event  occurred 
when  the  impact  area  was  struck  by  the  ion  at  impact  point 
4060  (impact  displaced  0.40  LU  in  the  X direction  and  0.60 
LU  in  the  Z direction  from  the  center  of  the  PK& , hence 
4060)  . Four  atoms  were  sputtered  at  this  point.  A 
single-impact  program  was  run  at  this  point  with  the  output 
listing  the  positions,  velocities,  kinetic  energy,  and 
potential  energy  of  each  atom  with  a total  energy  greater 
than  0.1  eV.  The  listing  was  printed  once  every  five 
timesteps.  A timestep  is  the  time  interval  between  dynamic 
calculations.  It  is  variable  and  determined  as  the  time 
required  for  the  fastest  moving  atom  to  move  0.1  LU  from  its 
previous  position.  In  real  time  the  timestep  is  on  the 
order  of  10-15  seconds.  Figure  25  illustrates  atom 
trajectories  in  the  sputtering  of  atoms  17,  32,  57,  and  77. 
In  cases  (3)  and  (4),  the  sputtering  mechanisms  were  the 
same  displaying  only  a slight  difference  in  the  ejection 
angles  and  the  velocities  of  the  sputtered  atoms. 


¥ 


36 


The  argon  ion  collides  with  atom  22  (PKA) 
and  atom  27  by  timestep  5,  pushing  them  down  and  under  atoms 
17  and  32  respectively.  The  ion  continued  downward  and 
collided  with  atom  67  by  timestep  15.  The  ion  reflected 
back  up  through  the  hole  vacated  by  atoms  22  and  27.  By 
timestep  15,  atoms  17  and  32  were  moving  upward.  Atom  57 
moved  down  colliding  with  atoms  107  and  112  by  timestep  20. 
Atoms  107  and  112  squeezed  [21]  under  atoms  102  and  117 
respectively.  Atoms  102  and  117  were  moving  up  by  timestep 
30,  and  struck  atoms  57  and  77  from  below  by  timestep  40. 
By  this  time,  atoms  17  and  32  had  moved  far  enough  above  the 
lattice  with  sufficient  upward  velocity  to  be  classified  as 
"sputs".  Atom  17  was  ejected  at  an  angle  of  33.56°  to  the 
surface  normal  with  an  upward  velocity  of  2924  m/sec  for 
case  (3)  , and  at  an  angle  of  33.36°  with  a velocity  of  2974 
m/sec  for  case  (4) . Atom  32  was  ejected  at  an  angle  of 
33.65°  with  an  upward  velocity  of  2926  m/sec  for  case  (3), 
and  at  an  angle  of  33.56°  with  a velocity  of  2918  m/sec  for 
case  (4) . Atoms  57  and  77  continued  upward  through  the 
holes  left  by  the  sputtered  atoms  17  and  32.  As  atom  57 
passed  through  this  hole  its  trajectory  was  bent  closer  to 
the  surface  normal  by  atoms  12,  16,  and  21.  Likewise,  atom 
77  was  reflected  upward  by  atoms  28,  33,  and  37.  3v 
timestep  75,  atoms  57  and  77  were  far  enough  above  the 
lattice  surface  with  sufficient  upward  velocity  to  be 
classified  as  "sputs”.  Atom  57  was  ejected  at  an  angle  of 
18.12°  with  an  upward  velocity  of  2967  m/sec  for  case  (3), 
and  at  an  angle  of  17.75°  with  a velocity  of  2882  m/sec  for 
case  (4) . Atom  77  sputtered  at  an  angle  of  25.04°  with  an 
upward  velocity  of  2799  m/sec  for  case  (3) , and  at  an  angle 
of  24.75°  with  a velocity  of  2859  m/sec  for  case  (4). 
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(4)  Copper  Ions^  The  results  when  100-eV 
copper  ions  impact  the  target  area  are  shown  in  figures  26, 
27,  and  28.  Figure  25  shows  the  sputtering  summaries  for 

(a)  Gibson-II,  resulting  in  a sputtering  ratio  0.944,  and 

(b)  KSE-B,  resulting  in  a sputtering  ratio  of  1.111. 
Gibson-II  simulations  with  Ar+  and  Cu+  were  the  only  runs 
that  had  atoms  sputter  from  the  second  layer.  Also,  these 
two  ions  resulted  in  the  highest  sputtering  ratios  for  both 
potential  functions.  These  results  further  indicate  a 
maximum  sputtering  ratio  occurs  near  an  ion/atom  mass  ratio 
of  unity.  Figure  27  illustrates  a continued  spread  outward 
of  the  ring  probabilty  for  (a)  Gibson-II  and  (b)  KSE-3. 
Figure  23  illustrates  the  angular  dispersion  results  for  (a) 
Gibson-II  and  (b)  KSE-B. 


(5)  Krypton  Ions_.  The  results  when  100-eV 
krypton  ions  impact  the  target  area  are  shown  in  figures  29, 
30,  and  31.  Figure  29  shows  the  sputtering  summaries  for 

(a)  Gibson-II,  resulting  in  a sputtering  ratio  of  0.694,  and 

(b)  KSE-B,  resulting  in  a sputtering  ratio  of  1.139.  Here, 
a maximum  in  the  sputtering  ratio  was  reached  for  KSE-B, 
while  the  sputtering  ratio  for  Gibson-II  dropped 
drastically.  Figure  30  illustrates  a continued  outward 
spread  of  the  ring  probabilities  for  (a)  Gibson-II  and  (b) 
KSE-B.  It  is  noticed  that  the  outward  spread  still  occured, 
even  though  the  Gibson-II  yield  decreased.  Figure  31 
illustrates  the  angular  dispersion  results  for  (a)  Gibson-II 
and  (b)  KSE-B. 
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(6)  Ion s_;_  The  results  when  100-eV  xenon 
ions  impact  the  target  area  are  shown  in  figures  32,  33,  and 
34.  Figure  32  shows  the  sputtering  summaries  for  (a) 
Gibson-II,  resulting  in  a sputtering  ratio  of  0.556,  and  (b) 
KSE-B,  resulting  in  a sputtering  ratio  of  1.028.  Figure  33 
illustrates  the  ring  probabilities  for  (a)  Gibson-II  and  (b) 
KSE-B.  The  outward  spread  of  the  ring  probability  peaks  is 
still  seen  for  KSE-B,  while  Gibson-II  remained  constant. 
Figure  34  illustrates  the  angular  dispersion  results  for  (a) 
Gibson-II  and  (b)  KSE-3. 

(7)  Gold  Ions..  The  results  when  100-eV  gold 
ions  impact  the  target  area  are  shown  in  figures  35,  36,  and 
37.  Figure  35  shows  the  sputtering  summaries  for  (a) 
Gibson-II,  resulting  in  a sputtering  ratio  of  0.528,  and  (b) 
KSE-B,  resulting  in  a sputtering  ratio  of  0.778.  Figure  36 
illustrates  the  ring  probabilities  for  (a)  Gibson-II  and  (b) 
KSE-B.  Gibson-II  has  settled  out  with  approximately  equal 
probability  of  sputtering  from  atoms  in  the  first  two  rings 
in  the  <101>  direction.  In  the  case  of  KSE-3,  a snift  from 
the  first  and  second  rings  in  the  <100>  direction  toward  the 
second  ring  in  the  < 1 0 1 > direction  is  noticed.  Figure  37 
illustrates  the  angular  dispersion  results  for  (a)  Gibson-II 
and  (b)  KSE-B. 

(8)  Mercury  Ions . These  Hg+  bombardment 
results  are  compiled  from  the  sets  of  simulation  runs  used 
to  investigate  the  "reflective  collision"  theory.  Figure  38 
shows  the  sputtering  summaries  for  all  four  cases,  in  which 
the  masses  of  the  PKA  and  the  four  atoms  underneath  are 
varied.  The  results  for  (a)  Gibson-II  were  identical  in  all 
four  cases,  resulting  in  a sputtering  ratio  of  0.528.  The 
sputtering  ratio  was  0.743  using  KSE-B  in  a lattice 
consisting  of  atoms  with  only  the  natural  mass.  The  results 
for  (b)  KSE-9  depended  on  the  order  of  the  heavy  over  light 
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interaction.  When  the  PKA  was  a heavy  isotope,  shown  in 
figure  38(b),  the  sputtering  ratio  was  0.7229  However,  when 
the  PKA  was  a light  isotope,  atom  30  sputtered  an  additional 
time  (denoted  by  ▲ ) and  the  sputtering  ratio  became  0.750. 
Figure  39  illustrates  the  ring  probabilties  for  (a) 
Gibson-II,  which  remained  relatively  unchanged,  and  (b) 
KSE-B,  which  still  displayed  a shift  to  the  second  ring  in 
the  101  direction.  Figure  40  illustrates  the  angular 
dispersion  results  for  Gibson-II.  Figure  41  shows  the 
angular  dispersion  results  for  K5E-3.  Again,  atom  30 
sputtered  an  additional  time.  The  additional  ejection  angle 
is  denoted  by  ▲ . 

b.  1-keV  Results 

The  1-keV  events  were  obviously  not  contained  in 
the  9x4x9  aicrocrystallite . However,  the  simulations  were 
run  to  compare  the  sputtering  ratios  with  experimental 
results.  Also,  the  dependence  of  the  sputtering  on  ion  mass 
was  investigated.  Useful  information  can  be  obtained  about 
the  ring  probabilties  and  isotope  angular  dipersion  at  this 
higher  energy.  The  162  atom,  9x4x9  microcrysrallite,  with 
the  correct  ratio  of  copper  isotopes  at  random  sites,  was 
used.  Even  though  all  atoms  were  used  in  the  calculations 
and  the  printout  provided  information  or.  all  atoms,  only  the 
top  layer  results  are  shown  in  the  figures.  Even  in  the 
worse  case  (Cu/Ar  + , Gibson-II),  more  than  86  per  cent  of 
the  sputtered  atoms  were  from  the  top  layer. 

The  sputtering  simulations  for  normally 
incident,  1-keV  ions  were  run  using  He+,  Ne+,  Ar*,  Cu+, 
Kr+,  Xe+,  Au+,  and  Hg+  ions.  Figure  42  shows  the  sputtering 
ratio  as  a function  of  the  ion  mass  for  both  potential 
functions.  As  expected,  the  sputtering  ratios  have 
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increased.  However,  they  are  no  longer  smooth  and 
continuous,  but  exhibit  a randomness,  as  if  each  was 
independent  of  the  others. 

All  the  1-keV  simulations  used  the  same  random 
isotope  sites.  The  letters  L and  H,  denoting  light  or  heavy 
isotopes,  appear  only  in  the  Cu/Ne+,  KSE-B  summary  (more 
different  atoms  sputtered  - but  not  the  highest  sputtering 
ratio)  . 

(1)  Hsliua  Ions . The  results  when  1-kaV  helium 
ions  impact  the  target  area  are  shown  in  figures  43,  44,  and 
45.  Figure  43  shows  the  sputtering  summaries  for  (a) 
Gibson-II,  resulting  in  a sputtering  ratio  of  1.914,  and  (b) 
KSE-B,  resulting  in  a sputtering  ratio  of  2.833.  In  both 
cases,  boundary  atoms  sputtered.  Figure  44  illustrates  the 
ring  probabilities  for  (a)  Gibson-II  and  (b)  KSE-B.  3oth 
display  a large  action-reaction  effect,  as  if  a bomb  were 
dropped  at  the  center  and  everything  out  from  the  center  was 
blasted  upward.  Figure  45  illustrates  the  angular 
dispersion  results  for  (a)  Gibson-II  and  (b)  KSE-B. 

(2)  Neon  lons^  The  results  when  1-keV  neon  ions 
impact  the  target  area  are  shown  in  figures  46,  47,  and  48. 
Figure  46  shows  the  sputtering  summaries  for  (a)  Gibson-II, 
resulting  in  a sputtering  ratio  of  3.686,  and  (b)  KSE-3, 
resulting  in  a sputtering  ratio  of  4.500.  Host  of  the 
surface  atoms  were  sputtered,  indicating  that  the  energy  ( 
or  momentum  ) spread  out  over  the  entire  lattice  surface. 
Figure  47  illustrates  the  ring  pr ooabilities  for  (a) 
Gibson-II  and  (b)  KSE-B.  The  KSE-3  simulation  indicates  a 
spreading  outward  of  the  sputtered  atoms.  The  Gibson-II 
simulation  not  only  exhibits  increased  sputtering  of  atoms 
further  from  the  center,  but  it  also  clearly  shows  a shift 
in  the  number  of  sputtered  atoms  from  those  at  2.83  LU  in 
the  <131>  direction  to  those  at  3.16  L'J  in  the  < 1 0 3 > 
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direction.  As  mentioned  earlier,  this  could  be  artificial 
because  the  ring  at  3.16  LU  is  the  first  ring  to  contain 
eight  atoms.  Figure  48  illustrates  the  angular  dispersion 
results  for  (a)  Gibson-II  and  (b)  KSE-B. 

(3)  Argon  Ions^  The  results  when  1-keV  argon 
ions  impact  the  target  area  are  shown  in  figures  49,  50,  and 
51.  Figure  49  shows  the  sputtering  summaries  for  (a) 
Gibson-II,  resulting  in  a sputtering  ratio  of  3.806,  and  (b) 
KSE-B,  resulting  in  a sputtering  ratio  of  4.611.  Figure  50 
illustrates  the  ring  probabilities  for  (a)  Gibson-II  and  (b) 
KSE-B.  The  Gibson-II  simulation  continued  to  increase  in 
the  outward  direction  at  the  expense  of  the  atoms  at  3.16  LO 
in  the  <103>  direction.  However,  the  KSE-3  simulation 
reversed  its  trend  and  the  ring  probability  for  the  close-in 
rings  (1.41  LU  and  2.00  LU)  increased  at  the  expense  of  the 
ring  at  2.83  LU  and  the  rings  further  out  from  the  center. 
Figure  51  illustrates  the  angular  dispersion  results  for  (a) 
Gibson-II  and  KSE-B. 

(4)  Copper  lons^  The  results  when  1-keV  copper 
ions  impact  the  target  area  are  shown  in  figures  52,  53,  and 
54.  Figure  52  shows  the  sputtering  summaries  for  (a) 
Gibson-II,  resulting  in  a sputtering  ratio  of  3.029,  and  (b) 
KSE-B,  resulting  in  a sputtering  ratio  of  3.972.  These  low 
sputtering  ratios  not  only  disagree  with  a monatonic 
increase  in  the  yield  with  mass,  but  are  drastically  low  for 
the  maximum  to  occur  at  a ion/atom  mass  ratio  of  unity. 
Figure  53  illustrates  the  ring  probabilities  for  (a) 
Gibson-II  and  (b)  KSE-B.  The  Gibson-II  simulation  showed  an 
ubexpected  change.  The  ring  probability  increased  for  the 
rings  at  2.00  LU  (<100>  direction)  and  3.16  LU  (<  1 0 3 > 
direction) . This  trend  was  also  displayed  in  the  KSE-B 
simulation,  where  it  starts  to  occur  in  the  Cu/Ar* 
simulation.  Figure  54  illustrates  the  angular  dispersion 
results  for  (a)  Gibson-II  and  (b)  KSE-3. 
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(5)  Krypton  Ions . The  results  when  1-keV 
krypton  ions  impact  the  target  area  are  shown  in  figures  55, 
56,  and  57.  Figure  55  shows  the  sputtering  summaries  for 

(a)  Gibson-II,  resulting  in  a sputtering  ratio  of  2.833,  and 

(b)  KSE-B,  resulting  in  a sputtering  ratio  of  3.917.  note 
that  the  sputtering  ratio,  for  both  potential  functions, 
continues  to  decrease  with  increased  ion  mass.  Figure  56 
illustrates  the  ring  probabilities  for  (a)  Gibson-II  and  (b) 
KSE-B.  Both  potential  function  simulations  continue  to 
exhibit  the  trend  to  establish  a peak  at  2.00  LU  in  the 
<100>  direction  and  at  3.16  LU  in  the  <103>  direction. 
However,  the  KSE-3  simulation  shows  significantly  larger 
ring  probability  at  3.16  LU.  Figure  57  illustrates  the 
angular  dispersion  results  for  (a)  Gibson-II  and  (b)  KSE-3. 


(6)  Xenon  Ions4  The  results  when  1-keV  xenon 
ions  impact  the  target  area  are  shown  in  figures  58,  59,  and 
60.  Figure  58  shows  the  sputtering  summaries  for  (a) 
Gibson-II,  resulting  in  a sputtering  ratio  of  2.333,  and  (b) 
KSE-B,  resulting  in  a sputtering  ratio  of  3.472.  Both 
potential  functions  continued  to  show  a decrease  in  the 
sputtering  ratio  with  an  increase  in  the  mass  of  the  ion. 
Figure  59  illustrates  the  ring  probabilties  for  (a) 
Gibson-II  and  (b)  KSE-B.  All  the  previous  trends,  except  an 
increased  ring  probability  at  3.16  LU  were  shattered. 
Figure  60  illustrates  the  angular  dispersion  results  for  (a) 
Gibson-II  and  (b)  KSE-B. 

(7)  Gold  Ions . The  results  when  1-keV  gold 
ions  impact  the  target  area  are  shown  in  figures  61,  62,  and 
63.  Figure  61  shows  the  sputtering  summaries  for  (a) 
Gibson-II,  resulting  in  a sputtering  ratio  of  2.178,  and  (b) 
KSE-3,  resulting  in  a sputtering  ratio  of  4.167.  A definite 
increase  in  the  sputtering  ratio  occured.  Figure  62 
illustrates  the  ring  probabilities  for  (a)  Gibson-II  and  (b) 
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KSE-B.  So  new  trends  were  noted,  except  that  the  ring  at 
3.16  L U continues  to  dominate  and  the  other  ring 
probabilities  fluctuate.  Figure  63  illustrates  the  angular 
dispersion  results  for  (a)  Gibson-II  and  (b)  KSE-3. 

* (8)  Mercury  I ons^  The  results  when  1-keV 
mercury  ions  impact  the  target  area  are  shown  in  figures  64, 
65,  and  66.  Figure  64  shows  the  sputtering  summaries  for 

(a)  Gibson-II,  resulting  in  a sputtering  ratio  of  2.778,  and 

(b)  KSE-B,  resulting  in  a sputtering  ratio  of  4.278.  These 
simulations  resulted  in  no  change  in  the  sputtering  ratio 
for  Gibson-II,  and  only  a slight  increase  for  KSE-B.  Figure 
65  illustrates  the  ring  probabilities  for  (a)  Gibson-II  and 
(b)  KSE-B.  Only  insignificant  changes  occured  compared  to 
the  gold  ion  results.  Figure  66  illustrates  the  angular 
dispersion  results  for  (a)  Gibson-II  and  (b)  KSE-B. 

c.  2-keV  Results 


The  2-keV  ion  simulations  were  the  earliest  runs 
made.  No  isotope  information  was  obtained.  As  a 
consequence,  the  only  results  are  the  sputtering  ratio 
versus  ion  mass  curve,  figure  67,  for  both  potential 
functions.  They  exhibit  an  even  more  erratic  behavior  than 
did  the  1-keV  simulations.  The  sputtering  summaries 
frequently  show  atoms  sputtered  from  the  third  layer,  and 
rarely  an  atom  sputtered  from  the  fourth  layer.  At  certain 
impact  points  as  many  as  25  atoms  sputtered.  The  small, 
9x4x9  microcrystallite  was  effectively  blasted  apart  by  the 
2-keV  ions. 


2 • .HI  CL  0 r ie a t ajt  1.0 n 1.0  0 eV 

The  effort  devoted  to  (110)  orientation  simulation 
runs  was  very  limited  due  to  the  higher  priority  assigned  to 
the  (100)  orientation  runs.  is  a result,  only  the  argon  and 
mercury  ion  simulations  at  100  eV  were  run.  Still  needing 
investigation  are  the  sputtering  ratio  as  a function  of  ion 
mass  and  sputtering  mechanism  traces. 


The  (110)  orientation  microcrystallite  was 
constructed  in  a lattice  consisting  of  123  copper  isotopes. 
In  the  random  isotope  site  locations,  the  PKA  (14)  was  a 
light  isotope,  as  were  three  (33,  36,  and  37)  of  the  four 
atoms  underneath  the  PKA.  Figure  4 shows  an  isometric 
projection  of  the  lattice  and  its  numbering  scheme.  Figure 
10  shows  an  expanded  view  of  the  grid  of  36  impact  points 
used.  A Sx6x9  lattice  was  used  with  no  increase  in  computer 
time  required  to  complete  a run. 


Two  sets  of  simulations  were  run  to  investigate  tne 
"reflective  collision"  theory.  Four  runs  were  made  in  each 
set  of  simulations  for  each  potential  function.  The  impact 
area  was  bombarded  by  100-eV,  araon  and  mercury  ions.  The 
four  runs  consisted  of;  (1)  the  random  isotope  sites  (PKA 
was  a light  isotope) , (2)  forcing  the  PKA  to  a heavy 
isotope,  (3)  forcing  the  PKA  to  a light  isotope  and  all  four 
atoms  underneath  to  a heavy  isotope,  and  (4)  forcing  the  PKA 
to  a heavy  isotope  and  all  four  atoms  underneath  to  a light 
isotope.  The  remainder  of  the  lattice  isotopes  were 
unchanged. 


a.  Argon  Ions 

The  results  when  100-eV  argon  ions  impact  the 
target  area  are  shown  in  figures  68-75.  Figure  63  shows  the 
sputtering  summaries  using  Gibson-II  for  all  four  cases. 
Case  (1)  and  case  (3)  , shown  in  figure  63  (a)  , resulted  in  a 
sputtering  ratio  of  0.889.  Case  (2)  and  case  (4),  shown  in 
figure  63(b),  resulted  in  a sputtering  ratio  of  0.361.  A 
predominance  of  sputtered  atoms  in  the  <100>  and  <101> 
directions  is  evident.  Figure  69  illustrates  the  ring 

probabilities  for  (a)  case  (1)  and  case  (3),  and  (b)  case 
(2)  and  case  (4)  using  Gibson-II.  Notice  that  an  atom  at 
the  boundary  ( in  the  <101 > direction  ) sputtered.  Figure 
70  illustrates  the  angular  dispersion  results  for  (a)  case 
(1) , and  (b)  case  (3)  using  Gibson-II.  Figure  71  shows  the 
same  information  for  (a)  case  (2)  , and  (b)  case  (4)  . 

Figure  72  shows  the  sputtering  summaries  using 
KSE-3  for  all  four  cases.  Case  (3)  resulted  in  a sputtering 
ratio  of  0.750,  while  case  (1)  sputtered  atom  26  one 
additional  time  (denoted  by  ▲ ) and  had  a sputtering  ratio 
of  0.773.  Case  (1)  and  case  (3)  are  shown  in  figure  72(a)  . 
Case  (2)  resulted  in  a sputtering  ratio  of  0.778,  while  case 
(4)  sputtered  atom  3 one  additional  time  and  had  a 
sputtering  ratio  of  0.305.  Case  (2)  and  case  (4)  are  shown 
in  figure  72(b).  Again,  it  is  evident  that  only  atoms  in 
the  <100>  and  < 1 1 0 > directions  sputtered.  Figure  73 
illustrates  the  ring  probabilities,  using  FSS-9,  for  (a) 
case  (1)  and  case  (3)  , and  (b)  case  (2)  and  case  (4)  . The 

close-in  trend,  similar  to  that  of  Gibson-II,  is  evident,  as 
are  the  sputtered  atoms  at  the  boundary  in  the  <101> 
direction.  Figure  74  illustrates  the  angular  dispersion 
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results,  using  KSE-B,  for  (a)  case  (1)  , and  (b)  case  (3). 
Figure  75  shows  the  same  information  for  (a)  case  (2) , and 
(b)  case  (4)  . 

b.  Mercury  Ions 

The  results  when  100-eV  mercury  ions  impact  the 
target  area  are  shown  in  figure  76,  77,  and  78.  Figure 
76(a)  shows  the  sputtering  summaries  for  all  four  cases, 
using  Gibson-II.  Case  (1),  case  (2),  and  case  (4)  resulted 
in  a sputtering  ratio  of  0.278.  Case  (3)  sputtered  atom  13 
one  additional  time  (denoted  by  ▲ ) and  had  a sputtering 
ratio  of  3.05.  Figure  76(b)  shows  the  sputtering  summaries 
for  all  four  cases,  using  KSE-B,  which  resulted  in 
sputtering  ratios  of  0.222.  Figure  77  illustrates  the  ring 
probabilities  for  (a)  Gibson-II  and  (b)  KSE-3.  Figure  78 
illustrates  the  angular  dispersion  results  for  (a)  Gibson-II 
and  (b)  KSE-B. 

No  significant  trends  were  noted  in  the  (110) 
orientation.  However,  the  ease  with  which  atoms  sputtered 
in  the  <101>  direction  is  evident. 

3 . (Ill)  Orientation  10  Oe  V 

Like  the  (110)  orientation,  the  (111)  orientation 
did  not  receive  extensive  attention.  Both  the  sputtering 
ratio  as  a function  of  the  ion  mass,  and  sputtering 
mechanism  traces  need  further  investigation. 

The  (111)  orientation  microcrystallite  is 
constructed  in  a 9x4x9  lattice  consisting  of  163  atoms.  In 
the  random  isotope  site  locations,  the  PKA  (22)  was  a light 
isotope,  as  were  three  (58,  63,  and  67)  of  the  four  atoms 
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underneath  the  PKA.  Figure  5 shows  an  isometric  projection 
of  the  lattice  and  its  numbering  scheme.  Figure  11  shows  an 
expanded  view  of  the  grid  of  21  impact  points  used. 

Two  sets  of  simulations  were  run  to  investigate  the 
"reflective  collision"  theory.  Four  runs  were  made  in  each 
set  of  simulations  for  each  potential  function.  The  target 
was  bombarded  by  100-eV  argon  and  mercury  ions.  The  four 
runs  consisted  of;  (1)  the  random  isotope  sites  (PKA  was  a 
light  isotope)  , (2)  forcing  the  PKA  to  a heavy  isotope,  (3) 
forcing  the  PKA  to  a light  isotope  and  all  four  atoms 
underneath  to  a heavy  isotope,  and  (4)  forcing  the  PKA  to  a 
heavy  isotope  and  all  four  atoms  underneath  to  a light 
isotope.  The  remainder  of  the  lattice  isotopes  were 
unchanged. 

a.  Argon  Ions 


The  results  when  100-eV  argon  ions  impact  the 
target  area  are  shown  in  figures  79-83.  Figure  79  shows  the 
sputtering  summaries  for  (a)  Gibson-II,  resulting  in 
sputtering  ratios  of  0.857,  and  (b)  KSE-B,  resulting  in 

sputtering  ratios  of  0.762  . For  KSE-3,  case  (2)  , atom  16 
sputtered  one  additional  time  (denoted  by  ▲),  resulting  in 
a sputtering  ratio  of  0.310.  Figure  80  illustrates  the  ring 
probabilities  for  (a)  Gibson-II  and  (b)  KSE-B.  It  was 
unexpected  that,  for  KSE-B,  all  of  the  sputtered  atoms  are 
in  the  < 1 0 1 > direction  in  the  ring  at  2.33  LU.  Figure  31 
illustrates  the  Gibson-II  angular  dispersion  results  for  (a) 
case  (1)  and  case  (3)  , and  (b)  case  (2)  and  case  (4)  . 
Figure  82  illustrates  the  KSE-B  angular  dispersion  results 
for  (a)  case  (1)  and  (b)  case  (3)  . Figure  83  shows  the  same 
results  for  (a)  case  (2)  , and  (b)  case  (4)  . 
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b. 


Mercury  Ions 


The  results  when  100-eV  mercury  ions  impact  the 
target  area  are  shown  in  figures  84-89.  Figure  84(a)  shows 
the  sputtering  summaries,  using  Gibson-II,  for  case  (1), 
case  (2)  , and  case  (4)  , which  resulted  in  sputtering  ratios 
of  0.476.  Case  (3)  had  a sputtering  ratio  of  0.381  and  is 
shown  in  figure  84(b) . Figure  85  shows  the  sputtering 
summaries  using  the  KS2-B.  Case  (3)  resulted  in  a 
sputtering  ratio  of  0.762,  while  case  (1)  sputtered  two 
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87  illustrates  the  Gibson-II  angular  dispersion  results  for 

(a)  case  (1)  , case  (2)  , and  case  (4)  , and  (b)  case  (3)  . 

Figure  83  illustrates  the  KSE-B  angular  dispersion  results 

for  (a)  case  (1)  and  (b)  case  (3)  . Figure  89  shows  the  same 

/ 

information  for  (a)  case  (2)  , and  (b)  case  (4)  . 
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mercury  ions  bombarding  a 13x6x13,  (110)  orientation 
microcrystallite  containing  255  atoms  are  shown  in  figure 
90.  Figure  90(a)  is  for  Gibson-II,  while  figure  90(b)  is 
for  KSE-B.  Comparison  of  figure  90  with  figure  76, 
confirmed  that  the  9x5x9  microcrystallite  did  not  contain 
the  100-eV,  mercury  event.  Based  on  the  original  definition 
of  containment,  the  13x6x13  microcrystallite  did  not  contain 
the  Gibson-II  simulation  run. 

In  the  (HI)  orientation,  the  sputtering  summaries  for 
mercury  ions  bombarding  a 13x4x13  microcrystallite 
containing  339  atoms  are  shown  in  figure  91.  Figure  91(a) 
is  for  Gibson-II,  while  figure  91  (b)  is  for  KSE-B. 
Comparison  of  figure  91  with  figures  84  and  85,  indicate 
that  even  though  the  9x4x9  lattice  had  no  boundary  atoms 
sputtered,  the  sputtering  ratio  increases  when  a 13x4x13 
lattice  is  used. 

1 • 1 00-eV  Containment 

Based  on  the  above  information,  the  definition  for 
"containment1*  was  changed  to:  "an  event  is  contained  when  an 
increase  in  the  microcrystallite  size  does  not  result  in  any 
further  increase  in  the  yield."  Earlier  simulations 
resulted  in  larger  yields  for  KSE-B.  For  this  reason,  KSE-B 
was  used  to  investigate  the  size  microcrystallite  needed  to 
contain  a 100-eV  event. 

The  sputtering  ratio  as  a function  of  the  ion  mass 
was  examined  for  the  (100)  orientation.  Figure  92  shows  the 
results  for  9x4x9,  11x4x11,  and  11x6x11  microcrystallites. 
Additionally,  simulations  were  run  for  copper  (largest 
sputtering  ratio)  and  mercury  ions  for  13x4x13,  13x6x13,  and 
17x4x17  microcrystallites.  Based  on  the  new  definition  of 
containment,  100-eV  helium  and  xenon  events  are  contained  in 
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a 9x4x9  microcrystallite  because  an  11x4x11  microcr ystallite 
produces  the  same  yield.  However,  the  argon,  copper,  and 
krypton  events  were  not  contained.  Further  increase  of  the 
mi'crocrystallite  size  demonstrated  that  the  100-eV,  copper 
event  is  contained  by  a 13x4x13  microcrystallite.  Figure  92 
also  illustrates  that  the  maximum  sputtering  ratio  occurs 
when  the  ion/atom  mass  ratio  is  approximately  unity.  Figure 
93  shows  the  same  information  in  a plot  of  sputtering  ratio 
versus  lattice  size  for  all  eight  ions. 

2.  1 ~xi  x 1 7 Mi crocrysta Hite 

Asuming  that  all  100-eV  events  are  contained  by  a 
17x4x17  microcrystallite,  the  effect  of  increasing  the  ion 
energy  was  investigated.  The  reason  for  this  investigation 
was  to  determine  how  energetic  an  event  could  be  contained 
by  a 17x4x17  microcrystallite.  And  to  compare  higher  energy 
events  with  experimental  results.  The  reason  for  not 
increasing  the  lattice  size  without  limit  was  obviously  to 
keep  the  computer  run  time  within  reason.  Again,  only  KSE-B 
with  mercury  and  copper  ions  were  investigated. 

a.  Copper  Ions 

The  results  when  100-eV  copper  ions  impact  the 
target  area  are  shown  in  figures  94,  95,  and  96.  Figure  94 
shows  the  sputtering  summary  with  a sputtering  ratio  of 
1.389.  Containment  for  the  100-eV  copper  event  had  already 
been  established.  Figure  95  illustrates  the  "atom 
sputtering  probability",  henceforth  called  the 
"probaoilit y".  The  "probability"  is  the  "ring  sputtering 
probability"  divided  by  the  number  of  atois  in  each  ring. 
The  dimensions  for  "probability"  is  sputs  per  atom  divided 
by  the  total  sputs  at  a fixed  distance  from  the  impact  area 
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reference  point.  Figure  95  indicates  an  enhancement  of 
sputtered  atoms  in  the  <100>  and  < 1 0 1 > directions.  Figure  96 
illustrates  the  angular  dispersion  results. 

Results  for  the  200-eV  copper  ion  simulations 
are  shown  in  figures  97,  98,  and  99.  Figure  97,  the 
sputtering  summary  with  a sputtering  ratio  of  2.889,  shows 
that  boundary  atoms  are  sputtered,  and  that  the  event 
definitely  is  not  contained.  Figure  98  illustrates  the 
"probability"  results  and  indicates  that  an  atom  farther  out 
had  increased  chances  of  sputtering.  Figure  99  illustrates 
the  angular  dispersion  results.  A bias  toward  light 
isoropes  closer  to  the  normal  is  evident. 

The  results  for  the  500-eV  copper  ion 
simulations,  with  a sputtering  ratio  of  4.028,  are  shown  in 
figures  100,  101,  and  102.  Again,  boundary  atoms  are 
sputtered  and  a slight  shift  outward  of  the  "probability"  is 
noticed.  Figure  102  illustrates  rhe  angular  dispersion 
results  which  showed  that  more  light  isotopes  are  ejected 
more  nearly  normal  to  the  lattice  surface. 
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Obviously,  the  200-eV,  500-eV,  and  1-keV  copper 
events  are  not  contained,  even  in  a 17x4x17 
micro crystallite . 
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b.  Mercury  Ions 

The  results  for  the  100-eV  mercury  ion 
simulations  are  shown  in  figures  106,  107,  and  108.  Again, 
the  11x4x11  lattice  has  already  been  determined  to  contain 
the  100-eV,  mercury  event.  Figure  106  shows  the  sputtering 
summary  with  a sputtering  ratio  of  0.833.  Figure  107 
illustrates  the  "probability",  and  like  the  copper 
simulation,  was  predominant  in  the  < 1 0 1 > direction.  Figure 
108  illustrates  the  angular  dispersion  results. 

The  200-eV  mercury  ion  results,  with  a 
sputtering  ratio  of  1.722,  are  shown  in  figures  109,  110, 
and  111.  Figure  109,  the  sputtering  summary,  shows  an 
outward  spreading  of  the  sputtered  atoms.  However,  no 
boundary  atoms  were  sputtered  and  the  200-eV,  mercury  event 
is  tentatively  classified  as  "contained".  Figure  1 TO 
confirms  the  outward  shift  in  "probability".  Figure  111 
illustrates  the  angular  dispersion  results,  which  indicate  a 
greater  number  of  heavy  isotopes  sputtering  than  the  isotope 
ratio  would  justify. 

The  500-eV  mercury  ion  results,  with  a 
sputtering  ratio  of  3.667,  are  shown  in  figures  112,  113, 
and  114.  Figure  112,  the  sputtering  summary,  indicates  that 
boundary  atoms  are  sputtered.  Figure  113  illustrates  the 
continued  outward  spreading  of  sputtered  atoms.  Figure  114 
illustrates  the  angular  dispersion  results,  and  the  number 
of  heavy  to  light  isotopes  closer  to  the  correct  ratio. 


The  1-keV  mercury  ion  results,  with  a sputtering 
ratio  of  6.139,  are  shown  in  figures  115,  115,  and  117. 
Figure  115  shows  that  boundary  atoms  are  sputtered.  Figure 

116  shows  the  "probability"  still  spreading  outward.  Figure 

117  illustrates  the  angular  dispersion  results. 


The  "probability"  figures  for  both  the  copper 
and  mercury  ions  are  remarkably  similar  and  display  the  same 
trends. 


c.  Sputtering  Mechanism  Traces 


The  containment  studies  indicated  that,  as  t he 
lattice  size  was  increased,  the  sputtering  ratio  increased. 
Investigation  into  the  similarities  of  the  17x4x17  and  9x4x9 
microcrystallites,  indicated  that  a*om  108  (17x4x17) 
sputtered  twice  in  36  events.  This  corresponds  to  atom  40 
in  the  9x4x9  microcrystallite,  which  did  not  sputter  in  36 
events.  Sputtering  mechanism  trace  runs  were  performed  to 
determine  the  sputtering  mechanism  involved  in  the 
sputtering  of  atom  108.  Atom  108  sputtered  when  ions 
impacted  at  points  50  ( impact  point  displaced  0.50  LU  in 
the  Z direction  only)  and  60  ( impact  point  displaced  0.6  1U 
in  the  Z direction  only)  . 

Figure  118  illustrates  the  a*om  trajectories  ir 
the  sputtering  of  atom  108  from  impact  point  50  ( the 
mechanism  is  the  same  for  impact  point  60) . The  copper  ion 
struck  atoms  74  (PKA)  and  91  forcing  them  down  by  timestep 
15.  As  atom  91  was  moving  downward,  atoms  9Q,  100,  and  108 
prevented  atom  91  flying  off  in  the  positive  Z iirection, 
and  channeled  it  downward.  By  timestep  25,  atom  91  had 
collided  with  atom  224,  reversing  its  direction.  At 
timestep  30,  atom  91  was  forcing  atom  108  upwari  from 
underneath.  Atom  91  transfered  almost  all  of  its  energy 
(momentum)  to  atom  108  and  had  less  than  0.1  eV  by  timestep 
45.  As  atom  108  was  going  upward,  it  also  had  a velocity  in 
the  positive  Z direction  of  greater  than  6000  a/sec.  Atoms 
116,  117,  and  125  bent  atom  108  upward  quite  radically.  The 
velocity  in  the  Z direction  had  been  reduced  to  less  “han 
700  m/sec  by  timestep  50.  3y  timestep  60,  aton  108  was  far 
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enough  above  the  lattice  surface,  with  sufficient  velocity 
to  be  classified  as  a sput.  Atom  108  was  ejected  at  an 
angle  of  19.23°  to  the  surface  normal  with  an  upward 
velocity  of  3206  rn/sec. 

Quite  evident  is  the  role  atoms  116,  117,  and 
125  have  in  channeling  atom  108  upward.  In  the  9x4x9 
lattice,  atom  40  ( a boundary  atom)  corresponds  to  atom  103. 
And  nothing  equivalent  exists  for  atoms  116,  117,  and  125. 
Investigation  was  then  made  to  determine  whether  an  atom 
corresponding  to  atom  103  could  be  sputtered  in  the  9x4x9 
lattice  by  using  the  same  impact  point  on  a different  PKA. 
In  this  simulation,  atom  31,  corresponding  to  atom  103  in 
the  larger  lattice,  was  expected  to  sputter.  Figure  119 
illustrates  the  atom  trajectories  when  the  copper  ion  struck 
the  new  PKA  (atom  13)  at  impact  point  50.  Atom  31  was 
ejected  at  an  angle  of  20.17°  with  an  upward  velocity  of 
3127  m/sec.  The  mechanism  is  the  same,  and  the  importance  of 
boundary  atoms  is  very  evident. 

Additional  runs  were  made  on  the  9x4x9  lattice 
bombarding  atom  13  at  impact  points  50  and  60.  At  both 
points,  comparison  is  made  between  runs  with  (1)  only  atoms 
of  the  natural  mass,  , and  natural  mass,  and  (2)  forcing 
atom  22  (91)  to  a heavy  isotope  and  atom  67  (224)  to  a light 
isotope.  In  all  four  cases,  the  sputtering  mechanism  is  the 
same.  The  differences  in  the  ejection  angles  and  upward 
velocities,  for  each  impact  point,  are  less  than  0.4°  and  75 
m/sec  respectively. 


3.  Ion/Atom  Mass  Ratio 

Sigmund  [28]  proposed  a theory  that  the  maximum 
sputtering  yield  increases  as  a function  of  energy. 
However,  results  of  this  investigation  strongly  indicate  a 
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dependence  on  the  ratio  of  the  ion  mass  to  the  target  atom 
mass.  A rearrangement  of  Sigmund's  theoretical  curves 
(figure  16b,  ref.  29)  to  show  the  sputtering  ratio  as  a 
function  of  the  ion/atom  mass  ratio  results  in  figure  129. 
The  curves  shown  by  Sigmund  are  predicted  for  209-eV, 
600-eV,  1-keV,  and  5-keV  ions  bombarding  an  aluminum  target. 
Figure  120  illustrates  that  Sigmund's  theory  would  predict 
the  sputtering  ratio  to  increase  as  a function  of  the 
ion/atom  mass  ratio.  Figure  121  shows  experimental  results 
obtained  by  Wehner  and  is  presented  b y Sigmund  ( figure  12, 
ref.  28  ) to  support  his  theory.  This  figure  has  beer 
replotted  to  show  the  sputtering  ratio  as  a function  of  the 
ion/atom  mass  ratio.  The  results  were  obtained  by 
bombarding  a copper  target  with  100-eV,  200-eV,  and  600-eV 
ions.  Experimental  results  plotted  in  this  manner  show  a 
decrease  in  the  sputtering  ratio  with  an  ion/atom  mass  ratio 
greater  than  unity,  in  disagreement  with  the  theory  proposed 
by  Sigmund.  Results  obtained  in  the  present  investigation 
completely  agree  with  the  shape  of  experimental  data  shown 
in  figure  121. 

C.  RANDOM  IMPACT  POINTS 


In  response  to  other  investigators'  suggestion  that  the 
target  impact  area  be  bombarded  randomly  rather  than  the 
systematic  grid  of  impact  points,  the  computer  program  was 
slightly  modified.  Since  additional  data  were  desired  for 
mercury  events,  Hg+  was  used  as  the  incident  ion.  The 
KSE-E  potential  function  was  used.  Earlier  results  revealed 
that  the  ICO-eV,  mercury  ion  event  was  contained  in  an 
11x4x11  microcrystallite. 


56 


Figure  122  illustrates  the  (100)  orientation  impact  area 
that  was  struck  at  the  144  randomly  selected  impact  points. 
A few  gaps  are  evident,  bur  the  randomness  is  very  apparent. 
Figures  123-128  show  a comparison  of  the  results  from  a 
randomly  bombarded  impact  area  and  the  results  from  the  grid 
of  impact  points.  Since  the  number  of  events  were  not  the 
same,  the  sputtering  summaries  indicate  a sputtering  ratio 
for  each  atom.  Figure  123  shows  the  normalized  sputtering 
summary  for  the  randomly  selected  impact  points.  Figure  124 
shows  the  same  information  using  the  grid  of  impact  points. 
The  four  groups  of  36  randomly  selected  impact  points  had  a 
mean  sputtering  ratio  of  0.812  with  a variance  of  0.0043. 
Additional  runs  ( not  shown  ),  combined  with  the  four 
previous  groups  to  give  a total  of  21  groups  of  36  randomly 
selected  impact  points  ( 756  shots  ) , resulted  in  a mean 
sputtering  ratio  of  0.782  with  a variance  of  0.016.  This 
compares  very  favorably  with  a sputtering  ratio  of  0.778 
using  the  grid  of  impact  points.  Figures  125  and  126  show 
the  "probability"  results  for  the  randomly  selected  impact 
points  and  the  grid  of  impact  points  respectively.  It  is 
guite  evident  that  they  are  almost  identical.  Figure  127 
illustrates  the  angular  dispersion  results  for  randomly 
selected  impact  points.  A definite  shift  of  the  lighter 
isotopes  toward  the  surface  normal  is  observed.  Figure  128 
illustrates  the  same  information  for  the  grid  of  impact 
points.  The  absence  of  the  shift  of  the  light  isotopes 
toward  the  normal  is  noted. 
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V. 


CONCLUSIONS  AND  33C0MM ENDATIONS 


A.  SIMULATION  MODEL 

The  more  sophisticated  a simulation  model  becomes,  the 
more  uncertain  are  its  results.  Some  of  these  uncertainties 
are  now  beginning  to  appear  in  this  model. 

The  lattice  generators  can  now  construct  randomly 
positioned  isotopes.  Symmetrical  unfolding  of  the  event 
results  is  the  reverse  of  the  folding  required  to  reduce  the 
representative  impact  area.  The  model  can  now  display  the 
ejection  angle  of  the  sputtered  isotopes  in  a histogram  for 
visual  comparison.  It  now  determines  the  distance  from  the 
surface  layer  center  to  each  ring  of  sputtered  atoms,  the 
number  of  atoms  and  the  number  of  times  each  atom  sputters 
in  each  ring,  and  the  probability  that  that  each  ring 
contained  an  atom  that  was  sputtered.  Other  modifications 
were  made  in  the  model  to  approximate  more  nearly  a 
meaningful  physical  decription.  These  include  changes  in 
the  distances  from  the  lattice  surface  at  which  calculations 
are  truncated,  the  depth  of  lattice  needed  to  prevent  any 
change  in  the  yield,  and  a systematic  determination  of  the 
distance  above  the  surface  at  which  an  atom  is  declared  to 
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CONTAINMENT 


The  containment  of  an  event  is  paramount  in  a complete 
description  of  the  sputtering  process.  dowever,  it  is  less 
significant  when  studying  mechanisms  and  trends. 

3ased  on  the  definition  of  containment  ( no  increase  in 
yield  with  an  increase  in  the  microcrystallite  size  ) , tne 
100-eV  results  showed  that  a 9x4x9,  (100)  orientation 
microcrystallite  contained  only  the  helium  and  xenon  events. 
The  100-eV,  mercury  event  is  contained  in  an  llxuxll 
lattice.  To  contain  a 100-eV  copper  event,  which  had  the 
largest  yield,  a 13x4x13  lattice  was  needed.  Sputtering 
mechanism  traces  revealed  the  importance  of  boundary  atoms 
in  deflecting  and  channeling  atoms  in  the  sputtering 
process.  3y  choosing  the  PKA  to  be  an  atom  other  than  the 
one  at  the  surface  layer  center,  the  effect  of  an  increased 
lattice  size  could  be  duplicated. 

Investigation  with  a 17x4x17  microcrystallite  revealed 
that  Gibson-II  potential  function  simulations  are  more 
difficult  to  contain  than  when  using  the  KSE-3  potential 
function.  Only  the  100-eV,  copper  event  was  contained  using 
Gibson-II,  while  using  K3E-3  the  100-eV  and  200-aV,  copper 
events  were  contained.  It  is  quite  obvious  that  to  contain 
any  higher-energy  events  would  require  a microcrystallite  so 
large  that  the  computer  run  time  would  be  prohibitive. 
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ANGULAR  DISPERSION 


Because  the  low-energy  events  resulted  in  so  few 
sputtered  atoms,  the  angular  dispersion  results  for  any 
particular  ion  are  inconclusive.  A trend  of  light  isotope 
enrichment  is  observed  when  the  results  from  numerous  runs 
are  combined.  For  example,  the  histograms  of  the  four  sets 
of  36  randomly  selected  impact  points  show  a definite  shift 
of  the  lighter  isotopes  to  smaller  ejection  angles.  In  his 
report  on  the  enrichment  of  lighter  isotopes  in  the  normal 
direction,  Wehner  [24]  noted  only  a 0.6%  enrichment  in  Cu 
(63).  The  number  of  sputtered  isotopes  was  not  sufficient 
to  detect  this  small  difference. 

The  "reflective  collision"  theory,  which  Wehner  suggests 
to  be  the  cause,  was  observed  to  have  very  little  effect  on 
the  sputtered  atoms.  Investigation  revealed  that  in  most  of 
the  simulation  runs,  the  same  atoms  sputtered  independent  of 
its  isotope  type  or  the  various  positions  of  the  lattice 
isotopes.  When  the  PKA  was  a light  isotope,  the  atoms 
that  did  sputter  were  ejected  more  nearly  normal  to  the 
surface.  Varying  the  isotope  type  of  the  PKA,  atoms  that 
were  known  to  sputter,  and  some  of  the  target  atoms  that 
they  collided  with  resulted  in  differences  in  the  ejection 
angle  normal  of  less  than  0.4°. 

D.  RING  FORMATION 

Results  definitely  show  that  the  sputtering  mechanisms 
do  produce  a formation  of  rings  of  sputtered  atoms  around 
the  impact  area.  The  results  indicate  that  the  transfer  of 
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momentum  along  the  <100>  and  <101>  directions  in  the  surface 
layer  is  relatively  easy.  The  predominance  of  sputtering  of 
next-nearest  neighbors,  in  the  < 1 0 1 > direction  is  evident. 
For  example,  at  100-eV,  all  simulations  using  an  ion  of  the 
mass  of  argon  or  greater  sputtered  the  same  atoms,  but  with 
different  freguencies.  The  particular  atoms  that  sputter  do 
not  appear  to  be  dependent  on  the  mass  of  the  incident  ion. 
However,  the  number  of  times  the  atom  sputters  does  appear 
to  be  dependent  on  the  ion  mass.  Ring  and  atom  sputtering 
probabilities  clearly  snow  the  sputtered  atoms  from  certain 
site  directions  and  the  distances  of  the  sites  from  the 
center  of  the  PKA. 

A trend  was  noted:  as  the  energy  and  mass  of  the  ion  is 
increased,  the  sputtered  atoms  come  from  sites  which  appear 
to  spread  outward  from  the  impact  area.  These  higher-energy 
events  give  the  appearance  of  a "liquid  splash",  with  the 
size  of  the  crater  dependant  on  the  ion  energy  and  mass. 

S.  ION/ATOM  MASS  RATIO 

The  results  of  this  investigation  do  not  support  a 
monatonic  increase  in  the  sputtering  ratio  with  increased 
ion  mass  as  implied  in  figures  of  Sigmund[28]  and  Ser  Misse 
[27].  All  the  plots  showing  the  sputtering  ratio  as  a 
function  of  the  ion  mass,  clearly  indicate  a maximum  near  an 
ion/atom  mass  ratio  of  one.  At  ion/atom  mass  ratios  greater 
than  one,  the  sputtering  ratio  appears  to  become  asymtoptic 
to  a lower  value.  This  may  be  an  artifact  of  the  computer 
model.  These  plots  also  highlight  the  containment  problem. 
The  behavior  of  the  sputtering  ratio,  for  larger  ion/atom 
mass  ratios,  is  erratic  because  the  microcrystallite  does 
not  contain  the  event. 
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The  general  shape  of  the  sputtering  ratio  curve,  when 
plotted  as  a function  of  the  ion/atom  mass  ratio,  is  in 
agreement  with  the  experimental  data  of  Wehner  and  his 
co-workers.  Results  conclusively  show  that  the  ion/atom 
mass  ratio  is  an  important  factor  in  the  determination  of 
the  sputtering  ratio. 

F.  RANDOMLY  SELECTED  IMPACT  POINTS 

The  sputtering  ratio  obtained  using  the  grid  of  impact 
points  is  extremely  close  to  that  obtained  when  randomly 
selected  impact  points  were  used. The  similarity  between  the 
sputtering  summaries  and  the  atom  sputtering  probabilities 
are  quite  evident.  The  quantity  of  data  obtained  with  the 
randomly  selected  impact,  points  was  required  to  confirm  that 
light  isotope  enrichment  was  present. 

Based  on  the  comparison  of  results  obtained  from  the 
randomly  selected  impact  points  and  from  the  grid  of  impact 
points,  the  systematic  method  of  "walking”  through  the  grid 
seems  completely  justified. 


G.  FUTURE  RESEARCH 
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The  (110)  and  (111)  orientations  need  to  be  studied  in 
depth.  The  few  sputtering  mechanisms  observed  were  the  same 
as  those  for  the  (100)  orientation.  However,  the  sputtering 
ratio  as  a function  of  the  ion  mass,  containment 
requirements,  ring  formation,  and  angular  dispersion  of 
isotopes  need  further  investigation.  It  is  highly 
recommended  that  the  effect  of  momentum  transfer  along 
certain  directions  be  investigated. 


/ 
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APPENDIX  A 


ISOTOPE  GENERATION 

I 

^ All  atoms  in  the  lattice  were  one  of  two  types;  type  1 - 

the  heavier,  least  abundant  isotope,  or  type  2 - the 
lighter,  most  abundant  isotope.  The  average  atomic  weight 
of  copper  is  63.54  amu.  The  isotope  weights  are  62.93  amu 
and  64.93  amu  with  natural  abundances  of  69.1%  and  30.9% 
respect ively . 

The  lattice  was  generated  with  either  the  isotope  masses 
and  their  corresponding  abundance  or  the  weighted  average  of 
the  isotope  masses  with  an  abundance  of  100%.  Initially, 
all  atoms  were  forced  to  the  type  2 atom  in  the  lattice 
generator.  If  the  isotope  lattice  was  desired,  the 
International  Mathematical  and  Statistical  Library  (I3ML) 
routine  GGUB  was  called.  This  routine  generates  a uniformly 
distributed  set  of  psuedo  - random  numbers  between  zero  and 
one.  Each  of  these  random  numbers  was  assigned  to  a 
particular  atom  and  was  tested  to  see  if  it  was  less  than  or 
equal  to  0.309.  If  the  number  was  less  than  this  test  value 
(the  abundance  of  the  least  abundant  isotope) , the  atom  was 
changed  to  a type  1 atom.  A seed  for  the  random  number 
generator,  for  the  particular  size  lattice  to  be  generated, 
was  chosen  by  trial  and  error  which  would  produce  the 
correct  number  of  atoms  above  and  below  the  test  value  to 
assure  the  correct  isotope  ratio. 


i 
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A feature  was  designed  to  allow  any  number  of  the  atom 
masses  (types)  to  be  selectively  changed  after  the  lattice 
had  been  generated.  This  ' feature  was  incorporated  to 
investigate  the  changes  and  effects  of  picking  heavy/light, 
heavy/heavy,  or  light/light  collision  atoms.  This  would 
facilitate  the  dynamic  investigation  of  the  "reflective 
collision"  theory  proposed  by  Winters  and  Sigmund  [101. 

The  dynamics  of  the  simulation  remained  unchanged  except 
that  all  calculations  that  were  dependent  on  target  atom 
mass  were  modified  to  use  the  correct  mass. 


The  simulation 
with  the  correct 
selected  sites. 


was  now  capable 
isotope  ratio 


of  generating  a 
and  masses  with 


1 attics 
randomly 


APPENDIX  B 
LATTICE  FOLDING 


As  mentioned  in  the  simulation  model  development,  a 
lattice  that  was  symmetrical  about  the  folding  axes  was 
required.  Once  this  was  achieved  by  the  lattice  generators, 
it  was  necessary  to  fold  the  lattice  in  a manner  such  that, 
after  the  folds,  the  reduced  representative  impact  area  was 
identical  to  the  complete  representative  impact  area.  After 
these  folds'  were  performed  correctly  and  in  the  right 
sequence,  the  results  were  identical  to  that  obtained  if  the 
complete  representative  impact  area  had  been  bombarded 
systematically. 

After  each  fold,  the  sputtering  count  for  all  atoms  that 
were  symmetric  about  the  folding  axis  were  summed.  The 
summing  process  started  with  the  first  atom  in  a layer  to 
the  highest  numbered  atom  in  that  layer  After  this  process, 
the  sputtering  count  for  each  atom  was  twice  what  the  sum 
was  since  the  sputs  were  summed  each  time  a symmetric  atom 
was  encountered.  This  doubling  effect  of  each  fold  was 
corrected  by  scaling  the  sputtering  count  of  each  atom  after 
all  the  folds  in  a layer  were  completed.  All  necessary 
folds  were  then  performed  on  each  successive  layer. 
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(100)  ORIENTATION  FOLDING 
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Figure  12  shows  the  (100)  orientation  reduced 
representative  impact  area  folding  geometry.  The  reduced  36 
impact  point  area,  when  folded,  represents  the  complete  238 
impact  point  area. 

Since  the  (100)  orientation  impact  area  is  a 45°-45° 
triangle,  three  folds  were  required  to  complete  the  folding 
process.  The  first  fold  was  along  the  diangon  The  first 
fold  was  along  the  diagonal  axis,  the  second  along  the 
horizontal  axis,  and  the  third  along  the  vertical  axis.  The 
diagonal  fold  process  for  all  layers  was  the  same.  However, 
the  horizontal  and  vertical  folds  were  modified  depending  on 
which  layer  was  being  folded. 

After  all  the  folds  in  a layer  for  the  (100)  orientation 
were  complete,  the  sputtering  count  for  each  atom  was 
divided  by  eight  to  correctly  scale  the  results. 


B.  (110)  ORIENTATION  FOLDING 

Figure  12  also  shows  the  (110)  orientation  reduced 
representative  impact  area  folding  geometry.  The  reduced  36 
impact  point  area,  when  folded,  represents  the  complete  144 
impact  point  area. 

Since  the  (110)  orientation  reduced  impact  area  is  a 
rectangle,  two  folds  were  required  to  complete  the  folding 
process.  The  first  fold  was  along  the  horizontal  axis  and 
the  second  fold  was  along  the  vertical  axis.  The  diagonal 
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fold  was  not  necessary.  The  horizontal  and  vertical  folds 
(exactly  the  same  as  for  the  (100)  orientation)  were 
modified  depending  on  which  layer  was  being  folded. 

After  all  the  folds  in  a layer  for  the  (110)  orientation 
were  complete,  the  sputtering  count  for  each  atom  was 
divided  by  four  to  correctly  scale  the  results. 


C.  (HI)  ORIENTATION  FOLDING 

The  (111)  orientation  was  by  far  the  most  difficult  to 
determine.  Figure  13  shows  the  (111)  orientation  reduced 
representative  impact  area  folding  geometry.  The  reduced  21 
impact  point  area,  when  folded,  represents  the  complete  126 
impact  point  area. 

Since  the  (111)  orientation  impact  area  is  a 30°-60° 
+ riangle  three  folds  along  the  diagonals  resulted  in  the 
overlap  shown  in  Figure  13.  However,  if  the  sat  of  three 
folding  sequences  were  summed  the  results  could  be  scaled  to 
represent  the  complete  impact  area.  The  number  one  fold  was 
along  the  positive  slope  diagonal  axis,  the  number  two  fold 
was  along  the  negative  slope  diagonal  axis,  and  the  number 
tnree  fold  was  along  the  vertical  axis. 

Before  each  sequence  of  folds,  the  sputtering  count  for 
each  atom  was  reinitialized  to  the  original  value.  The 
first  sequence  (1,2,3)  was  performed  and  upon  completion  the 
sputtering  count  for  each  atom  was  divided  by  eight  to 
correctly  scale  the  results.  The  second  sequence  (1,2,1) 
was  performed  and  again  the  sputtering  count  for  each  atom 
was  divided  by  eight  for  scaling.  The  third  sequence 
(1,3,2)  was  performed  and  again  the  sputtering  count  for 
each  atom  was  divided  by  eight  for  scaling. 
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Finally,  the  results  obtained  from  each  sequence  was 
summed  and  the  final  sputtering  count  for  each  atom  was 
divided  by  four  to  correctly  scale  the  results  due  to  the 
summing  process. 


Since  the  geometry  of  the  (111)  orientation  was 
complicated  and  not  fully  developed,  the  folding  process  was 
performed  only  for  the  first  layer. 


APPENDIX  C 


GLOSSARY 


ABUN1,  ABDN2 

The  abundance  of  the  target's  isotopes. 

AC 

A distance  measurement  used  in  impact  point 
generati on . 

ALFA,  ALFA2 

The  cosine  of  the  angle  between  vectors  31, R2, 
sq  uared. 

ALPHA 

Input  Morse  potential  parameter. 


ANGLE 

Th<-  angle  between  the  direction  of  motion  of  an 
ejected  atom  and  the  surface  normal  in  degrees. 

AVE 

The  sputtering  ratio. 

AXT,  AYT,  AZT 

Components  of  the  random  thermal  displacement. 

BENRGY 

The  lattice  binding  energy  in  ev. 
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BIND 


Negative  of  the  total  potential  energy  (TPOT)  at 
time  zero.  Retained  for  energy  conservation 
checks . 

BK  AS 

The  mass  of  the  ballet  in  amu. 

BQZ 

A fee  crystal  constant  -(1/ROOT (3)). 

BOL 

The  bullet  heading  Hollarith. 

BX,  BZ 

Unsealed  x,z  coordinates  of  the  impact  point. 

BXNAX , BXHIN 

X direction  boundaries  of  the  target  impact  area. 

BZHAX,  BZHIN 

Z direction  boundaries  of  the  target  impact  area. 

CELS 

A frictional  force  multiplier  when  frictional 
forces  are  included  in  the  program. 

CGB1 , CGB2 

Morse  potential  parameters. 

CGD1 , CGD2 

Morse  potential  parameters. 

CGP1 , CGF2 

Morse  force  function  parameters. 
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CFO,  CF 1 , CF2 

Force  parameters  of  the  cubic  fit  between  Morse 
and  Born-Mayer  functions. 


CO,  COHF 

"Phonon"  displacement  correlation  parameter  used 
as  an  option  in  HARM. 


CORN 


1 - CORF  . 


COX,  COY,  COZ 

Directional  cosines  of  the  incoming  ion's  velocity 
vector,  relative  to  the  surface  normal. 

COXI,  COYI , COZI 

Negative  values  of  COX,  COY,  COZ. 

CPO,  CPI,  CP2,  CP 3 

Potential  parameters  of  cubic  fit  between  Morse 
and  Born-Mayer  functions. 

CVBD 

A ratio  used  to  avoid  repeated  division. 

CVD 

Distance  conversion  factor  - angstroms  to  meters. 

CVE 

Energy  conversion  factor  - ev  to  joules. 


A ratio  used  to  avoid  repeated  division,  CVF/CVD, 
used  to  convert  exponential  potential  functions  to 
forces . 
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J 


CVED 


cvn 

f 

Mass  conversion  factor  - amu  to  kg. 


CVR 

Conversion  factor  - lattice  units  to  angstroms. 

\ 

DBXHAX,  DBZMAX 

Upper  boundary  limits  for  the  target  impact  area. 

DCON 

Input  Morse  potential  parameter. 

■>  DFP 

: m . 

The  distance  closer  than  ROE  that  an  atom  is  to 
the  bullet.  Used  when  forces  are  turned  on  or  off. 

DIST 

The  distance  between  any  two  atoms. 

DNN  (I) 

Distance  to  the  Ith  ring  in  lattice  units. 

**** 

DRX,  DRT,  DRZ 

Components  of  DIST. 

DT 

Length  of  each  timestep  in  seconds. 

DTI 

Number  of  lattice  units  the  most  energetic  atom 
can  move  in  one  timestep. 

DTOD 

DT  over  distance  - a ratio  used  to  avoid  repeated 
division . 

Ii 
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DT0M1,  DT0B2 

DT  over  target  isotope  mass  - a ratio  used  to 
avoid  repeated  division. 


DTOMB 

DT  over  bullet  mass  - a ratio  used  to  avoid 
repeated  division. 


DX(I)  , DT  (I)  , DZ  ( I) 

Change  in  the  position  of  the  Ith  atom. 

EMAX 

The  velocity  squared  of  the  fastest  atom.  Used  to 
make  ETEST. 


ETEST 


The  kinetic  energy  of  the  most  energetic  atom. 


ETHR 


Energy  threshold  to  escape  the  surface. 


EV 

Primary  energy  of  the  bullet  in  ev. 

EVR 

Primary  energy  of  the  bullet  in  kev.  R indicates 
a read-in  and  saved  value. 


EXA,  EX3 

Input  Born-dayer  potential  function  for  the  target 
atom  - target  atom  interaction. 


FA 

The  temporary  component  force  increment  on  an 
atom . 
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PIC 

The  decrement  size  used  in  SETUP  in  locating  the 
ion  impact  point. 

FINISH 

The  elapsed  program  running  time  in  seconds. 

FLAT 

A lattice  unit  in  angstroms. 

FM 

A small  number  used  in  checking  potential  energy 
zero  point. 

FM  AX 

The  maximum  frequency  of  occurrence  in  the 
histogram  intervals. 

FORCE 

A numerical  value  of  the  force  function  for  a 
specific  pair  of  atoms. 

FPTC 

The  corrective  force  at  ROE. 

FBC 

Numerical  value  of  the  target  force  function  at 
FOE. 

FRAN 

A real  variable  for  number  of  passes  completed  at 
termination . 

FRANR 

The  reciprocal  of  FRAN. 


75 


FHEQ1,  FREQ2 

Humber  of  events  in  each  interval  of  the 
histogram. 

FSCAL 

A histogram  scale  factor. 

FTDI 

A parameter  used  to  determine  DT  by  maximum  energy 
method. 

FXA 

4 

A Born-Mayar  force  function  parameter. 

?X(I)  , FY  (I)  , FZ  ( I) 

Components  of  the  total  force  on  the  Ith  atom. 

HBHAS 

Half  Bullet  Mass  a ratio  used  to  avoid  repeated 
division . 

HDTOD 

Half  DT  over  distance  - a ratio  used  to  avoid 
repeated  division. 


HDTOM 

Half  DT  over  most  abundant  isotope  mass  - a ratio 
used  to  avoid  repeated  division. 

HDTOMB 

Half  DT  over  bullet  mass  - a ratio  used  to  avoid 
repeated  division. 

HTHA§ 1 , HTHAS2 

Half  target  isotope  mass  - a ratio  used  to  avoid 
repeated  division. 
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IBOX 1 , IBOX2 

Ejected  atom  angle  tally  keeper. 

ICOUNT 

Flag  used  in  sequencing. 

IDEEP 

The  first  frozen  layer.  in  the 
microcrystal . 

IET 

Elapsed  time  in  2.6  microsecond  blocks, 
determine  actual  elapsed  time. 

IFFH 

The  atom  number  of  any  target  atom  that 
mass  constrained  due  its  isotope  ratio. 

IFLAG 

Flag  used  in  sequencing. 

IFM 

The  isotope  type  (heavy  or  light)  of  any 
atom  that  has  its  mass  constrained  due 
isotope  ratio. 

IFT 

A number  used  as  a constant  when  gen 
pseudo-random  numbers. 

IH1 

Alphanumeric  array  of  the  program  heading. 

IH2 

Alphanumeric  array  of  the  Morse  parameters. 


target 


Used  to 


has  its 


target 
to  its 


erating 
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IHB(1)  , 


IBB  (4)  , 


IHB  (2)  , IBB  (3)  , 


IHB  (5)  , IHB  (6) 


Alphanumeric  array  describing  the  bullet. 


IHBPOT 

The  type  of  ion-atom  potential. 

IHBOLL 

The  bullet  elemental  symbol  or  name. 

IHS 

Alphanumeric  array  for  the  type  and  orientation  of 
the  crystal. 

IHT 

Alphanumeric  array  describing  the  target. 

IHT  ( 1)  , IHT  (2)  , IHT  (3),  IHT  (4)  , IHT  (5)  , IHT  (6) 

Alpha  numeric  array  describing  the  target. 


IHTARG 

The  target  elemental  symbol  or  name. 


IHTPOT 

The  type  of  target  atom  - target  atom  potential. 


IL  ATT 

A flag  used  in  the  lattice  generators. 


HATER 

A variable  used  to  identify  wnich  layer  is  being 
folded . 


ILL 


The  number  of  the  last  atom  in  the  first  layer. 
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IMP 

The  mass  of  each  target  atom  with  the  correct 
isotope  ratio. 

IN 

The  number  of  intervals  in  the  histograms. 

INDEX 

A variable  used  in  determining  phase  in  a dynamics 
cycle  step. 

INT 

The  width  of  the  intervals  in  the  histograms. 

IOH 

The  target  plane  of  incidence. 

IOOT 

Flag  used  in  sequencing. 

IPASS 

The  number  of  passes  (shots)  completed. 

IPNCH 

The  impact  point  number. 

ipnons 

Scaled  impact  point  number. 

IPROG 

Input  variable  used  to  skip  attractive  potential. 

IRAN 

Random  number  seed  used  in  NORMAL  subroutine  which 
generates  pseudo-random  numbers. 
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IRANGE 

The  upper  limit  in  the  histograms. 

IRBAD 

The  number  of  atoms  which  have  the  isotope  type 
constrained  externally. 

ISEEDA,  ISEEDB,  ISEBDC 

Seeds  used  in  the  uniform  random  number  generator 
which  determines  the  isotope  ratio  and 
distribu  tion . 

I SET 

Flag  used  in  sequencing. 

ISHOT 

A parameter  used  to  shut  down  the  program. 

ISNUH 

The  atom  number  of  a sputtered  atom. 

ISOGEN 

A flag  used  as  an  option  to  generate  the  target 
atom  masses  with  isotope  or  single  element  values. 

ISPX,  ISPZ 

Initial  start  point  x,z  coordinates. 


IT 

Unsealed  x coordinate  used  in  lattice  generation. 

ITEHP 

The  lattice  temperature  in  degrees  Kelvin. 

ITIME 

Flag  used  in  sequencing. 
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ITT 

An  odd-even  integer  used  to  determine  atom  site 
location . 

IX,  II,  IZ 

The  number  of  x,y,z,  planes  in  the  crystal. 

IXP,  IYP,  IZP 

Crystal  dimensions  in  x,y,z. 

JN 

The  numerical  value  of  the  histogram  intervals. 

JOOT 

A histogram  sequencing  variable. 

JSCAL 

A histogram  scale  factor. 

JT 

Unsealed  y coordinate  used  in  lattice  generation. 

JTS 

Variable  used  to  establish  atom  site. 

JTT 

Variable  used  to  establish  atom  site. 

K2S0M 

The  square  of  the  sputtering  count  which  is  used 
in  calculating  the  variance. 

KATOH 

The  number  of  atoms  in  a ring. 
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A sequencing  variable. 


KFIRST 


KHING 

The  highest  numbered  ring  in  a layer. 

KSPOT 

The  number  of  sputs  from  all  atoms  in  a ring. 

KSUM 

The  total  sputtering  count  before  folding. 

KT 

Unsealed  z coordinate  used  in  lattice  generation. 

LCUT  (I) 

Used  to  identify  an  I th  atom  which  can  be  omitted 
from  the  force  and  potential  calculations  by 
changing  LCUT(I)  from  zero  to  one. 

LD 

The  highest  numbered  atom  in  a mobile  layer. 

LL 

The  highest  numbered  atom  in  tne  entire  crystal. 

LL2,  LL  3 

Twice  and  triple  the  value  of  LL . 


F 


k 

► 


LNOMI<I) 

The  number  of  times  the  Ith  atom  sputters  after 
all  shots  but  before  folding  the  summary. 

LS 

Sum  of  the  Miller  index  integers. 


LSS 


Variable  used  to  identify  the  type  of  surface. 


LSOH 

The  lattice  surface  type. 

LTYPE 

The  type  of  lattice  (i.e.  fee  ) . 

MAX 

A height  scale  factor  used  in  the  histograms. 

HCBO 

The  number  of  simultaneous  equations  solved  by 
CHOSYM. 


N 1 

NL AY  - 1,  a sorting  variable. 

NCT 

Increments  the  timestep  interval. 

s 

NCTT,  ND 

Data  output  interval  in  number  of  timesteps.  Also 
the  timestep  interval  between  potential  energy 
calculat ions . 
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NFIRST 


NLAY 


NL  AY  1 


NLINE 


NOTH 


NPAGE 


NR  AN 


A sequencing  variable. 


The  number  of  the  first  atom  in  a gi.ven  layer. 


The  highest  numbered  atom  in  a given  layer. 


The  highest  numbered  atom  in  the  first  layer, 


The  line  number  in  the  output  block. 


A Hollar  it h consisting  of  blanks. 


A page  numbering  variable. 


Variable  for  number  of  passes  (shots)  completed  at 
termination . 


A restart  flag. 


The  next  timestep  number  at  which  output  will  be 
printed. 


HSHOTS 


The  number  of  passes  (shots)  to  be  completed  in  a 
program.  Used  to  break  up  an  extremely  large, 
many-shot  program  into  a few  small  programs  with 
manageable  computer  run  time. 


NSPOT  (I) 

The  number  of  atoms  that  sputter  for  the  Ith 
impact  point. 


NSR 

The  first  timestep  where  output,  is  printed. 

NSS 

Flag  to  print  final  output  status  of  current  shot. 

NT 

The  timestep  number. 

NTT 

Timestep  number  limit  before  shut  down. 

PAC 

Parameter  for  bullet  force  function  correction. 

PBMAS 

Bullet  mass  in  kg. 

PEXA,  PEXB 

Input  Born-Mayer  potential  function  parameters  for 
the  bull et- target  interaction. 


PFIV 


Point  fi ve  (0.5). 
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PFPTC 

Ballet  force  function  evaluated  at  HOE. 


PFRC 

Bullet-atom  repulsive  force  interaction  evaluated 
at  ROE. 
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PFX& 

Bullet  force  function  parameter. 

PKE  (I) 

The  kinetic  energy  of  the  Ith  atom. 

PNEG 

Total  negative  potential  energy  in  the  lattice. 

PNOU 

Alphanumeric  variable  for  the  impact  point  number. 

POT 

The  potential  energy  between  two  atoms. 

PPE(I) 

The  potential  energy  of  the  Ith  atom. 

PPOS 

Total  positive  potential  energy  in  the  lattice. 

PPTC 

Bullet  potential  function  evaluated  at  ROE. 

PBOB 

The  probability  that  an  atom  will  sputter  from  a 
given  ring. 


La 
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PTC 

Target  atom  potential  function  evaluated  at  ROE. 

PTE  (I) 

The  total  energy  of  the  Ith  atom. 

QM 

A small  number  used  in  checking  the  kinetic  energy 
zero  point. 

QOIT 

Coutoff  variable  checked  against  the  total 
potential  energy  which  can  be  used  as  a shut  down 
criterion. 

HI,  HIS 

Distance  from  (0,0)  impact  point  to  the  initial 
bullet  position. 

HI,  R IT  , R 1 Z 

The  x,y,z  components  of  R1S. 

H2 

Magnitude  of  vector  from  impact  point  to  the  first 
atom  hit  by  the  bullet. 

B3 

Magnitude  of  vector  from  impact  point  to  bullet 
start  position. 

HI,  R2 

Scalar  product  of  vectors  HI,  R2. 

R2SQ 

R2  squared. 
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RAD 

The  angle  between  the  direction  of  an  ejected  atom 
and  the  surface  normal  in  radians. 


RADDE6 


Conversion  factor  - radians  to  degrees. 


RAN  (I) 


\ list  of  normally  distributed  random  numbers. 


RANGE 

A scaling  variable  used  in  the  histogram. 

RBX,  RBZ 

Unsealed  x,z  coordinates  of  the  impact  area 
reference  point. 

RDIST 

The  distance  between  an  atom  and  the  center  point 
of  its  layer. 

RE 

Input  Norse  potential  parameter. 

RINCR 

The  RANGE  increments. 


RHIDPT 

The  midpoint  of  each  interval  in  the  histogram. 

ROE,  R0E2 

Nearest  neighbor  distance  - squared. 

ROEA 

The  maximum  cutoff  for  the  Born-Mayer  potential. 
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BOEB 

The  minimum  cutoff  for  the  Morse  potential. 

BOEC,  BOEC2 

The  maximum  cutoff  for  the  Morse  potential 
squared. 

BOEM 

Defines  the  region  where  the  repulsive  force 
function  must  be  modified. 

BX  (I)  , BT  (I)  , RZ  (I) 

X,y,z  coordinates  of  the  I*h  atom  at  any  time. 

HXBL,  RYBL,  BZBL 

Maximum  negative  x,t,z  coorlinates  an  atom  aav 
move  before  it  i3  LCUT. 

R X BH  , RYBH,  R ZB  M 

Maximum  positive  t,/,z  cooriina**s  an  a*’?i  aay 
move  before  it  is  LCUT. 

SIP,  RYP,  RZF 

5inal  x,y,z  coordinates  of  the  Ith  a^om  w sen  i * s 
LCUT. 

RXI  (I)  , RYI  (I)  , RZI  (I) 

The  Ith  atom's  initial  x,y,z  coordinates. 

RXK  (I)  , BYK  (I)  , RZK  (I) 

Temporary  x,y,z  coordinates  of  the  Ith  atom  luring 
the  force  cycle. 

BIS,  BYS,  BZS 

Scaled  x,y,z  coordinates  of  the  impact  point. 
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SAVE 

One-half  POT. 

SBX,  SBZ 

Scaled  x,z  coordinates  of  the  impact  point. 

SCALE 

A scaling  variable  for  the  frequency  of  events  in 
the  histogram. 

SCI,  SCT , SCZ 

The  x,y,z  lattice  scale  factors  used  to  convert 
tfe  gent-rated  lactic*  to  ma*-ch  a real  crystal. 

SCXP,  SCZR 

The  inverse  of  SCX  and  SCZ  - a ratio  used  to  avoid 
repeated  livision. 

SD  1ST 

The  square  of  t he  distance  from  the  center  of  a 
layer  to  an  a*om. 

SETIHE 

A library  subroutine  used  to  keeD  track  of  run 
time  . 

SLOW 

A cutoff  velocity  variable  checked  against  a long 
DT. 

SPDX , SPDZ 

The  fraction  of  the  target  impact  area  that  the 
impact  points  are  incremented. 
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SPZ 


SPX, 

The  x,z  components  of  the  distance  from  the  impact 
area  reference  point  to  the  impact  point. 

sscz 

A z scale  factor  used  in  the  (111)  orientation 
lattice  generation. 

STD 

The  standard  deviation  of  the  sputtering  ratio. 

STOPER 

This  variable  stops  calculating  and  saves  60 
seconds  for  compilation  and  output. 

SOB 

Alphanumeric  variable  describing  the  surface. 


TA,  TABP 

The  thermal  amplitude  of  the  lattice. 

TIB 

Alphanumeric  variable  describing  the  target 

element . 

TARGET 

Alphanumeric  variable  describing  the  target 

element . 

TAS 

A scale  factor  - TAC1P/R00T  (3  .0)  . 

TASF , TASN 

Random  thermal  amplitude  scale  constants. 
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TE 

The  total  kinetic  and  potential  energy  of  all  the 
atoms. 

TEMP 

The  temperature  of  the  lattice  in  degrees  Kelvin. 

TESTE 

A termination  test  of  energy. 

TFAC 

A time  factor  ratio  used  to  determine  DT  by 
maximum  force  methods. 

THERM 

Mean  thermal  energy  of  a lattice  atom.  Dsed  in 
shut-down  tests. 

TIME 

The  elapsed  calculating  time  in  seconds. 

TITLE1 , TITLE2 

The  titles  for  the  histogram  outputs. 

TH AS  1 , TMAS2 

The  target  isotope  masses  in  amu. 

TPKE 

The  total  kinetic  energy  of  all  the  atoms. 

TPOT 

The  total  potential  energy  of  all  the  atoms. 

THUN 

The  time  the  program  is  to  run  in  minutes. 


I 
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f 


f. 

I 


TEL 

The  velocity  calculation  - SQET (ev/HBMAS) . 

vss 

A temporary  storage  variable  for  the  velocity 
componen  ts. 

VTEST 

A test  velocity  used  to  determine  if  the  atom  has 
sufficient  velocity  (energy)  to  escape  the 
surface. 

7TS 

A scale  factor  - 0.8*VTEST. 

VX(I)  f VY(I)  , VZ(I) 

The  x,y,z  components  of  the  Ith  atom's  velocity. 

VXS,  VYS , VZS 

The  x,y,z  components  of  the  bullet's  initial 
velocity . 

XDIST 

The  x component  of  the  distance  from  the  center  of 
a layer  to  the  atom. 

XH AX , X MIN 

The  x planes  to  determine  whether  two  atoms  are 
close  enough  to  interact. 

YLAX  (I) 

The  crystal  relaxation  of  the  Ith  layer  in  the  -v 
direction. 
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THAI , 


ZDIST 


ZE 


ZM  AX  , 


ZSTOP 


THIH 

The  y planes  to  determine  whether  two  atoms  are 
close  enough  to  interact. 

The  z component  of  the  distance  from  the  center  of 
a layer  to  the  atom. 


Absolute  zero  (0.0). 

ZHIN 

The  z planes  to  determine  whether  two  atoms  are 
close  enough  to  interact. 


Used  to  print  out  a warning  header  that  a pass  was 
terminated  before  completion. 
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APPENCIX  D 


COMPUTER  PROGRAM 


GLD3  IS  OLD!  MCCIFIEC  TC  INCREASE  SPEED 
ROLD  IS  GLD3  WITH  'WARMEC'  LATTICE 
RING  IS  RPF  MODIFIED  TC  FCLD  LATTICE  ANC 
STUDY  RING  FORM ATI  ON /SPUTTERING 
RMAS  IS  RING  MOCIFIEC  TC  STUDY  TARGET  ISOTCPES 
RANG  IS  RMAS  MOD.  TO  STUCY  ANGULAR  DISPERSION 


DIMENSION  DX(600),CY(600),DZ(60C) 

DIMENSION  RXK16 JO)  ,RYK(6JJ)  ,RZK(600) 

DIMENSION  PKE< 600 ) ,PTE (600) 

DIMENSION  NSPUT( 50)  , I PNUMSt 50  ) 

DIMENSION  I $NUM(2C)  , ISEG1600  ) 

DIMENSION  RXF(6GC),RYF(6G0),RZF(6CC) 

DIMENSION  180X1(20),  IBGX2(20),  RANNUM(250) 

DIMENSION  FPEQ1 180) ,FREG2(30) ,RANGE(20)  * T I T LE  1 ( 24  ) , 

X T ITLE2 (24  ) 

DIMENSION  P CTOM ( 2 ) 

DIMENSION  IFFM(60C)  ,IFM(600  ) 

CCMM0N/CCM1/RX(60C),RY(600) ,RZ(60C),LCUT(60C),LL,LD 
CCMMCN/CCM2/R0E, RCE2, RGEM,AC,PAC  , PPTC , PT C , P FPTC , F PTC , 

X FM,PFIV,TPCT 

COMMON/C GM 3/E XA,EXB,FXA,PEXA,PEXS,PFXA, THERM 
CCMMCN/C0M4/IX,IY , I Z , IXP,IYP,SCX,SCY,SCZ,IDEEP, 

X R E X , RE  Z 

COMMON/ COM4 A/  LNUM(6J0),LNUMI(600) ,LNUMF (600) 
CCMMCN/C0M5/FX( 60  0 ) ,FY(60Q)  ,FZ ( 6CC ) , PP5 ( 600  ) 
CCMMCN/C0M6/C0XI , COY  I ,CCZI ,RXS,RYS,RZS,FAC 
CCMMCN/C CM 7 /R 1, LSS ,SPX,SPZ, COX, C CY, COZ, YLAX  (20), ILL 
CCMMGN/CQM8/VX ( 60C i ,VY(600),VZ(6GC) 
CCMMCN/CGM8A/VXS(60G) ,VYS(600) , V Z S ( 600 ) 

COMMON/ COM!  I/RXI (6GC),RYI(600)  ,RZ  I (600) 

COMMON/ COM  13/  RAN(18C0),IRAN,TASN,TASF,LL2,LL3 
COM MON/ COM  14/  AX T (6 J 0 ) , A YT ( 60 J ) , A ZT ( 6 0 0 ) 

COMMON/COM  16  / EM  AS  , TMAS  1 ,T  MA  S2  , A E UM  , ABUN  2 , I PNOM 
CCMMCN/CQM16A/3X, BY ,3Z,EDI S P * C VR 

CCMMCN/C0M17/TIME,NT,DT,DTI , TEMP, CELS, EV,N PAGE, NL I NE 
COM  MON/COM  1 8/  I TEMP  , TAMP  .CORF 

CCMM.CN/COM2  0/IH1  (20),IHS(10),IHT(6),IH3(6),IHTARG(21, 

X IhBUL  L ( 2 ) 

COMMON /COM21/IHBPOT (2 ) , IHTPCT ( 2 ) , SUR ( 2 ) , IH2 ( 10  ) , I OR  , 

X LSLR 

COMMON/ COM2  2/ ROE A ,RCE8,RCEC,R0EC2,CP0,CP1,CF2,CP3,CF0, 
X CF1,CF2,CGC1,CGD2,CGB1,CG3  2,CGF1,CSF2, I PROG 
CCMMCN/C0M23/  SPDX,SPCZ,DBXMAX,CBZMAX,SSCZ 
COMMON/ COM2  4/  DNN (100),KATOM(LOO),KSPuT(100> , FROB ( 100  ) 
CCMMCN/ C0M24A/  KSLM,KRING,RGIST(6CC) 

COMMON/ COM 2 5/ IMF ( 60  0) , ISOGEN 
C CM  MON/ COM  3 0/  ISEEDA,  I SEEDE , ISEEDC 
COMMCN/CCMA/  A(4,5),MCRC 
90  J 5 FORMAT  (20A4) 

90  1C  FORMAT ( 20A4  ) 

902C  FCRMAT(10A4,I3,2A4,I2,3I3,2F3.0,2F5.2,I2) 

9030  FORMAT  ( 2 A 4 , A4 , F 8 . 3 , 2 A 4 , 2 F 1 0 . 5 , F 8 . 5 , 3F 7 . 4 , I 3 ) 

9035  FORMAT  (F5.2,I5) 
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904C  FORMAT! 10A4.3F5. 2 ,3F8. 5,11  ) 

905C  FCRMAT  ( 2X , 13 , 15 ( 1 3 , IX , 1 1 , 2X) ) 

90 8 C FCRMAT(5X,3F5.3,4I5,2F8.6,2X,A4,3F6.4) 

9 1 C 5 FORMAT!  11X,F5.0) 

96CC  FCFMAT ( //  ) 

96C1  FORMAT!  1H1) 

96 2C  FORMAT ( 1X,3(I5,4X, I1,4X,3F6.2, 8X)  ) 

96  3 0 FCPMAT(135X,4HPAGE,I3,/,1H1  ) 

964C  FORMAT ( ' BINOING  ENERGY  =',F6.2,',  TEST  SPEED  = ', 

X F 10 . 2 » ' M/SEC' .SAFETY  FACTOR  SPEED  =',F1C.2, 

X • M/SEC ' » / »1H1 ) 

96 8C  FORMAT ( 2X . ' ATOM* , 2X , 1 LOUT • , 9X , • VX • .9X , ' V Y ' , SX  , • VZ  1 , 

X 8X,'RXF',8X,'RYF',8X,'RZF',8X,'PPE',8X,'PKE',8X, 

X 'PTE'/) 

9681  FORMAT! IX, 'ATOM  DX • , 6X, • D Y ' , 6X , • DZ • , 7 X , • V X • , 7X , • VY ' , 

X 7X,'VZ',7X,'KE',8X,'PE',8X,'TE',7X,'RX',6X,'RY',6X, 

X 'RZ', AX, 'LOUT  I SEO  1 , / ) 

96  8 2 FCRMAT ( I 5 , 3 F8 . 3 , 3 F9 . 1,3F10.4, 3F8.3, 14, 3X  , 14  ) 

96  65  FCRMAT(2X,I3,3X,I3,5X,3(F8.1,3X),6(F8.3,3X)  ) 

9690  FCRMAT!/4X, F10.3  , 'EV,  TOTAL  KINETIC  ENERGY  F10 .3  , 

X 'EV,  TOTAL  POTENTIAL  ENERGY F 1 0 . 3 ,' EV , TOTAL', 

X ' ENERGY, ' ,/9X,  ' TOTAL  POSITIVE  POTENTIAL  ENERGY  IS', 
X F 10 .3 , ' TOTAL  NEGATIVE  POTENTIAL  ENERGY  IS', FLO. 3) 
9695  FCRMAT (3X, ' ATOM  8 POTENTIAL  SUR F • ,5X , ' V X ' , 7 X , ' VY ' , 7X  , 

X • VZ * , 10X , 'PKE* ,6X, 'PPE ' ,6X,  ' PTE  KEV  OCX  COY',//) 
9701  FCRMAT ( • IFT  =',110,',  IRAN  =',110,',  LL=',I5,//) 

98  0 0 FORMAT! IX , 2 1 5 ,2 A4 , 14 , A4 , 1 X , 3F9 . 1 , 2 F9 .4 , F 5 . 1 ,1X,2F6.3) 
98C5  FORMAT!  5X,  ' SPUTTERED  ATOM  / ANGLE  WITH  NORMAL'/) 

9810  FORMAT!//, • AT',F7.1,'  EV,  THE  FOLLOWING  ATOMS  wERE', 

X • SPUTTEREC:'/ ,2014,//) 

9315  FORMAT! 1X,5!5X, I3,3X,F7.4) ) 

982C  FORMAT!  • THERE  ARE  NO  SPUTS  AT  ',F8.1,'EV') 

98 AC  F0RMAT(25X,  'SUMMARY  OF  TARGET  IMPACT  AREA'/,25X, 

X 'FOLDED  TO  REPRESENT  TOTAL  COVERAGE'/) 

9845  FCRMAT  (IX, 'ATOM  NUMBER  / TIMES  SPUTTERED  ') 

985C  FORMAT!///,  • SPUTTERING  RATIO  WAS',F6.3, 

X ',  THE  VARIANCE  WAS',F8.3,',  FOR  ',14, 

X ' IMPACT  POINTS,  THE  S^D  DEV  OF  THE  MEAN  WAS', 

X F6.3  ,/  / ) 

986C  FORMAT! 2X,  • STATISTICAL  SUMMARY  FOR • , 14 , ' I M P AC T POINTS' 
X ,'  AT  ' , F 5 • 1 , 'KEV'  ,///) 

9865  FORMAT!'  NUMBER  OF  TIMES  EACH  ATOM  SPUTTERED'/) 

98  7 G FORMAT  (10!  16, 14)  ) 

9875  FORMAT!//,'  BINDING  ENERGY  WAS',F6.2,'  EV,  TEST  SPEED' 
X ,'  WAS'  , F 1 0 . 2 , ' M/SEC,  SPUTTEREC  ATOM  COUNT  IS', 15) 
9880  FORMAT! IX  ,3  ( • RING  DISTANCE  ATOMS  SPUTS  PROB  ')) 
9885  FORMAT!  IX, 3!  ' NN ( • , I 2 , ' ) ' , 2X , F6 . 3 ,3X , 12 , 5X , 1 3 , 3X  , 

X F6.3.1X)  ) 

98SC  FORMAT (25X, 'PROBABILITY  THAT  A RING  OF  ATOMS  SPUTTERS' 
X //  ) 

9895  FCRMAT(25X, 'SPUTTERED  ATOMS  ANGULAR  DENSITY'/) 

99  C C FORMAT! 2X,  • INCREMENTS  ARE  IN',F4.2,'  DEGREE  INTERVALS' 

X /) 

9905  FCRMAT(2X,  'ANGLE  (DEGREES)  / NUMeER  OF  ATOMS'/) 

99  1C  FCRMAT  ( 2 X , 19 ( I 3 , 3X ) , 1 3 ) 

9915  FCRMAT  (3X,I2,2(  1X.F4.2)) 

9920  FORMAT!  • 00-05  05-10  10-15  15-20  20-25  25-3C  30-35', 

X ' 35-4C  40-45  45-50  50-55  55-6C  60-&5  65-7C  70-75', 

X • 75-80  30-85  85-90' /) 

9925  FORMAT!  • CO-LO  10-20  20-30  30-40  40-50  5C-6C  60-70', 

X • 70-8  C 80-90'/) 

99 3 C FORMAT!'  NUMBER  OF  ATOMS  CHANGED  WAS', 14/,'  ATOMS', 

X • CHANGEC  AND  ISOTOPE  TYPE  ARE') 

9935  FORMAT  ( IX , 13 ( I 3 , 2X  , I 1 , 3X  ) ) 

9940  FCRMAT  l 1 X , ' SCX = ' , F7 .4 , • SCY  = ',F7.4,'  SCZ=',F7.4) 

9945  FORMAT! 3! I 10, 10X)  ,15) 

995C  FORMAT!  • ISEED  = ' ,110) 

9955  FORMAT! IX,' THIS  SHCT  WAS  TERMINATED  BEFORE  IT', 

X • WAS  COMPLETED  • /,4X , ' ZST0P  = 1 THE  CURRENT', 

X • TIMESTEP  EXCEEDED  THE  L I M I T ' / , 1 4X , • AN Y SPUTS', 

X ' ARE  REMOVED  FROM  THE  SUMMARY') 

996C  FORMAT! IX, ' ***  THIS  SHOT  WAS  TERMINATED  3EFCRE  IT', 
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X • WAS  COMPLETED' / , 14X  , ' ISHUT=1  TIME  RAN  GUT'/, 

X 14X , ' ANY  SPUTS  ARE  REMOVED  FR CM  THE  SUMMARY ' ) 

9965  FORMAT ( 4X , ' IMPACT  PCINT  • , 1 5 ) 

9970  FCRMAT(4Xf,BX=' , F 6 .4  , 'BZ=' , F6 . 4 ) 

9975  FCFMAT(8X» 'NOT  IN  TARGET  AREA',/) 

CALL  CVFLOW 
CALL  SETIME 
READ  (5,9005)  TITLE! 

READ  ( 5 , 90 C 5 ) TITLE2 
REAC(5,91C5)  TRUN 
STOPER=T RUN *6 0.0-60.0 
FINISH=0.0 
I LATT  =0 
NCRC=4 
NCTT=10 
TESTE=2 .0 
E I ND  = 0 . 0 
CV8=0.774E-9 
6 ENR  C-Y  = 2 .0 
I SFIT=0 
LC  =0 

FM  = 1 . OE-1 0 
PFIV  = 0 .5 
CMC. 01 
S LCW  = 1 . OE—  14 
to  P I TE ( 6 , 960  1 ) 

CC  1 1=1,20 
I ECX 1 ( I ) =0 
I ECX2 ( I ) = 0 

1 CONTINUE 

DC  2 1=1, 60C 
LNUM ( I ) =0 
LCLT ( I ) =0 
RXK( I )=C.0 
R YK ( I ) =C  .0 
RZK(I)=C.O 
VX (I )=0  .0 
V Y ( I ) =0 . 0 
V2 ( I ) =0 .0 
PXU  )=0.0 
RY(I)=0.0 
R Z ( I ) =0 .0 
FX ( I ) =0  .0 
F Y ( I ) =0 . 0 

f zm-o.o 

PTE ( I ) =0.0 
PKE(I)=C.0 
PPE(I 1=0.0 
R X I ( I ) = 0.0 
R Y I ( I ) = 0.0 

2 RZI(I)=C.O 

REAO  ( 5, 9C10)  I H 1 

REAC(5, 9020  )IHS , IQR,SUR,LSS ,IX,IY,IZ,RBX,RBZ,CTI,RCE2, 
X LTYPE 

READ  (5,9030  I HB ULL , BUL , BMAS , I F E POT , PEX A , P EX E , XX XX 
READ (5, 9030  JIHTARG, TAR , TMAS 1 , 1 HT POT , EXA , EXB  , F LAT , TMAS 2 
X ,ABUN1  ,ABUN2,  ISOC-EN 

R c AC  (5,9040  I H2 »ROE A , ROES , ROEC , CCCN , A L PFA , RE , I PRCG 
5 R£AC(5,9080)cVR,SPX»SPZ ,NTT , NSR , NC , I DEEP , CO > , COY , I TEMP 
X ,7A ,CC ,ETFR 

R E AC ( 5 » 9945  ) I SEE D A , I SEED B , I S EE CC  , NSHOT S 
READ  (5,9915)  I R ANC-E , RM IDPT  , R INCR 
REAO  (5,9050  I RE  AD  , ( I FFM ( J ) , I F w ( J ) , J= 1 , I R E AD ) 
IF(EVR.LE.O.O)  GO  TO  9999 
IF ( ISOGEN. NE.O)  GC  TO  33 
7MAS2=63. 54 
AEUN1  = C .000 
AeLN2  = i .000 
33  CCNTINUE 
T AMP  = TA 
CCRF=CC 
I PASS=1 
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NCTT=NO 

IFl ILATT.NE  .0)  GQ  TO  10 
C 

C MAKE  HEADING  FCLLARITH 

C 

1 1-  E ( 1 ) = IHBU  LL ( 1 ) 

I h B ( 2 ) = IHBULL( 2) 

I FE (3 ) = IHTA  RG ( 1 ) 

IFE<4)=  IHTARG<2) 

IFB(5)=IHBPCT(1) 

IFB(6)=IH8PCT(2) 

IhT  ( 1)  = IHB  ( 2) 

I l-T  (2  ) = IHB  ( A ) 

I HT  <3 ) = IHQ ( 3 ) 

II-TU)-IHBU) 

I l-T  ( 5 ) = IHT  PCT  ( 1 ) 

I l-T  ( 6 ) = IHT  PCT  ( 2 ) 

C 

C MAKE  CONSTANTS 

C 

RACCEG=S7. 2957795 
RCE=SQRT(RCE2 ) 

EV=EVR*1C0C.0 

NS=NSR 

FAC  =0. 1 
RCEM=RCE-CTI 
CVR=0 .5  *FLAT 
CVE-1.602L6-L9 
CVM=1.6604E-27 
CVC=CVR*1.0E-10 
C VEC  = C V E/C  VC 
R XB  L=-R  C E 
RVBL=-J  .8 
R ZB  L=-R  C E 
IXP=(  IX  + D/2 
IVP=( IY+1J/2 
IZP=(IZ+l)/2 
PEMAS=BMAS*CVM 
HEMAS=0.5*P8MA$/CVE 
PTMAS1=TMAS1*CVM 
PTNAS2=TMAS2*CVM 
HTMAS 1=  0. 5*PTMAS 1/CVE 
HTMAS2=  0 * 5 *PTMA  S 2/C VE 
VTEST=—  SQRT  (2.0*ETHR*CVE/PTMAS2  ) 

V T S=0 . 8 *VTE  ST 
C 

C MAKE  REPULSIVE  POTENTIALS 

C 

FXA  = ALOG(-EXB*CV  EC  )+EXA 

PFXA=ALCG(-FEXB*CVED)+PEXA 

AC=ALOG(CVEC)+EXA 

PAC=AC-EXA+FEXA 

PTC^EXPtEXA+EXB^RCE) 

FRC=EXP ( FXA+EXB*RCE) 
FPTC=EXF(AC+EX8*RCE) 
PFPTC=EXP(PAC+PEX8*R0E  ) 
PFRC=EXP(PFXA+PEXB*ROE) 
PFTC=EXP(PEXA+PEXE*ROE  ) 

C 

C GENERATE  THE  LATTICE 

C 

IF  (I0R.EC.100)  GC  TO  6 

IF  (I0R.EC.110)  GC  TO  7 

IF  UOR.EQ.lll)  GC  TO  8 

6 CALL  KF 100 
GC  TO  9 

7 CALL  KF  110 
GC  TO  9 

E CALL  K F 1 1 1 
GC  TO  9 
9 CONTINUE 

IF(IOR.EQ.IOO)  THETA=1.0000 
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IFU0R.EQ.113)  TH  6T  A = C . 0000 

1 F ( IOR . EG  . 1 11 ) THETA= ( SORT ( 1.5)/3.0)/(SQRT(2. 01/2.0) 
IF( IS0GEN.EC.0)  GC  TO  12 
DC  12  1 = 1,1  READ 
J = I FFM ( I ) 

IMF(J)=IFM(I) 

12  CCNTINUE 

DC  65  1=2, LL, 3 
K = I + 1 
J = 1 + 2 

65  WR ITE (6 ,962C  ) I , IMF  ( I)  , RXI  ( I J ,R  Y I ( I ) , RZ I ( I ) , K , I MF  ( K ) , 
X RXI (K) ,RYI <K) ,RZI  (K) , J,IMF ( J)  , RXI  l J) ,RYI  ( J ),RZI  ( J) 
HR  ITE  { 6, 96G0 ) 

WRITE  { 6 , 993  J ) IRE AD 

WRITE  ( 6,9925)  ( I FFM l J ) , I FM ( J ) , J = 1 , 1 READ  1 
WRITE  (6,9600) 

WRITE  (6,9940)  SCX,SCY,SCZ 

LL2=2*LL 

LL2=3*LL 

R 1S=R1 

I LATT= 1 

N FAGE  = 1 

0X=RBX* SCX 

EZ=REZ*SCZ 

sex=BX 

SEZ=ez 

RXEM=(IX-1 )*SCX+RGE 
RYEM=(IY-1)#SCY+RCE 
RZBM=(IZ-1)*SCZ+RCE 
10  CCNTINUE 

CCZ=1. J-CCX*CCX-CCY*COY 
CCZ=A BS (CQZ  ) 

CC  Z = SQR  T ( CC  Z ) 

CCXI=-CCX 
CC Y I=— C CY 
CCZI=-CCZ 
I FT  = IT  I ME ( I I ) 

IRAN=2*!FT+1 

I F ( IRAN.LT.O)  I RAN=-IRAN 
EXMIN=BX 
BXMAX=EX  +DB  XMAX 
EZMIN=EZ 
EZMAX  = E Z+QB  ZMAX 
I F ( IOR. NE. 100)  GO  TO  25 
IF(SFX.EQ.J.J)  GO  TG  25 
eZMAX=BZMAX-10.0*$PX*SPCZ 
EZMIN  = BZMIN  + 10.0*SPX=!‘SPCZ 
25  IF ( IGR.NE. L 11)  GO  TC  30 
I F ( SPX . EQ  . 0 .0  ) GO  TO  30 
EZMAX=BZMAX-5.0*SFX#SPCZ 
EZMIN=EZMIN+5.0*SPX*SPCZ 
2C  CCNTINUE 

BX=(RBX+SPX  ) *SCX 
EY=(RYI (4J+C.001 ) *SCY 
EZ=(RBZ+SPZ  )*SCZ 
SCXR=  1.0/SCX 
SCZR  = 1 • C/SC  Z 

I F ( IOR.EQ.l  11)  SCZR=  1 .0/SSC  Z 

PREPARE  TC  WARM  LATTICE 

IFT=ITIME (II) 

I RAN=2J>  IFT  + 1 

IF( IRAN.LT.O)  I RAN=-IRAN 
WRITE  (6,9600) 

WRITE(6,97  01)  IFT,IP.AN,LL 
WRITE(6,9640)  E THR , VTE ST , VTS 
IF(TAMP.EQ.C.O)  GC  TO  40 
TAS=TAMP/SQRT (3.0  ) 

CCRN=1. C-CCPF 
T ASN  =TA  S*C  C RN 
TASF=TAS*CCPF 
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C 

C EACH  REPLICA  STARTS  HERE 

C 

AC  CCNTINUE 

IF(IPASS.GT.NSHOTS)  GO  TO  1000 
IF(TAMP.GT.C.O)  CALL  WARM 
50  CCNTINUE 
TIME=O.C 
NS=NSR 
N T = C 
N SS  = 0 
NCT =0 
I NCEX  = 0 
Z S TCP  =0 . 0 
S PX=  ( BX-SBX  )*SCXR 
SPZ=(6Z-S8Z)*SCZR 
ISPX=($PX+0.001)*100.0 
ISPZ=(SPZ+0.001)*1GC.0 
IFNUM1C0*ISPX+ISPZ 
I F NUMS ( IPASS)  = IPNUM 
R 1=R1S 
C 

C RECONSTRUCT  LATTICE 

C 

I S EC  ( I) =0 
CC  60  1=2, LL 
LCLTl I) =0 
ISECt I ) =0 
R X(I)=RXICI  ) 

RY(I)=RYI (I) 

R Z ( I ) =R  Z I ( I ) 

PTE(I)=0.0 
PPE ( I ) =C  .0 
VXdl  = 0.0 
VV ( I ) = o.c 
VZ(I)  = o.c 
6C  CCNTINUE 

IF (TAMP  .NE.C.O)  CALL  WARM 

R XS  = EX 

R Y S =BY 

RZS  = BZ 

R 1=R 1+F AC 

CALL  SETUP 

LCL  T ( 1 ) =0 

RXI ( 1 ) = RX ( 1 ) 

RYI ( 1)  = RY<  1 ) 

FZI(1)=RZ(I) 

PFE(1)=0.J 
VEL=SQRT( EV/HBMAS  ) 

VX(1)=VEL*CCX 
V Y ( 1 ) = V EL*C  CY 
VZ(1)=VEL*CCZ 
FCTI=DTI*CVC 
C T=FDT I / VE  L 
C 

C BEGIN  TIMESTEP  CYCLE 

C 

ICC  CTCD=DT/CVD 

HCTCD  = 0.5  *Q  TGD 

CTCM1=DT/PTMAS1 

CTCM2=CT/PTMAS2 

HCTCMC  l J =0 . 5*DT0M  I 

FCTCM(2)=0.5*DT0M2 

CTCMB=DT/P6NAS 

FCTCMB=C. 5*DT0MB 

2CC  CALL  STEP 

IF( INDEX. GT.JI  GC  TO  250 
C 

C FIRST  HALF  OF  TIMESTEP 

C 

21C  INDEX=  I \ 

1 = 1 
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I F ( LCUT ( I ) . EQ.2 ) GC  TO  215 
I F ( LCUT ( I ) . NE . 0 ) GG  TO  220 
R XK  ( I > =RX  ( I ) 

RYK  (I ) = RY ( I ) 

RZK(IJ=RZ(I) 

RX< I )=RX (I)  +DTOQ=MhDTCM8*FX  ( I J + VX  ( I ) ) 

RY ( I ) = RY ( I )+QTOD*(HDTCM8*FY ( I ) + VY  ( I ) ) 
RZ(I)=RZ(I) +DTQD# (hOTGMB*FZ(I)+VZ(I ) ) 

GC  TO  220 

215  RX ( I ) =RX ( I)+DTOD*VX(I  ) 

RY(I)=RY(I  )+DTOQ=VY(I  ) 

R Z ( I ) =R  Z ( I >+QTOQ*VZ(I  ) 

2 20  DC  240  1=2,  LL 

IF(LCUT(I ) .EQ.2)  GC  TO  235 
I F (LCUT ( I ) . GT. 0 ) GO  TO  240 
J= I NF ( I ) 

RXK(I )=RX( I ) 

R YK ( I ) = P Y ( I ) 

RZK ( I ) =RZ(  I ) 

PX(I  ) =R  X ( I )+DTOD*(HDTCM(J ) *FX  ( I ) + V X ( I ) ) 

RY  ( I ) =R  Y ( I )+OTOO*(HDTCM(J)*FY(  I )+VY(  I ) ) 

R Z ( I ) =R  Z ( I ) +DTOC* (HOTCM(J) *F  Z ( I )+VZ(I )) 

GC  TO  240 

235  RX(I)=RX< Ii+DTOD*VX(I  ) 

P Y ( I ) =R  Y ( I ) +DTOD* VY ( I ) 
PZ(I)=RZ(I)+DTOD*VZ(I  ) 

24C  CONTINUE 
GC  TO  2C0 
C 

C SECOND  HALF  GF  TIMESTEP 

C 

250  CONTINUE 
INCEX=0 
EMX=O.C 
T INE=TI PE+DT 
NT  = NT  + 1 
1 = 1 

I F ( LCUT ( I ) . NE . 0 ) GO  TO  260 
VSS=VX( I ) 

VX(  I)=VSS  + bCT0M8*FX(I  ) 

RX  ( I )=RXK<  I ) + ( VX ( I ) +VSS ) *HDTOQ 
\SS  = VY(  I ) 

VY(I)=VSS+HDTOMB*FY(I  ) 

RY(I)=RYK(I )+ (VY ( I )+VSS  )*HDTOD 
VSS=VZ< I ) 

VZ( I )=VSS  + bCTOMB*FZ  (I ) 

RZ(I)=RZK( I )+ ( VZ ( I ) +VSS )*HDTOD 
FX(li=u«0 
F V ( 11=0.0 
FZ ( 1)=0 .0 

PKE< I ) = VXl  I )*VXl I )+VY l I )*VY( I ) +VZi  I )*VZ (I) 
EI*AX  = PKE(  I ) 

K = I 

PKE(I)=PKE(I)*H8MA$ 

260  CC  290  1=2, LL 

IF (LCUT ( I ) .NE.O)  GC  TO  29  0 
J = I t*F  ( I ) 

VSS=VX(  I ) 

VX( I ) =VSS+HCTQM( J )*FX<  I ) 

RX(I)=RXK(I  ) + ( V X ( I ) +VSS )*HDTQD 
VSS=VY ( I ) 

V Y ( I )=VSS+HDTOM( J )*FY( I ) 

RY(I)=RYK(I  ) + ( V Y ( I ) + VSS  )*HDTQD 
\.$S  = VZ(  I ) 

VZ(I)=V$S+hDTCM(J  )*FZ<  I ) 

RZ(I)=RZK(I ) + ( V Z ( I ) +VS S ) = HDTOO 
FX ( I )=0 .0 
F Y ( I )=0 .0 
F Z ( I )=0.0 

FKE ( I ) =VX ( I )*VX(  I ) +VY ( I )*VY (I )+VZ( I )*VZ(  I ) 
IF(PKE(  I) .LT.EMAX  ) GO  TC  280 
EMX  = PKE  ( I ) 
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280  I F ( IMF  ( I ) • EC.  1 ) GC  TO  282 
P KE ( I ) = FTMA S2*PKE ( I ) 

GC  TO  290 

282  PKE(  I )=FTMAS1*PKE  (I } 

29C  CONTINUE 

TEST  FOR  PRINT  CR  TERMINATION 

CALL  GET  IME  ( IET) 

CALL  SETIME 

F IN  ISH=  IET  *2 . 6E-05  +FINISH 
IF< FINISH. GT.STOPER)  GO  TO  370 
1 F ( IMF ( K ) • E C • 1 ) GO  TO  293 
E7EST=HTMAS2*EMAX 
GC  TO  294 

293  E7EST=HTMAS1*EMAX 

294  CCNTINUE 

IF(ETEST.LE.TESTE)  GO  TO  38C 
NCT=NCT  +1 

IF(NCT.NE.NCTT)  GC  TO  360 
NCT  =0 

CC  305  1*1,  LL 

I F ( LCUT ( I ) .GT.O)  GO  TO  305 
P F E ( I )=C.O 
3 C 5 CONTINUE 

CALL  ENERGY 
3C7  CC  350  1*1  , LL 

I F ( LCUT ( I ) .GT.O;  GC  To  350 

IF(RY( I ) .LT  .RYBL  . AND.VY ( I ) . LT.O  .0  ) GO  TO  32C 
I F (RX  ( I ) . LT  .RXBL.CR.RX(I).GT.RXBM)  GO  TO  31C 
I F ( R Z ( I ).LT  .RZBL.CR.RZ(I) .GT.RZSM)  GO  TO  3 1 C 
IF(RY< I ) .LT.RY8M)  GO  TO  350 
3 1C  LCUT ( I ) =1 

RXF(I)=RX(I) 

R Y F ( I ) = RY ( I ) 

P2F< I ) = PZ(  I ) 

PTE( I)  = PKE(I)+PPE  (I) 

GC  TO  350 

32C  I F ( I.EQ.l)  GO  TO  310 

IFiPPEU)  .GT.O.  0001)  GO  TO  350 
PTE(I)  = PKE(  I ) +PPE (I) 

IF  <PTE<  I)  .LT.O.O)  GO  TO  350 
L CUT ( I) =2 
RX F ( I ) =RX ( I ) 

RYF  ( I ) = RY  ( I ) 

R 2 F < I ) = RZ  < I ) 

350  CONTINUE 

TEST  FOR  CUT  FUT  OF  CURRENT  STATUS 

3 6C  CT  = FDTI/SQRT(EMAX  ) 

IF(NT.GE.NTT)  GO  TO  380 

GO  TO  NEXT  TIMESTEP 

IFINS.NE.NT  ) GO  TC  100 
GC  TC  4C0 
37C  I SFUT  = 1 
GC  TO  4CD 
3EC  N S S = 1 

IF(NT.GE.NTT)  ZSTCP=1.0 

CALCULATE  POTENTIAL  ENERGY  AND  PRINT  STATUS 
OF  THE  CURRENT  T1 HE  STEP 

4GC  CONTINUE 
TFCT  = 0. 3 
PFCS=0. C 
FNEG=O.C 
CC  450  1=1, LL 

I F ( I SEQ ( I ) • EQ . 1 ) I SE 0 ( I ) =2 
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I F ( LCUT ( I ) . NE • 0 ) GC  TO  450 
PFE(I)=C.O 
PTE ( I ) =0.0 
45C  CONTINUE 

CALL  ENERGY 

555  PTE{i)=PKE (1) +PPE  (1) 

TFKE  = PKE(  1 ) 

6CC  CC  625  1=2, LL 

T FKE  = T FKE+PKE ( I ) 

IF ( PPE ( I ) .GT.O.O)  GO  TO  615 
FNEG=PNEG+PFE(  I) 

GC  TC  620 

615  PFCS  = PPCS+PPE( I ) 

62C  PTE(  I ) =PKE  ( D+PPS  (I) 

625  CONTINUE 

TE=TPKE+TFGT+BIND 

IF ( (NSS .NE.O)  .ANC.  ( ZSTOP . EQ .0 . 0 ) ) GO  TC  55J 
CALL  KPPINT 
WRITE  (6,9661) 

733  CC  750  1=1, LL 

IF  (PTE ( I ) .LT.O. 1 ) GO  TC  750 
IFdSEQd  ).EQ.O)  I SEQ ( I 1=1 
CX(I)=RX( I)-RXI (I  ) 

CY ( I )=RY(  I ) -RY I ( I ) 

CZ( I)=RZ( I)-RZI ( I ) 

WRITE (6, 9682)  I , 0 X ( I ) , D Y ( I ) , DZ ( I ) , VX ( I ) , VY ( ! ) , VZ ( I ) , 
XPKEd  ) , PPE  (I)  ,PTE(  I)  ,RX(I  ) ,RY  (I  ) ,RZd  ) , LCUT  (I  ) , ISEC(I) 
75C  CONTINUE 

WRITE  (6 ,9650  T?KE  , TPCT  ,TE  , PPOS  , FNEG 
IF< ISHUT.NE.O)  GO  TO  972 
IF(ZSTOP.NE.  0.0)  GO  TC  973 
WRITE(6 ,9630)  NPAGE 
NFAGE=NF AGE+1 
N S=NS+NC 

GO  TO  NEXT  TIMESTEP 

GC  TO  100 
95C  CONTINUE 

TERMINATE  REPLICA  ANC  TFEN 
PRINT  FINAL  STATUS  OF  CURRENT  PASS 

J = 1 
K = C 
I H =0 

CC  560  1=1, LL 

IF  (LCUT  ( I ) .LT.2  ) C-G  TC  96  0 
IF ( VY ( I ) .GT  .VTEST  ) GO  TC  96C 
K =K  + 1 

ISNLMt J)=I 
LNLM ( I )=LNUM(I )+l 
J=J  + 1 

5 6 C CONTINUE 
K K = K 

NSFUT(IPASS)=KK 
CALL  KPPINT 
IF(KK.EC.O)  GO  TO  972 
WRITE (6 ,9813)  EV,  ( I SNUM ( I ) , I = 1 , K K ) 

WPITE(6  ,9600 
WRITE (6, 9655) 

CO  963  11=1, KK 
I ~ I SNUM  (II  ) 

WRITE  (6  ,9800  IPNUM,  1 ,8UL  ,TAR,  I OR  ,SUR  (1),VX(I),VY(I),VZ 
X ( I ) ,PKE ( I)  ,PPE(  I ) ,PTE  (I ) ,EVR ,CCX ,CCZ 
966  CONTINUE 

CC  569  1=2, LL 

IF(LCUT(I).EQ.2  .ANC.  V Y ( I ) . L T . VTS  ) W R I TE  (6 , 930 0 ) 

X IPNUM  ,1  , SUL, TAR  , IOR,SUR(  1 ) ,VX(  I ) ,VY(I  ) , V Z ( I i , PK E 1 1 ) , 
X PPE( I J ,PTE( i)  ,EVR 
965  CONTINUE 
C 
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C 


DETERMINE  SPUTTERED  ATOM  ANGLE 

IF(KK.EQ.O)  GO  TO  572 
WRITE  ( 6 » 96C0 ) 

WRITE(6,98C5) 

CC  570  I = L * KK 
L = I SNUM ( I ) 

RAC  = ATAN(SGRT(VX(L)^#2  + VZ(LJ^*2)/(-VY(Li  ) ) 

ANGLE=PAC*RADDEG 

J = I F I X ( (ANGLE /RINCR)+1.0) 

I F ( IMF  ( L)  .EC.l  ) I BOX  1 ( J ) = I BOX  1 ( J ) +1 
1 F ( IMF ( L ) • E C.2 ) ieGX2( J ) = I 80X2  < J I +1 
WRITE  (6,9815)  L, ANGLE 
97C  CONTINUE 

1F< (ISHUT.EQ.l)  .CR.  (ZSTOP .EG  . 1 .0) ) GO  TC  572 
WRITE<  6,9630  NPAGE 
NPAGE=NFAGE+1 
GC  TO  950 

972  IF( ISHUT.EQ  .0)  GO  TO  973 
WRITE (6 ,9965)  IPNUV 
WRITE(6,960C) 

WRITE  ( 6 ,9960 
N SPOT ( I F A S S ) =0 
I F A SS  = I FASS-1 
NSHCT$=N3HC7S-1 
GC  TO  975 

972  IF(ZSTCF.EC.O)  GO  TO  974 
WRITE(6, 99  6 5 ) IPNUM 
WRITE (6  ic  ) 

W R I TE ( 6 ) 

N S FLT  ( =0 

IFASS 

NSI-C  -1 

GC  T 

974  WRIT. ,^d2C)  EV 
97  5 WRITE(6  ,9630  NPAGE 
NFAC-E=NPAGE  + 1 

DEFINE  NEW  REPLICA  OR  TERMINATE 

95  C 1PASS  = I FASS  + l 

IF ( ISHUT.NE  .0)  GO  TO  LOGO 
EZ=BZ+SPDZ 

IF (EZ.LE.SZMAX)  GO  TO  40 

ex=ex+spox 

IF(BX.GT.aXMAX)  GO  TO  1000 
I F ( I0R.EQ.1C0)  BZMIN=BZMIN  + SPDZ 
I F ( ICR. EC. Ill)  BZMIN=BZMIN+SPDZ/2.0 
E Z = E ZM I N 

IF ( I0R.EQ.1C0)  BZMAX=eZMAX-SPDZ 
IF ( IOR. EC. Ill ) BZMAX=EZMAX-SPDZ/2.0 
C-C  TO  4 0 

PRINT  SUMMARY  OF  ALL  SHCTS 

100C  CONTINUE 

CALL  KPRINT 
WRITE  (6,9600) 

N R AN  = I P ASS  - 1 

WRITE  ( 6 ,9860  NR  A N , EVR 

WFITE(6,987C)  ((  IPNLMS<  I)  ,NSPUT(  I))  ,1=1,  NR  AM 

KSLM=0 

K2SLM=0 

DC  1005  1 = 1 ,NRAN 
KSLN=K5LM+NSPUT< I ) 

K2SLM=K2SUN+NSpUT  (I )**2 
IOC 5 CONTINUE 
FRAN=NRAN 
FR  ANR  = 1 .0/ F RAN 
AVE  = FLOAT  ( K SUM  ) =*F  R ANR 
VAR=FLOAT U2SUM) *Fh ANR-AVE**2 
STC=VAR  *FRANR 
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10C7 

10C8 

IOCS 


1011 

10  12 

10  12 
1014 


10  15 


1020 


1025 

102C 

104C 


STC=SQR  T ( S TC ) 

WRITE  (6,9600) 

WRITEU, 9875)  ETHR  , VTEST , KSUM 
WRITEU  ,98  50)  AVE  , VAR, IMRAN,  STD 
WRITE  (6,9600) 

WRITE  (6,9865) 

WRITEU  ,9870  ( I , LNUM ( I ) » I = 1 » L L ) 

WRITE  (6,9600) 

WRITE  (6,9920)  IREAO 

WRITE  (6,9925)  ( I FFM ( J ) , I FM ( J ) , J = 1 , 1 READ  ) 

IF  ( I0R.NE.100)  GO  TC  1007 
WRITE(6  ,9950  ISEEDA 
IF(IOR.NE.llO)  GO  TO  1008 
WRITE  (6, 9950  ISEEC8 
I F ( I0R.NE.111)  GO  TO  1009 
WF  ITE  (6  ,9950  ISEECC 
CONTINUE 

WRITEtd, 9630  NPAGE 
NPAGE=NPAGE+1 
CALL  KPNT 

FOLD  FINAL  SUMMARY 

IF  (I0R.EQ.100)  GC  TO  1011 

IF  (I0R.EG.110)  GC  TD  1012 

IF  ( I OR  . EG  . 11 1 ) GC  TO  1013 

CALL  KFCLDA 

GC  TO  1 C 14 

CALL  KFCL08 

GC  TO  1014 

CALL  KF  CLDC 

WRITE  (6,9600) 

WRITE  (6,9840) 

WRITE  (6,9845) 

WRITEU  ,9870  ( I , LNUM ( I ) , I = 1 , LL ) 

WRITE  ( 6 , 9 6 CO ) 

DETERMINE  RING  SPUTTERING  PROBABILITY 

CALL  KPPCB 
WRITE (6  ,9890) 

WRITEU  ,9880 
CC  1020  1 = 1 ,KR I NG , 3 
J = I + 1 
K = I +2 


WRITEU,  98  85)1, DNN(I)  , K ATOM  (I),KSPUT(I),PROE(I),J, 

X DNN(  J j , K A T 0 R ( J ) , KSPUT ( J ) , PROS ( J ) , K , ONN ( K ) ,K ATOM  < K ) , 
X KSPUT (K) , FROB (K  ) 

CCN.TINUE 
WRITE  l 6 , 96C0 ) 

WRITE  (6,9620)  NPAGE 
NFAC-E=NPAGE+1 
CALL  KP NT 
WRITE  ( 6 , 9 6 CO ) 

CRAW  ISOTOPE  HISTOGRAMS 

CC  1030  1=  1 , 1 RANGE 
FRECKI  ) = I B CX 1 ( I ) 

F F EC2 ( I ) = I B CX2 ( I ) 

CONTINUE 
RANGE ( 1 )=RM  IDPT 
CC  1040  1=2 , IRANGE 
J = I - 1 

R A N C-  E ( I )=  RANGE  (J  )+RINCR 

CONTINUE 

IN  =-IRANGE 

CALL  H I STC ( IN,FREC1,RANGE,TITLE1  ) 

WRITE  ( 6 , 96C0) 

WRITE  (6,9895) 

WRITE  (6,9900)  RINCR 

IF  ( IRANGE .EQ. 9 ) WRITE  (6,9925) 
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IF  ( IRANGE  .EQ.I8 ) WR I TE ( 6 , 9920  ) 

WRITE  (6,9910)  ( I BCX1 ( J > ,J= 1 , IRANGE > 
WRITE  ( 6 , 9 6C0 ) 

WFITE  ( 6 , 96C0 ) 

CALL  HI STQ { IN , F REC2 , R ANGE , T ITL E 2 ) 
WRITE  (6,9600) 

WRITE  (6,9895) 

WRITE  (6,9900)  R I NCR 

IF  ( IRANGE. EQ. 9)  WRITE  (6,9925) 

IF  ( IRANGE. EQ. 18)  WR I T E ( 6 , 9 920  ) 

WRITE  (6,9910)  ( I 80X2 ( J ) ,J=  1 , IRANGE ) 
WRITE  ( 6 , 96C0 ) 

WFITE  ( 6,9620)  N PAGE 
NFAGE=NF AGE+1 
IF(ISHUT.EC.O)  GO  TG  5 
9999  STCP 
ENC 
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SLEFOUTINE  SETUP 


THIS  S/R  FINCS  THE  STARTING  POINT  FCR  EACH 
IMPACT  POINT  IN  THE  TARGET  AREA 


CCMMCN/C0M1/RX( 600 ,RY(600),RZ(60C),LCUT(60C),LL,LD 
CCMMCN/CCM2/RCE,RCE2,RCEM,AC,PAC,FPTC,PTC,PFPTC,FPTC, 
X FN,PFIV,TPCT 

CCMMCN/C0M6/C0XI , COY  I ,CCZI , RXS,RYS,RZS,FAC 
CCMMCN/CCM7/R1, LSSt SPX, SPZ,COX,CCYfCCZt YLAX (2  3 ) » I LI 
CCMMCN/COMI 1/RXI  ( 600  ,RYI  (600)  ,RZI  (600) 

1 = 1 

100  R 1=R1-FAC 

IF (Rl.LT.O.C)  GO  TO  300 
R 1X=R1*CCXI +RXS 
R 1 Y =R  l*CCY  I+RYS 
F 1Z=R1*CCZI+RZS 
X MN  = R 1 X- ROE 

Y N IN  = R1Y-R0E 
ZM  IN  = R1Z-P0E 
X M AX  = R1X  + RGE 

Y M A X = R1Y+R0E 
ZNAX  = P1Z+P0E 
CC  195  J = 2 * LL 

IF  (LCUT  ( J ) .C-T.O  ) GC  TO  195 
IF(RYfJ).GE.YMAX)  GC  TO  195 
IF (RY( J ) .LE.YMIN)  GO  TO  195 
IF (PX(J ) .GT.XMAX)  GO  TO  195 
I F < R X < J J .LE.XMIN)  CC  TO  195 
IFtRZIJ ) .GE.ZMAX)  GO  TQ  195 
I F ( R Z ( J ).L5.ZMIN)  GO  TO  195 
C F X =RX  ( J )-  R IX 
CRY=RY(  J )-RlY 
CRZ=RZ( J ) -R 1Z 
0 I S T =0 R X **  2 +D RY**  2 + DR  Z * *2 
I F ( C I ST .GE.  RQE2 ) GO  TO  195 
GC  TC  200 
195  CONTINUE 
C-C  TO  ICO 
2 0 C CONTINUE 

R2SQ=  (R  X I ( J )-RXS  ) **2+ (RYI ( J ) -RYS ) **2+ (RZ I ( J )-RZS) **2 
R 2 =SCRT (R2SC) 

R IF  2=  ( P 1X-RXS)* ( RXI ( J )-RXS) +(RLY-PYS )*( RYI  ( J ) -RYS  ) + 

X (R1Z-RZS)*(RZI (J)-RZS) 

ALFA=R1F2/(P1*R2) 

ALFA2=ALFA*ALFA 

R2=R2*ALFA+SQRT(R2SQ*ALFA2-R2SQ+RCE2> 

RX ( 1 )=R3*CCXI+RXS 
R Y ( 1 ) =R  2*C  C Y I +R  Y S 
RZ(1)=R2*CCZI+RZS 
RETURN 

3CC  RX ( 1 ) =RXS 
PY ( 1 )=RYS 
R Z ( 1 ) =R  Z S 

97CC  FORMAT ( 4 1 H THIS  QIC  NCT  UOCATE  A PROPER  START  POTNT) 
WRITE  ( 6,9700) 

RETURN 

ENC 


SLEROUTINE  WARM 


THIS  IS  A MCOEU  3 WARMER 

EACH  ATOM  IS  DISPCACED  AN  AMOUNT  CETERMINEC 
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BY  ITS  THERMAL  AMPLITUDE 


CCMMON/COM1/RX(60J)  » R Y ( 600 ) , RZ ( 600  ) , LOUT ( 60 C ) , LL , LD 
CCMMCN/C0M1 3/  RAN  (1800) , I R AN , T A SN , TASF , L L2  , LL3 
CCMMGN/C0M14/  AXT  (6C0)  , AYT(600)  , AZT(600) 

A XT ( 1 ) = 0 
AYT (1 )=C 
A 2T ( 1 ) =C 

CALL  NCFMAL  (IRAN, RAN, LL3) 

CC  1100  1=2, LL 
A X T ( I ) =RAN ( I) 

A YT ( I ) =P  A N ( I + LL ) 

AZT ( I ) =R AN ( I+LL2  ) 

RX ( I ) =RX ( I ) +TASN* AXT ( I ) +T AS  F*AX  T ( I- 1 ) 

R Y ( I ) =R Y ( I )+TASN*AYT( I ) +T ASF*AYT ( I — 1 ) 

RZ  (I )=RZ( I )+TASN*AZT( I ) +TAS F*AZT ( 1-1 1 
1100  CONTINUE 
R ETLRN 
t N C 


SLERCUTINE  KF100 


LATTICE  GENERATOR  FCR  A FCC  (100)  ORIENTATION 
LATTICE  IS  DEVELOPED  IN  THE  X DIRECTION,  THEN 
ThE  Z DIRECTION,  THEN  FINALLY  THE  Y DIRECTION 


CCMMCN/CCM1/RX(60C) ,RY (600 ),RZ(60C),LCUT(60C),LL,LC 
CCNMCN/CCM2/RCE, RCE2 , ROEM , AC , P AC , P PTC, P TC, PFP TC ,F PTC , 

X FM,PFIV,TPCT 

C C MNCN/ COM4 /IX, I Y ,IZ  ,IXP,IYP,SCX,SCY,SCZ,IDEEP, 

X PEX,ReZ 

CCMMCN/CGM7/R1,LSS,SPX,  SPZ,CCX  , C C Y ,COZ  , Y L A X (20  , I LL 
CCMMCN/C0M1 1/RXI (600),RYI(600),PZI(600) 

CCMMCN/CQM23/  SP D X , SPDZ  ,08X N AX , C E ZMAX , SSC Z 
CCMMCN/CCM25/IMF(600) , ISOGEN 
CCMMCN/CCM3C/  ISEECA,  ISEEOB,  ISEECC 
DIMENSION  R (600 ) 

9001  FC  FMAT ( ' I S EED= ' , 1 10// ) 

9690  FC?MAT( 25X , 'SUMMARY  OF  ATOMS  FCC  (100)  ORIENTATION* 

X / , 3 5X  , ' MA  SS  1 IS  THE  HEAV  I ER /L E A ST  ABUNDANT  ISOTOPE' 
X / / , 3 ( ' ATOM  MASS  POSITION  '),/) 

L C = C 
N =600 

CC  3 1=2, N 
I N F ( I ) = 2 
2 CONTINUE 

I S EED= I SEED  A 

WRITE  (6,9001)  I5EED 

WRITE  (6,9690) 

IF  (ISOGEN. EQ.O)  GO  TO  4 

GENERATE  RANCOM  MASSES 

CALL  GGU3 ( I SE ED , N , R ) 

CC  4 1=2, N 

IF  (R(I)  .IE.C.3C9)  IMF(I)  = 1 

4 CONTINUE 

CC  5 1 = 3,  IY 

5 YLAX( I ) =0.0 

DETERMINE  ORIENTATION  PARAMETERS 

SCX=1.0 
S C Y = 1 . 0 
SCZ=1.0 
SFCX  = 0. 1 
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SFCZ=0.  1 
C EXHAX=G .501 
C EZNAX= 1.001 
YLAX<1)=-0.12 
YLAX(2) =-0.03 

GENERATE  THE  LATTICE 

N = 2 
J T =0 

Y =-SCY 

CC  60  J = l,  IY 

Y = Y + SCY 
K T =C 

Z = -SCZ 

CC  59  K =1 » I Z 
Z=Z+SCZ 
IT  = 0 
X=-SCX 

CC  58  1 = 1, IX 

x=>+$cx 

I TT  = IT+ JT  + KT 

I F ( ITT- C ITT/2)*2  ) 57,30,57 
2C  R > I ( M ) = X 

R Y I ( M ) = Y + YL AX ( J ) 

R Z I ( M ) = Z 
N = N + 1 

57  CONTINUE 
IT  = I T+ 1 
56  CCNTINLE 
K T = KT  + 1 
55  CONTINUE 

IF(JT.EC.O)  I L L = M - 1 
JT  = JT  + 1 

IF(JT.EC.I DEEP ) L C = M - 1 
6C  CONTINUE 
L L = H— 1 

IF(LD.EC.O)  LD=LL-1 
IF(CCY.EQ.O)  GO  TC  120 
GC  TC  ( 70,80,90)  ,LSS 
7 0 R1=PCE/A0S (COY) 

GC  TO  ICO 
ac  I S F X = RS  X 
I5FZ=R8Z 

N=( IX/2 )*ISFZ+ISPX/2+3 
IF < IX-( IX/2  )*2)  82,85,83 
62  N=N+ISPX/2 
65  LCLT ( N ) = 1 
P 1 =2 . 0 
LCLT ( 2 ) =-I 
GC  TC  100 
50  V = 2 

R X I (M)=R8X*SCX 

RY  I (M )=-SCY+YLAX ( 1) 

RZI <M)=(R8Z+1.0)*SCZ 
R 1= (ROE  + SCY  )/A35 ( CO Y ) 

ICC  CONTINUE 
11C  PETURN 
12C  R 1=  1 . 5 
RETURN 
ENC 


SLEPOUTINE  KF110 


LATTICE  GENERATOR  FOR  A FCC  (110)  CPIENTATICN 
LATTICE  IS  DEVELOPED  IN  THE  X CIRECTICN,  THEN 
THE  Z DIRECTION,  THEN  FINALLY  THE  Y DIRECTION 
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CCMMCN/CGM1/RX(60C ) ,RY(600) ,RZ ( 6CC ) , LCUT ( 6 0 C ) t LL , LD 
CCMMCN/CCM2/RCE»RCE2,R0EM»AC»PAC,PPTC,PTC»PFPTC,FPTC» 

X FN,PFIV,TPCT 

CCMMCN/C0M4/IX, I Y , I Z , I X P , I Y P , SC  X , SCY  , SC  Z , 1 0 EE  P , 

X PEXfRBZ 

CCMMCN/C0M7/R1,LSS, SPX , SPZ , COX  ,CC Y ,C0Z , YLAX  (20) , I LL 
CCMMCN/COM11/RXI (600)  , R Y I ( 6 00  ) , R Z I ( 6 0 0 ) 

CCMMCN/CGM23/  SPOX , SPDZ , DBX MAX  , DBZ VAX , SSC Z 
CCMMON/COM25/IMF( 6CC) , ISOGEN 
CCMMCN/ CGM30/  ISEEDA,  ISEEDB,  ISE60C 
C IMENS I CN  R (600 ) 

90G1  FCPMATC  IS EED=  • , 1 10// ) 

969C  FCPMATl 25X, 'SUMMARY  OF  ATOMS  FCC  (110)  CR  IENT  AT  I GN  ' 

X /»35X  , ' MASS  I IS  THE  HEAV I ER /L EAST  ABUNDANT  ISOTCPE' 
X / / » 3 ( ' ATOM  MASS  POSITICN  •),/) 

LC  = C 
N=6C0 

CC  3 1=2, N 
IMF(I)=2 

3 CONTINUE 
I S E ED= I SE5C  2 
V>  R I TE  ( 6 , 90C I ) ISEED 
fcPITE  (6,9690) 

IF  ( ISOGEN. EQ.O)  GC  TO  4 

GENERATE  RANCOM  MASSES 

CALL  GGL B ( ISEED, N, R ) 

CC  4 1=2, N 

I F ( R ( I ) .LE.C.309)  IMF ( I ) = 1 

4 CONTINUE 
DC  5 1=3, IY 

5 YLAX(I)=O.C 

DETERMINE  ORIENTATION  PARAMETERS 

PC  = 1.0/SQRT  (2.0) 

SCX=RQ 
S C Y =R  0 
SCZ=I.O 
SFCX  = 0. 2*SC  X 
S PC  Z=0 . 2 *SC  Z 
CEXMAX=SCX+0.01 
C EZMAX=SCZ  + C.Ol 

Y LAX ( I ) =-0 . 15 

Y L A X ( 2 ) =-0.03 

GENERATE  THE  LATTICE 

M = 2 
J T =0 
Y=-SCY 

CC  60  J = I,  IY 
Y=Y+SCY 
KT  = C 
2 = -SC  Z 

CC  59  K = I , I Z 
Z=2+SCZ 
I T =0 
X=-SCX 

CC  58  1 = 1, IX 

x = x +sr  x 

I F 4 1 T— ( IT/2)  *2)  21,11,21 

11  I F ( JT- ( JT/2  )*2  ) 57,12,57 

12  IF  (KT-(KT/2  )*2)  57,3C,57 

21  IF(jT-( jT/2)*2)  22,57,22 

22  IF(KT-(KT/2  )*2  ) 3C, 57, 3C 
20  R X I ( M ) = X 

RYI (M)=Y+YLAX( J) 

RZI(M)=Z 
► a N +1 
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57  CCMINUE 
IT  = I T+  I 

58  CONTINUE 
KT=KT+1 

59  CCNTINUE 

JT  = JT  + 1 

I F ( JT.EC.IOEEP)  LC=M-1 

60  CCMINUE 
LL=M-1 

IF(LD.EC.O)  LD=LL- 1 
I F (COY . Efl  .0  ) GO  TO  120 
GC  TO  ( 70,80  ,9  3)  ,LSS 
70  R1=RCE/ABS(C0Y) 

GC  TO  100 
80  I$FX=SPX 
ISPZ=SPZ 

N = { IX+1  )/2*ISPZ/2+ISPX  + 2 
LCLT(N) =1 
R 1 = RCE/A8S  (COY) 

LCLT ( 2 ) =-l 
GC  TC  100 
90  M = 2 

RXI (M)=(SPX+1.0)*SCX 
RYI(M)=-SCY+YLAX(  1) 
RZI(M)=(SPZ+1.0)*SCZ 
R1=<RGE+SCY)/ASS(C0Y) 

ICC  CCNTINUE 
11C  RETURN 
12C  R 1 = 1 . 5 

RETURN 

ENC 


SLERCUTINE  KFlll 


LATTICE  GENERATCR  FCR  A FCC  (111)  ORIENTATICN 
LATTICE  IS  DEVELOPED  IN  TPE  X DIRECTION,  THEN 
THE  Z DIRECTION,  THEN  FINALLY  THE  Y DIRECTION 

COMMON/ COM 1/RX( 60 C ) ,RY(600),RZ(6CC),LCUT(60C),LL,LD 
CCMMCN/CCM2/RCE,RCE2,RCEM,AC,PAC,PPTC,PTC,PFPTC,FPTC, 

X F M ,PFIV»TPCT 

CCNMCN/C0M4/IX, I Y , IZ , I X F , I Y P , SC  X , SCY  , SC Z , I C EE P , 

X REX,R ez 

CCMMCN/CCM7/R1,LSS,SPX,SPZ,CCX,CCY,CCZ,YLAX(2C),ILL 
CCMMCN/COM11/RXI (600),RYI(600),RZI(60J) 

CCMMCN/CCM23/  SPDX,$PDZ,DBXMAX,CeZNAX,SSCZ 
CCMMCN/COM25/ IMF ( 6CC) , ISOGEN 
CCMMCN/CCM33/  ISEECA,  ISEED8,  ISEECC 
DIMENSICN  R ( 600 ) 

9001  FORMAT ( ' I S EED= • , 1 10/ / ) 

9690  FCRMAT(35X, 'SUMMARY  OF  ATOMS  FCC  ( 1 L 1 ) ORIENTATION' 

X / ,3  5X  , « MASS  1 IS  THE  HE AV  I ER / L E AST  ABUNDANT  ISOTOPE' 
X / / , 3 ( ' ATOM  MASS  POSITION  '),/) 

LC  =C 
N =600 

DC  3 1=2, N 
IMF ( I ) = 2 
3 CCNTINUE 
ISEEC=I SEECC 
WRITE  (6  ,9001)  I SEED 
WRITE  (6,9690) 

IF  ( ISOGEN. EQ.O)  GC  TO  4 

GENERATE  RANDOM  MASSES 

CALL  GGL  B ( I SEED , N v R ) 

CC  4 1=2, N 


non  oon 


IF(R(I)  .LE.C.309)  IMF ( I ) = 1 

4 CCNTINUE 

DC  5 1 = 2, IY 

5 YLAX(I>=0.0 


DETERMINE  ORIENTATION  PARAMETERS 


SCX=1.0/SQRT(2.0! 

SCY=2.0/SQRT<3.0) 

SCZ=SQRT(I .5) 
SSC7=2.0*SCZ/3.0 
SPCX=0. 2*SCX 
SFCZ=0. 2*SSCZ 
Ce>MAX=SCX+C.01 
CEZMAX=SSCZ+0.01 

Y L AX ( I) =-0.C5 

Y L A X ( 2 ) =-0 .Cl 

GENERATE  ThE  LATTICE 

N=2 
JT  = C 
Y=-SCY 

DC  60  J=l,  IY 

Y = Y + SCY 

J 1 S = JT+ JT/ 3 
Z=-SCZ 

CC~59  K = 1 , IZ 
Z=Z+SCZ 
I 7 = C 
X=-SCX 

DC  58  1=1, IX 

x=x+scx 

IN=IT+JTS+KT 

IF  ( I A- < IN/2  )*2)  57,30,57 
2C  PXI(M)=X 

RYI  ( M ) = Y + Y L AX ( J ) 

IF(JT-r3*(  JT/3)  ) A 1 ,45 , 41 

41  JTT= JT 

42  0 1 T = J TT-3 
IF(JTT)  43,45,42 
JTT  =JTT  +3 

Z F = JTT 

RZI (M)=Z+ZP$SSCZ 
GC  TO  5 C 
45  P Z I ( M ) = Z 
5C  N=N+1 
57  CCNTINUE 
I T = I T+  1 
5 £ CCNTINUE 
KT  =KT  + 1 
55  CCNTINUE 
JT= JT+1 

IF(JT.EC.IDEEP)  LC=M-1 
6C  CCNTINUE 
LL=M-1 

IF(LD.EC.O)  LD=LL-I 
IF(COY.EQ.O)  GO  TO  120 
GC  TC  ( 70,80,90) ,LSS 
7 C R 1 = RCE/ ABS (COY) 

LCUT  C 2 ) =- 1 
GC  TO  100 

ec  i s fx=spx 

ISFZ=SPZ 

N = I X*  ( I SPZ/2J  + ISPX  + 2 
LCUT ( N ) = 1 
R1=RCE/A6S(C0Y) 

GC  TC  100 
GC  M = 2 

RX  I (M)  = SPX*SCX 
RYI  (M)=~SCY  + YLAX(  1) 
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RZI (M)=(SPZ+1.0)*SCZ 
R 1= (ROE+SCY  )/ABS ( COY ) 
ICC  CCMINUE 
11C  RETURN 
1 2 C P 1=  1 . 5 
RETURN 
ENC 


SL6RCUT  INE  ENERGY 


THIS  S/R  CALCULATES  THE  MUTUAL  POTENTIAL  ENERGIES 
THIS  S/R  ALLOWS  ATTRACTIVE  POTENTIALS 
THIS  IS  THE  MUNDKUR  HIGH  SPEED  VERSION 


CCMMCN/CCMl/RX(60C},RY<600)tRZ(60C  ) ,LCUT (60C  1 ,LL,LC 
CCMMCN/CCM2/RCE, RCE2,R0EM,ACtPAC  tPPTC,?TC,  PFPTC  ,FPTC  , 

X FM  » PF IV  » T FCT 

CCMMCN/COM3/EXA ,EXB,FXAfPEXA,P5XE,PFXAfTFERN 
COMMON/ COM5/FX ( 60C 1 ,FY( 600)  , FZ < 60C ) , PPE ( 6CC  ) 
CCPMCN/C0M22/R0EA,RCEBtRCEC  » R0EC2fCP3tCPl» CF2»CP3,CF0 , 
X CFltCF2fCGDl»CGC2tCGBl»CGE2tCGF»IPRQG 
1 = 1 

I F ( LCUT ( I ) . NE . 0 ) GC  TC  200 
I F = I + 1 

INITIAL  POTENTIAL  ENERC-Y  INTERACTION 
MAXIMUM  DISTANCE  BOUNDARIES 

X N AX  = RX ( I ) +ROE 
YMAX  = RY  ( I ) +ROE 
ZM  AX  = RZ1D+RCE 
XMIN  = RX ( I )-ROE 
YMIN  = RY  (I  1-ROE 
ZM  IN  = R Z ( I l-ROE 

TEST  TO  DETERMINE  IF  BULLET  AND  ATOMS 
ARE  CLOSE  ENOLGH  TO  INTERACT 

CC  195  J=IP,LL 
I F ( LCUT  ( J ) .C-T.D)  GC  TO  195 
IF  (RY< J I .GE  .YMAX ) GO  TG  195 
IF  (RY ( J l.LE.YMIN)  GO  TO  195 
IF <RX< J ) .GT.XMAX ) GC  TO  195 
IF  { *5  x C J ) . LE  .XMIN  1 GO  TO  195 
IF(RZ(J).GE.ZMAX)  GO  TO  195 
IF  (RZ(J  ) . LE.ZMI N)  GO  TO  195 

TEST  TO  DETERMINE  IF  THE  ATOMS  ARE  INSIDE  A 
CIRCLE  OF  RADIUS  ROE  FROM  THE  BULLET 

C RX  =RX ( J ) -R  X ( I) 

C R Y =RY ( J 1 -R  Y ( I ) 

CRZ=RZ(J)-RZ(  I) 

D IST=CRX**2+0RY**2+DRZ**2 
I F(CIST.GE.R0E2)  GC  TO  195 
CIST  = SOPT (CIST) 

CALCULATION  OF  THE  POTENTIAL  ENERGY 

PC  T = EXP (PEXA+PEXB*OIST l-PPTC 
IF(PCT.LE.C.OCOOl)  POT=0.0 
TPCT=TPCT+PCT 
SAVE  = PF  IV*PCT 
PPE<  I ) = PPE { Il+SAVE 
PPE(J)=PPE( Jl+SAVE 
195  CONTINUE 
2CC  CONTINUE 
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CC  300  I = 2 , LD 

IF1LCUT (I) .GT.O)  GC  TO  300 
IF=I+l 

XMAX  = R X ( I J + ROE 
X M I N = RX ( I )-R06 
Y MX  = RYIIJ+ROE 
YMIN  = RY  ( I )-R0E 
ZMAX  = RZ  ( I ) + R0E 
ZMIN  = RZ  ( I J-ROE 

TEST  TO  DETERMINE  IF  THE  ATOMS 
ARE  CLOSE  ENOUGH  TO  INTERACT 

CC  295  J=IP,LL 
IF(LCUTIJ)  .GT.O)  GC  TO  295 
IF ( RY ( J ) . GE  .YMAX ) GO  T0  295 
IF <PY(J ) .LE.YMIN)  GC  TO  295 
IF (RX( J ) .GE  .XMAX  ) GC  TO  295 
IF { R X ( J ) .LE.XMIN)  GC  TO  295 
IF <RZ(J) .G£  .ZMAX)  GC  TO  295 
I F ( RZ < J ) • LE  .ZM I N ) GO  TO  295 
CRX=RX ( J ) -R X ( I) 

DRY=RY l J )-RY(  I) 

CRZ=RZ( J)-RZ(I) 

D IST=0RX**2+CRY**2+CRZ**2 
IFICIST .GE.ROE2 > GC  TO  295 
D IST  = SQRT ( D 1ST ) 

CALCULAT I CN  CF  THE  FOTENTIAL  ENERGY 

FCT  = EXP IEXA  + EX6*0  1ST  )-PTC 
IF(POT .LE.O.OOCOI ) PGT=0.0 
28C  TFCT=TPCT+PCT 
SAVE=PF I V*FCT 
PPE  ( I ) = PPE  ( D+SAVE 
PFE (J)=PPE( J) + SA  V E 
295  CONTINUE 
3CC  CONTINUE 

DC  313  1=1, LL 

I F ( LCUT ( I ) . EQ.O ) GC  TC  310 
T FCT  = TP CT  + P PE ( I ) 

3 1C  CONTINUE 
RETURN 
END 


SLERCUTINE  STEP 


THIS  S/R  CALCULATES  THE  FORCES  BETWEEN  ATOMS 
ATTRACTIVE  FORCES  ARE  ALLOWED  BETWEEN  ATOMS 
THIS  IS  THE  MLNDKUR  hIGH  SPEED  VERSION 


CCMMCN/CCM1/RX(60C)  ,RY(600)  ,RZ ( fcOC ) , LCUT ( 60 C ) ,LL,LD 
COMMON/ CCM2/RCE, RCE2 ,RCEM, AC, PAC ,F FTC, PTC, PFPTC ,FPTC  , 

X FM  , PF  I V , T PCT 

CCMMCN/C0M3/EXA  , EXB , FXA , P EX A , P t X E , P FXA , T HER X 
CCMMCN/C0M5/FX160C)  ,FY(600)  , F Z ( 60C ) , PPE ( 600  ) 

COMMON/ COM 2 2/ ROE A ,RCEB , ROEC , RGEC 2 , C PO ,CP1,CF2,CP3,CF0, 
X CFl,CF2tCC-Dl,CGC2,CGBl,CGE2,CGFf  IPRQG 
1 = 1 

I F (LCUT ( I ) . Nt . 0 ) GO  TC  200 
IF  = 1 + 1 

INTIAL  FORCE  DISTANCE  BOUNDARIES 
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ZHAX  = RZ  ( I )+ROE 
XMN  = RX  (I  )-R0£ 

KCIN  = RY ( I J-ROE 
ZMN  = RZdJ-ROE 

TEST  TO  DETERMINE  IF  BULLET  AND  ATGMS  ARE 
CLOSE  ENOUGH  TO  INTERACT  AND  WHICH  FORCE 
FUNCTION  TC  USE 

DC  195  J=IP,LL 
IF(LCUT(J)  .GT.O)  GC  TO  195 
IF(RY(J).GE.YMAX)  GO  TC  195 
IF  (PY( J J.LE.YMIN)  GC  TG  195 
IF(RX(J).GT.XMAX)  GC  TO  195 
IF(RX(J J.LE.XMIN)  GO  TO  195 
IF  (RZ(J) .GE.ZMAX  ) GC  TO  195 
IF  (RZ(  J ) .LE.ZMIN)  GC  TO  195 

TEST  TO  DETERMINE  IF  ATOM  IS  INSIDE  A CIRCLE 
CF  RADIUS  PCE  FROM  THE  EULLET 

CRX  = RX ( J ) — R X ( I ) 

CPY=RY( J)-RY( I ) 

CRZ=RZ( J l-RZ( I ) 

DI$T=DRX~*2+DRY**2+CRZ**2 
IF(DIST.GE.R0E2)  GC  TO  195 
C IST  = SGRT( C 1ST ) 

IF(0IST  .GT.flCEM)  GC  tc  lb5 

ACTUAL  FORCE  FUNCTION  DETERMINED  FR CM  THE 
BCRM-MEYER  BULLET-ATOM  POTENT  IAL 

FCRCE=EXP(PFXA+PEXB*DIST) 

GC  TC  180 
165  CFF=ROE-DIST 

IF(CFF.LT.l.OE-lO)  GO  TC  195  f 

AN  AVERAGE  FORCE  FUNCTION  FOR  THE  BLLLET-ATCM 
INTERACTION  AT  TURNCN 

FCPC£=(EXP(FAC+PEXB*OIST»-PFPTC)/CFF 
IF(FM.GT.FCRCE)  GC  TO  195 
18C  FCC  =FCRC E/CIST 
FA  = FCD*DRX 
F X ( J ) = FX ( J ) +FA 
FM  I)  = F X ( I ) -p  A 
FA  = FQD*DRY 
F V ( J ) = FY ( J ) + F A 
FY(  1)  = FY ( I)-FA 
FA  = FCD*QRZ 
F Z ( J ) = FZ  ( J ) +F  A 
F Z < I ) = F Z ( I ) - FA 
195  CONTINUE 
2CC  CONTINUE 

CC  200  1=2, LD 

I F ( LCUT { I ) .GT.O)  GC  TO  200 
I F = 1 + 1 

XMN  = RX  ( I )-ROE 
X N A X = RX ( I ) +ROE 
YMN  = RY  (I  )-ROE 
Y N A X = R Y ( I ) + ROE 
ZMN  = RZ(I)-ROE 
ZN  A X = RZ(I)+ROE 


TEST  TO  DETERMINE  WHETHER  ATOMS  ARE  CLOSE 
ENOUGH  TO  INTERACT  AND  WHICH  FORCE  FUNCTION 
TO  USE. 

DC  295  J= I P ,LL 
I F (LCUT  (J  ) .GT.O)  C-C  TO  295 
IF (RY ( J ) .GE  .YMAX ) GC  TO  295 
IF  (RY( J ) .LE  .YMIN)  GC  TO  295 
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GC  TQ  295 
GO  TO  295 
GO  TO  295 


IF(RXU).GE.XMAX) 

I F ( RX ( J ) .IE  .XMIN) 

IF  (RZ( J ) .GE.ZMAX > 

IF(RZlJ).LE.ZMIN)  GO  TQ  295 
CRX  = RX ( J )-R  X ( I ) 

CPY=RY ( J ) -R  Y ( I ) 

CRZ=RZ( J )-RZ( I ) 

D IST=DRX**2+0RY**2  + DRZ**2 
I F ( C I ST  .GE.  R0E2 ) GO  TO  295 
DIST=SQRT(DIST) 

IFICIST.GT.ROEM)  GC  TO  265 

ACTUAL  FORCE  FUNCTION  DETERMINED  FR CM  THE 
CU-CU  (GIB-II)  POTENTIAL  FUNCTION 

FCRCE=EXP( FXA+EX9*QIST ) 

GC  TO  280 
265  C FF=RQE-D 1ST 

I F ( CFF . LT  .1  .0  E-LC)  GO  TO  295 

AN  AVERAGE  FORCE  FUNCTION  FOR  THE  ATOM-ATOM 
INTERACTION  AT  TURN  ON 

FCRCE=(EXP(AC+EXB*CIST  )-FPTC)/DFF 
IF(FM.GT.FCRCE)  GC  TO  295 
28C  FCC  = FGRCE/D  1ST 
FA  =FOD=!'CRX 
FX(J)=FX(J)+FA 
FX(I)=FX( I )-FA 
FA=FCO*CRY 
FY ( U ) =F  Y ( J ) +FA 
FY(I)=FV( I ) -FA 
F A=FOD*ORZ 
FZ(J)=FZ<J)+FA 
FZ(  I )=FZ(I )-FA 
295  CCNTINLE 
3CC  CONTINUE 
RETURN 
END 


SLEROUTINE  KPRINT 


THIS  SUBROUTINE  PRINTS  THE  PAGE  FEACING 


CCMMCN/COMl/RXI  600 ) , RY ( 600 ) , R Z ( 6CC ) t LOUT ( 60 C ) , LL  , LD 
CCMMCN/CGM3/EXA, EXB,FXA,PEXA,PEXfi,PFXA,TPERM 
CCMMCN/C0M4/IX, I Y , IZ, I X F , I Y P , SC  X , SCY  , SC  Z , I C EE  P , 

X REX  R 6 2 

CCMMCN/CCM7/R1 ,LSS , SPX,SPZ, CCX , CCY,COZ,  YLAX (20) , I LL 
CCMMCN/COMl  1/RXI ( 6C0 ) , R Y II 6 00 ) , R Z I (600) 
CCMM0N/C0M16/BMAS  ,TMAS 1 ,TMAS2, ABUN1,ABUN2,  I PNUM 
CCMMCN/C0M16A/BX, EY  ,BZ, EDISF,CVR 

COMMON/ COM  17/TIME ,NT,DT,DTI,TEMF,CELS,EV,NPAGE,NLINE 
COMMCN/COM  18/  I TEMP , T AMP , CO RF 

CCMMCN/C0M2C/IH1 ( 20 ) » I HS( 10 ) , I HT ( 6 ) , IHB ( 6 ) , IHTARG ( 2 ) , 

X I F8UL L ( 2 ) 

COMMON/ COM2 1/ IHBP CT ( 2 ) , IHT POT ( 2 ) , SUP ( 2 ) , Ih2  ( 1 0 ) , I OR  , 

X L SLR 

COMMON/ C0M22/R0EA ,RCEB,ROEC ,R0EC2,CP0,CP  l,  CF2tCP3  ,CFC  . 
X CF1,CF2,CGC1,CG02,CG81,CGB2,CGF1,CGF2,  IPRCC- 
96 IC  FORMAT ( 2 CX , 20A4  , / ,4CX, 10A4, /// ) 

964C  FORMAT ( ' ( • , 13 , • » PLANE  ',2A4,»  SURFACE,  PRIMARY', 

X • ENERGY  = ' ,F7. 1 » ' EV,  CRYSTAL  SIZE  (',12,'  X ',12, 

X • X ',12,'  ),  IMPACT  POINT  ',19/) 

9645  FORMAT ( • TARGET  -',2A4,20X,'  BULLET  -',2A4,7X, 

X 'LATTICE  UNIT  =',F7.4,'  ANG’) 

965C  FORMAT ( ' MASSI  = ',  F7 . 2 , 2X ,' ABUNDANCE  =',F5.2,6X, 
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X ' MASS  = • ,F7. 2 ,9X,  ' LATTICE  TEMPERATURE  = • , A4  , 

X • THERMAL  AMPLITUDE  =',F7.4) 

9652  FORMAT ( ' MASS2  = • , F7 .2 , 2X , ' A6UNC ANCS  =',F5.2, 

X 29X, 'CORRELATION  =',F5.3,//) 

96  5 5 FCFMAT ( 12X,2(2A4,1X)  , 2A4, 3X , 5hP EX A = , F9 . 5 , 2X , 5HPEX B = , 

X F9.5,2X,5HPFXA=,F9.5/) 

966C  FCFMAT ( 12H  POTENTIAL  , 2 ( 2A 4 , IX  ) , 2A4 , 3X , 5H E X A =,F9.5, 

X 2X.5HEXB  =,F9.5 ,2X,5HFXA  =tF9.5) 

9665  FORMAT ( 20H  TARGET  POINT  ON  CRYSTAL  X =,F7.3,*,  Y = ', 
X F7 . 3 , ' , Z = ' , F7  .3 , 13H , COS  TO  X =,F8.5, 

X 12H,  CCS  TO  Z = » F8 . 5 ) 

967C  FCFMAT ( 20H  PRIMARY  START  POINT  ( LU  ) X =,F7.3,5H,  Y =, 
X F7 .3  » 5H  , Z =,F7.3/) 

9675  FCFMAT ( / ' TIMESTEF  ',I4,40X,'  ELAPSED  TIME  (SEC)  = ', 

X E 10 .4  * ' , LAST  TIMESTEP  WAS  = • , E 10 . 4/ ) 

9725  FCFMAT ( • CUT-OFF  AT',F5.2,',  WHEN  R =',F6.3,'  LU, MORSE 
X POTENTIAL  PARAMETERS  ARE  ' ,10  A4  , // , 10X  , • CC-D1  =', 

X F8.4, ' , CGD2  =',F8.4,',  CGB1  =',F8.4,',  CGB2  =', 

X F8.4,  ' , CGF 1 = ' ,F8.4,  ',  CGF2  =',F8.4,/i 
9 7 3 C FCFMAT ( ' WHEN',F8.4,'  <R<»,F8.4,»  THE  MATCHING  POTENTI 
XAL  PARAMETER  ARE',//,3X,'  CPC  =',E12.4,',  CPI  =', 

X E 12 . 4 , '8  CP2  = ' , E 12 . 4 , ' , CP3  = ' ,E12 .4  , /7X  , ' CFO  = ', 

X E 12 .4  , ' , CF1  = ' » E 12. 4 » ' » CF  2 *',E12.4,//) 

974C  FORMAT ( 11 h THE  FIRST  ,I3,22H  ATOMS  ARE  BLCCKED,  ,12, 
X 24H  LAYERS  ARE  FREE  TC  MOVE,/) 

WRITE  ( 6,9610)  IFl.IHS 

WRITE  ( 6,9640)  I CR , SUR , E V , IX  , I Y , I Z , I PNU M 

WRITE  ( 6,9645)  I H T ARG , IH BU LL , C VR 

WRITE  (6,9650)  TM A S 1 , ABUN 1 , BMAS  , I TEMP , T AMP 

WRITE  (6,9652)  TM AS2 , ABLN2 , COR F 

WRITE  ( 6,9660)  I HT  , E X A , EXB  , FXA 

WRITE  ( 6,9655)  I H B , PEX A , PEX B , P FXA 

WRITE  ( 6,9665)  B X , BY , B 2 , CO  X , CCY  , CCZ 

WRITE  ( 6,9670)  RX I < 1 ) , R YI ( 1 ) , R 2 I ( 1 ) 

IFUPRCG.NE.4)  GO  TO  130 

WRITE  ( 6,9730)  R CE A , ROEB , C PO , C P 1 , CP2 , C P 3 , C FC  , CF1  , CF2 
WRITE (6,9725)RQEC,RCES»IH2,CGD1,CGC2,CG31,CGB2,CGF1, 

X CGF  2 

IOC  WRITE  ( 6,9675)  NT, TIME, DT 
N L I NE= 1 
RETURN 
ENC 


SUBROUTINE  KPNT 


THIS  SUBROUTINE  PRINTS  THE  HEADING 
FOR  OUTPUT  THAT  IS  FILEC 


CCMMCN/C0M4/IX, I Y,I2, I X P , I YP ,SC X , SCY , SC Z , I C E E P , 

X F EX  , RB  Z 

CCMMCN/COM  16/BMAS  ,TMAS1  ,TMAS2  , ABLM  ,ABUN2 , 1 FNUM 
CCMMCN/C0M16A/BX, B Y , BZ , ED  I $ P , C VR 

COMMON/ COM  17/TIME ,NT,DT ,DTI , TEMP, CELS, EV,NPfiGE,NLINE 
CCMMCN/COM20/IHK20),  IHS(IO)  , I HT  ( 6 ) , IHB  ( 6 ) , IHTARG  ( 2 ) , 
X IHBUL  L ( 2 ) 

CC MMCN/C0M2 1/ IHBPCT (2) , IHTP0T(2 ) , SUR ( 2 ),  I H2  ( 1 0 ),  I CR , 

X LSLR 

96  4 C FORMAT ( ' (',13,')  P LANE  • , 10X , ' I ON  ENERGY  =',F7.1,10X, 
X 'CRYSTAL  SIZE  C,I2,'  X',I2,'  X ',  12 ,')',// ) 

96 A 5 FORMAT ( ' TARGET  -',2A4,20X,'  BULLET  -*,2A4,6X, 

X 'POTENTIAL  ' , 2 ( 2A4  , IX ) , 2 A4 ) 

96  5 C FORMAT ( ' MASS  1 = ' , F 7 . 2 , 2X  , • ABU N C ANCE  =',F5.3,4X, 

X 'MASS  =' ,F7.2,16X,2(2A4,1X) ,2A4) 

9652  FORMAT  ('  MASS  2 =', F7 .2 , 2 X ,' ABUNDANCE  =',F5.3,//) 
WRITE  (6,9640)  IOR , E V , I X , I Y , I Z 
WRITE  (6,9645)  I H T ARG , I HBUL L , I HT 
WRITE  (6,9650)  TM A S 1 , ABUN l , EM AS  , I HB 
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WRITE  (6,9652)  TMAS2,A6UN2 

RETURN 

END 


SL6RGUT INE  HIST0( IN , FREC1 , R ANGE  , T ITL E 1 ) 


THIS  SUBROUTINE  CRAWS  THE  HISTORAM  (FREQUENCY 
VS.  EJECTEC  ANGLE)  FOR  THE  TYPE  1 HEAVY  ATOMS 


DIMENSION  FREQL (20  ) , RANGE ( 20 ) , JOLT ( 20 ) , T ITLcl (24 ) 
CIMENSICN  POUT ( 2 0 > 

CATA  NOTH/'  '/,«/•****»/ 

2 FORMAT ( *0  INTERVAL ',2X,19(F5. 1, IX) ,F5.1) 

4 FCFNAT(//,8X,24A4/) 

£ FORMATt  'OFREQUENC Y • ,2016) 

6 FORMAT ( 'OINTERVAL ' , 2X , 19( F5.1 » IX ) , F5. 1) 

7 FCFMAT ( • -•  , 12  8 ( ' -•  ) ) 

WRITE  (6,4)  TITLE  1 
JN=  IN 

INT  = 0 

IF(JN.GT.O)  GO  TO  10 

I NT  = 1 

JN=-JN 

10  CC  12  1=1, JN 
12  JCLT(I) = FREC1 ( I ) 

WRITE  (6,5)  ( JOUT(I ) ,1  = 1, JN) 

WRITE  (6,7) 

FINC  LARGEST  FREQUENCY 

FMAX=0. 

CC  20  1=1, JN 

I F (FREQ1 ( I ) .GT.FMAX)  FM AX =F PEQ 1 ( I ) 

2C  CONTINUE 

SCALE 

JSCAL=1 

I F ( FMAX  .GT  .60. ) J SC AL= ( FMAX +59 . )/60. 

F SC AL  = J SCA  L 
CC  50  1=1, JN 
SC  JCLT ( I ) =NOTH 
MAX=FMAX/FSCAL 
CC  8 J 1=1, MAX 
X =M AX- ( 1-1 ) 

CC  70  J = 1 , J N 

IF(FREQ1<  J)/FSCAL.C-E.X)  JOUT(  J)=K 
7 C CONTINUE 
IX=X*FSCAL 

£0  WRITE  (6,2)  I X , ( JCUT ( J ) , J = 1 , JN) 

2 FCFMAT ( 16, 4X,20(2X,A4) ) 

WRITE  (6,7) 

IF(INT.EQ.l)  GO  TC  16 
WRITE  (6,3)  (RANGE (J  ), J = 1,JN) 

GC  TO  15 
16  CC  51  1=1, JN 
51  FCLTd  ) =RANGE  ( I ) 

WRITE  (6,6)  (POUT (I ) ,1=1, JN) 

15  CONTINUE 

WRITE  (6,100) 

ICC  FCPMATt  IX, 'MIDPOINT'/) 

RETURN 


SLEPOUTINE  I-  IS  TO  ( IN , FREC2  , R ANGE  , T ITLE2  ) 


C 
C 

C THIS  SUBROUTINE  CRAWS  THE  H ISTORAM  (FREQUENCY 

C VS.  EJECTEC  ANGLE)  FOR  THE  TYPE2  LIGHT  ATOMS 

C 
C 

CIMENSICN  FFEQ2 ( 2 J ) , R ANGE < 2 0 ) , JOUT ( 2 J ) , T ITL E2 ( 24 ) 

C IMENS I ON  POUT ( 2Q  ) 

CATA  NOTH/'  •/,«/•****•/ 

2 FCFMAT(  I6,4X,2J(2X,A<t)  ) 

3 FORMAT ( *0  INTERVAL' , 2X , 19 ( F5 . 1 , 1 X ) , F5 . 1 ) 

4 FCFMAT(//,8X,24A4/) 

5 FORMAT ( 'OFREQUENCY' ,2016) 

6 FORMAT ( '0  INTERVAL  *,2X,19(F5.1,1X),F5.1) 

7 FC F MAT ( • » 128  ( •-•  )) 

WRITE  (6,4)  TITLE  2 
J N = I N 
I NT  = J 

IF(JN.GT.O)  GO  TO  10 
I N T=  1 
JM-JN 

1C  CC  12  1=1, JN 
12  JCLT(I)=FREC2(I) 

WRITE  (6,5)  ( JQUT(I ) ,1=1, JN) 

WRITE  (6,7) 

C 

C FINC  LARGEST  FREQUENCY 

C 

F M A X=  0 . 

CC  20  1=1, JN 

IF (FREQ 2(1  ) .GT.FMAX)  FMAX  =F REQ2 ( I ) 

2C  CONTINUE 
C 

C SCALE 

C 

JSC AL=1 

I F(FMAX  .GT.60. ) J SCAL= ( FM AX +59  . ) /6 0 . 

FSC AL  = J SCA  L 
CC  50  1=1, JN 
5 C JCLT ( I ) =NOTH 
MAX=FMAX/FSCAL 
CC  80  1=1, MAX 

X = MAX-(  1-1  ) " 

CC  70  J = 1 , J N 

IF(FREQ2 (J  )/FSCAL  .GE.X  ) JOUT( J)=K 
7 C CONTINUE 
IX=X*FSCAL 

80  WRITE  (6,2)  IX, ( JCUT ( J) , J=1 , JN) 

WRITE  (6,7) 

I F ( INT. EQ.l  ) GO  TO  16 

WRITE  (6,3)  ( RANGE ( J ) , J = 1 , J N ) 

GC  TO  15 
16  CC  51  1=1, JN 
51  POUT  II ) =R ANGE l I ) 

WRITE  (6,6)  (POUT  (I),I=1»JN) 

15  CONTINUE 

WRITE  (6, ICC) 

ICC  FCPMAT(  IX, • MIDPOINT*/ ) 

RETURN 

ENC 


SLEROUTINE  KFOLDA 
C 
C 

C THIS  SUBROUTINE  STATISTICALLY  FOLDS  THE 

C TARGET  IMPACT  AREA  FOR  THE  FCC  1100)  CRYSTAL 

C LATTICE  TO  REPRESENT  FULL  COVERAGE 

C 
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CCMMON/CCM1/RX(600) ,RY(63J)  , RZ ( 6 0 0 ) , LC J T ( 6 0 : ) ,LL,LC 
CCMM0N/C0M4/IX, I Y,IZ, IXP,IYP,SCX,SCY,SCZ,ICEEP, 

X PEX.ReZ 

CC  N M0N/C0M4  A/  LNUM6J0),LNUMI(600)  ,LNUMF(6JC) 
CCMM0N/C0M21/IHBPCT(2)  , I HT P CT  ( 2 ) , SUR  ( 2 ) , I F2  ( 10  ) , I CR  , 
X LSLR 
I LA YER=  1 
N F IRST  = 2 

DETERMINE  THE  FIRST  ATOM  AND  THE  LARGEST 

NUMBERED  A TGM  IN  EACH  LAYER 

10  GC  TO  ( 20,30,40, 5C,6C,70) , ILAYER 

20  N LAY  =( < LL- 1 ) / ( IY/2)  )/2+2 
K =NFIRST 
NLAY 1 = N LAY 
GC  TO  8 C 

30  N L AY= ( ( LL-1 )/( IY/2 ) )+l 
K=NFIRST  +IX/2 
GC  TO  80 

40  N LA Y = NLAY 1 + NLAY  -I 
K = NF IRST 
GC  TO  80 

SC  NLAY  = NLAY1  + NLAY  -2 
K=NFIRST  +IX/2 
GC  TO  80 

6C  NLAY  = NLAY1  + NLAY  -3 
K = N F I RST 
GC  TO  80 

7 C NLAY  = NLAY1  + NLAY  -4 
K=NFIRST  +IX/2 
GC  TO  8 C 

FIRST  (DIAGONAL)  FOLC  - ALL  LAYERS 

80  DC  90  J = NF I RST , NL  A Y 

LNUM  J)  =LNLY(  J)  + LNUM(K) 

IF(J.EQ.K)  LNUM( J) =2*LNUM( J ) 

LNLM( K) =LNUV( J ) 

K =K  + 1 X 

IF (K.GT .NLAY)  K=K-NLAY 1+2 

9 C CONTINUE 

DC  100  J = NF  IRST, NLAY 
LNUM ( J ) =LNUM( J ) /2 
1JC  CONTINUE 

SECOND  ( HCR I ZONTAL  ) FOLD  - ODD  LAYERS 

105  ICCL'NT=1 
I S ET=G 

GC  TO  ( 110, 160,  11C, 160, 110, 160)  , ILAYER 
11C  K = N LA Y - I X / 2 

DC  150  J=NF  IRST, NLAY 
LNLM ( J) =LNUM( J) + L NON ( K ) 

IF(J.EQ.K)  LNUM (J)=2*LNUM(J) 

LNLM(K)=LNUM(J) 

ICCUNT= ICOLNT+1 
IF ( ISET .EQ.  1)  GO  TC  130 
IF ( ICCUNT.GT.  iX/2  +1)  GO  TO  120 
K=K  + 1 
C-C  TO  150 
1 2 C I S ET  = 1 
K =K - 1 X + 1 
I C CUNT  = 1 
GC  TO  150 

13C  IF  (ICCLNT.GT.  IX/2)  GO  TO  140 
K =X  + 1 
GC  TO  150 
1 4 C I S E7=0 
K=K-IX+1 
I CCLNT  = 1 
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1 5 C 
16C 


170 

180 

150 

2 0 C 
2 1C 

22C 

22C 


24C 

250 

26C 

27C 

2 ac 


250 


CONTINUE 
GC  TC  210 

SECOND  (HORIZONTAL)  FOLC 

K=NLAY-IX/2+l 
DO  200  J = NF  IRST *NLAY 
LNUN ( J ) =LNUM ( J ) + LNUM (K  ) 
IF(J.EQ.K)  LNUM( J)=2*LNUM(J) 
LNLM(K)=LNUN( J) 

I COUNT  = ICOUNT+1 

IF(ISET.EQ.l)  GO  TC  180 

IF ( ICOUNT.GT.IX/2  ) GO  TC  170 

K = K + 1 

GC  TO  2C0 

I S ET= 1 

K=K-IX+1 

I COUNT  = 1 

GC  TC  2C0 

IF  ( ICCUNT.GT.IX/2+1)  GC  TO  190 

K=K  + 1 

GC  TO  2C0 

I 'ET  = 0 

K = K— I X + 1 

I COUNT  = 1 

CONTINUE 

OC  220  J = NF  IRST, NLAY 
LNLM  J)=LNUM(  J) /2 
CONTINUE 

THIRD  (VERTICAL)  FCLD  - 

I COUNT  = 1 
I S ET=0 

GC  TO  ( 23C,  280, 23C, 280,230, 280)  , 
K =NF IR  S T+ 1 X/2 
CC  270  J=NF  IRST, NLAY 
LNLM(J)=LNUM(J)+LNUM(K) 
IF(J.EQ.K)  LNUM ( J ) =2*LNUM ( J ) 
LNLM(K) =LNUM( J) 

I CCUNT=  ICCUNT+1 
IF(ISET.EQ.l)  GO  TC  250 
IF(ICCUNT.GT.IX/2  +1)  GO  TO  240 
K =K—  1 
GC  TO  270 
I SET  = 1 
K = K + IX-  1 
I C CLNT= 1 
GC  TO  270 

IF  ( ICCUNT.GT.IX/2)  GO  TO  260 

K =X- 1 

GC  TO  270 

I S ET  =0 

K = K + IX-  1 

I COUNT  = 1 

CONTINUE 

GC  TO  330 

THIRD  (VERTICAL)  FOLD  - 

K=NFIRST+IX/2-l 
CC  320  J=  NF  IRST , N L AY 
L N(J N ( J ) =LNUM(  J ) + LNUM  (K  ) 
IF(J.EQ.K)  LNUM ( J ) =2* LNUM ( J ) 

LNLM ( K ) =LNU M ( J ) 

I COUNT  = ICCUNT  + 1 

IF(ISET.EQ.l)  GO  TC  300 

I F ( ICOUNT.GT.IX/2)  GO  TC  290 

X = K - 1 

GC  TO  320 

I 5ET=1 

K =K  + I X-  1 


- EVEN  LAYERS 


CCC  LAYERS 
I LAYER 


EVEN  LAYERS 
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I CCUNT  = 1 
GC  TO  32  0 

30C  IF(  ICGUNT.GT.IX/2  + 1)  GO  TO  310 
K = K — 1 
GC  TO  320 
3 1C  I SET  = 0 
K=K+IX-1 
ICCUNT=  1 
32C  CONTINUE 

33C  CC  340  J=NF  IRST  , NLAY 
LMH(  J)  =LNUM(  J)/2 
34C  CONTINUE 

I LA Y ER  = ILAYER+1 
NFIRST=NLAY+1 

IF ( ILAY  ER.GT.IY)  GC  TO  350 
GC  TC  1C 
35C  RETURN 
ENC 


SIEROOTINE  KF0LD3 


THIS  SUBROUTINE  STATISTICALLY  FOLCS  THE 
TARGET  IMPACT  AREA  FOR  THE  FCC  (UC)  CRYSTAL 
LATTICE  TC  REPRESENT  FULL  COVERAGE 


CCMMCN/CCM1/RX(60C)  ,RY(60J),RZ(60Q),LCUT(60C),LL,LC 
C CMMCN/ COM4 /IX, IY,IZ,IXP,IYP,$CX,SCY,SCZ,ICcEP, 

X REXtRBZ 

CCMMC N/C0M4A/  LNUM(60J),LNUMI(6:0)  , LNUM F (600 
CCMMCN/C0M21/  IH6PCT (2 ) , IHTPGT ( 2 ) , SUR(2)  , I F2  ( 1 0 ) , I CR  , 
X L S UR 
I LA YER=  1 
NF  IRST  = 2 

DETERMINE  THE  FIRST  ATOM  AND  THE  LARGEST 
NUMBEREO  ATOM  IN  EACH  LAYER 

1C  GC  TO  (20,3C,40,3C,40,3C) ,ILAYER 
2C  NLAY=( ( IX+1  )/2)*(  ( IZ  + l)/2)  + l 
N L A Y 1 =N  L AY 
GC  TO  50 

3C  NLAY  = NLAY  + ( IX/2  )* ( IZ/2 > 

GC  TC  50 

4 C NLAY  = NLAY  + NLAY1-1 
5 C I C CUNT=  1 

GC  TO  ( 6C, 90, 60,  90,60,90 , 1 LAYER 

SECOND  (HCRIZONTAL)  FCLC  - ODD  LAYERS 

60  K =N LAY- I X/ 2 

CC  80  J=NFIRST,NLAY 
LNLM( J) =LNUM( J j +LNUM(K ) 

IF(J.EG.K)  LNUM ( J ) = 2*LNUM(J) 

LNLM( K) = LNU  M J ) 

I COUNT  = ICCUNT  + 1 

IF  (ICCUNT.GT.  IX/2  +1)  GO  TO  70 
K =K+  1 
GC  TO  3 C 
70  K = K- 1 X 
I CCUNT=  1 
8C  CCNTINUE 
GC  TO  120 

SECOND  (HCRIZONTAL)  FOLC  - EVEN  LAYERS 

SC  K = N LAY-IX/2  +1 

CC  110  J = NF  IRST , NL AY 
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LNUM  J)=LNUM(  J)+LNUM(K  ) 

IF(J.EQ.K)  LNUM ( J ) = 2*LNUM(J) 

L NLM  ( K ) =LNU M ( J ) 

I CCUNT  = ICCUM+1 
I F ( ICOUNT  .GT.IX/2  ) GO  TG  LOG 
K =K  + 1 
GC  TO  110 
ICC  K = K- 1 X+  2 
I CCUNT  - 1 
11C  CONTINUE 

12C  DC  130  J=NF IRST .NLAY 
LNU  N ( J ) = LNUM l J 1/2 
120  CONTINUE 
ICCUNT=  1 

GC  TO  < 140,  150,  140, 19C,  140,  ISO)  , ILAYER 

THIRD  (VERTICAL)  FOLD  - CCD  LAYERS 

14G  K=NFIRST+IX/2 

CC  160  J = NF  IRST, NLAY 
LNUM  J)=LNUN(  J)  >LNUM(K) 

IF(J.EQ.K)  LNUM ( J ) = 2*LNUM(J) 

LNLM( K) =LNUM( J ) 

ICCUNT= ICCUNT+1 

IF  (ICCLNT.C-T.  IX/2+1)  GO  TC  150 

SPil  160 

15C  K =K  + I X 
I C C UNT  = 1 
1 6 C CONTINUE 
GC  TO  220 

THIRD  (VERTICAL)  FOLD  - EVEN  LAYERS 

1 5 G K=NFIRST  +IX/2  -1 

CC  210  J = NF  IRST, NLAY 
LNUM  J)=LNUM(  J) +LNUM(K  ) 

IF(J.EG.K)  LNUM ( J ) = 2*LNUM(J) 

LNLM( K) =LNUM( J ) 

ICCUNT=ICCUNT+1 
IF  ( ICCUNT.GT  . I X / 2 ) GO  TO  200 
K =K- 1 
GC  TO  210 
2 CC  K=K  + IX-2 
I CCUNT  = 1 
210  CONTINUE 

22C  CC  230  J=NF IRST, NLAY 
LNUM  J)=LNUM(J)/2 
2 2 C CONTINUE 

ILA YER= ILAYER+1 
N F IRST=NLAY  + 1 

IF( ILAYER.GT.IY)  GO  TO  240 
GC  TO  10 
24C  RETURN 
ENC 


SLERCUTINE  KFCLDC 


THIS  SUBROUTINE  STATISTICALLY  FOLDS  THE 
TARGET  IMPACT  AREA  FOR  THE  FCC  (111)  CRYSTAL 
LATTICE  TC  REPRESENT  FULL  COVERAGE 
***  VALID  ONLY  FOR  THE  FIRST  LAYER  AND  IX  £,  I*  =9 


CCMMCN/CCM1/RX(600 ) , RY ( 600 ) , RZ ( fcCC ) , LCUT ( 6C C ) , LL , LD 
CUNCN/C0M4/IX, I Y , I Z , IX F,  I Y P , SC X , SC Y , SC Z , ICEEP, 

X R EX , RE  Z 

CCN  NCN/C0M4 A/  LNUM (600 ), LNUM I (600 ) ,LNUMF ( 6 GC) 


' 


! 

5 


; 


T 


A 


'4 


i 
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CCPPCN/C0M21/IH3PCT(2)  , IHTPCT  ( 2 ) ,SUR(2)  ♦ IH2  ( 10) , ICR , 

X LSLR 
CC  1 J = 2 , L L 
LMJMF  (J  ) = 0 
LNUPH  J )=LNUM(J  ) 

1 CCNTINUE 
I CLT=0 

I T I PE  = 1 

NLAY=(  ( LL-2  )/(  IY/2)  )/2  + 2 

2 NF IPST  = 2+  I X/2 
I C CUNT=  C 

I SET=0 
K=NLAY-IX/2 
J =NF I RS  T 
5 CCNTINUE 

FIRST  ( + SLOPE  DIAGONAL)  FOLD  - FIRST  LAYER 

LNLP(  J)  =LNUM  J)  + LNUP  (K  ) 

IF  (J.EC.K)  LNUM< J)=2*LNUM( J) 

LNLM( K ) =LNUP( J ) 

ICCLNT=ICCLNT+1 

IF  (ISET.EQ.O  .AND.  ICOUNT.EQ.l)  GO  TO  10 

IF  (ISET.EQ.O  .ANC.  IC0UNT.EQ.2)  GO  TO  15 

IF  (ISET.EQ.O  .AND.  I COUN T . E Q . 3 ) GO  TO  20 

IF  (ISET.EQ.2  .ANC.  ICOUNT.EQ.l)  GO  TO  22 

IP  (ISET.EQ.2  .ANC.  ICOUNT  .EG  .2  ) GO  TO  15 

IF  (ISET.EQ.2  .AND.  ICOUNT.EQ.3)  GO  TO  10 

IF  (ISET.EQ.2  .ANC.  ICGUNT.EG.4)  GO  TO  35 

I F ( ICGUNT.EG.A-  .OR.  I COUNT  . E G . 8 .OR.  ICCUNT  .EG. 13 
X .CR.  ICOLNT.EQ  .17 ) GO  TO  25 
IF  (ICCUNT. EQ.  20)  GO  TC  30 
J=J  + 1 
K = K +5 
GC  TC  5 

1C  J=J  + < IX-3  )/2 
K = K— I X 
GC  TO  5 
15  J=J  + 1 

K = K + ( IX+1 )/2 
C-C  TO  5 

2G  ISET=ISET+1 
I C CLNT  = 0 
22  J = J+2 
K=K-14 
GC  TO  5 
25  K =K  — 19 
J=J  + 1 
GC  TO  5 

3C  ISET=ISET+1 
I CCLNT=  0 
J =0  + 1 
K =K+5 
GC  TO  5 

35  IF  (IOUT.EQ.1)  GO  TC  150 
GC  TO  (50,50,100) ,ITIPE 

SECOND  (-  SLOPE  DIAGONAL)  - FIRST  LAYER 

5 C CCNTINUE 
NFIRST=2 
I F L AG=  0 
ICCLNT=C 
I S ET  = 0 
I CL  T = 0 
K =N  LA Y 
J =NFIRST 
55  CCNTINUE 

LNLM( J)=LNUP( J)+LNUM(K) 

IF  (J.EC.K)  LNUM( J)=2*LNUM( J) 

LNLM(K)  = LNUM  J) 

ICCUNT= ICOUNT+1 


nnn 


L 
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60 


65 


7 C 
72 


75 


77 

ec 

85 

95 


IF 

( I SET  .EQ  .0 

.ANC. 

ICGUNT.EC.l ) 

GO 

TC 

6 C 

IF 

( ISET.EC.O 

• AND. 

ICQUNT. EQ. 2) 

GO 

TO 

65 

IF 

(ISET.EC.O 

.AND. 

ICQUNT. EG. 3) 

GO 

TC 

7 C 

I F 

( I SET  .EQ  .2 

.AND. 

I COUNT . EQ . 1 ) 

GO 

TC 

72 

IF 

( ISET.EQ  .2 

.AND. 

ICQUNT. EQ. 2) 

GO 

TO 

65 

IF 

(ISET.EG.2 

.AND. 

ICOUNT.EC  .3  ) 

GO 

TO 

6 C 

IF 

(ISET.EQ. 2 

.AND. 

ICQUNT. EQ. 4) 

GO 

TQ 

85 

ir  .«ntu.  ii^uuiN  ) icy  «*r  i ou  iu  03 

IF  (ICCUNT.EQ.4  .OR.  ICQUNT . E Q . 1 7 ) GO  TC  75 

IF  (ICCUNT.EQ.8  .OR.  ICGUM.EG.I3)  GO  TC  77 
IF  ( ICQUNT . EQ .20 ) GO  TO  80 


J = J + 1 
K=K-4 
GC  TC  55 
J = J + 5 
K =K-5 
GC  TO  55 
J = J+i 
K = K-4 
GC  TO  55 
I $ET=ISET+I 
I C CUNT=  0 
J = 0 + 3 
K =K  + 3 
GC  TO  55 
0 = J+2 
K = K +7 
C-C  TO  5 5 
J = J+1 
K = K + 1 1 
GC  TO  55 
I SET  = I S ET+ 1 
I CCUNT  = 0 
J=J  + 1 
K = K-4 
GC  TO  55 

GC  TO  (103,55,150) ,ITIME 

ICLT=1 

GC  TO  2 


THIRD  (VERTICAL)  RQLD  - FIRST  LAVER 


ICC 


1 1C 

12c 

1 20 

140 

15C 


CONTINUE 

I CCUNT  = 1 

I S ET  = 0 

NFIRST=2 

K=NFIRST+IX/2 

DC  140  J=  N F IRST , N L A Y 

LNUN(  J)  =LNUN( J)+LNUM(K) 

IF  (J.EQ.K)  LNUM< J)=2*LNUM( J ) 
LNUN(  K)  =LNUM  J ) 

I CCUNT = ICCUNT  + l 

IF  (ISET.EC.l)  GO  TC  120 

IF  (ICCUNT.C-T.  IX/2  +1)  GO  TO  110 

K =K  - 1 

GC  TC  140 

I SET=1 

K=K+IX-1 

I CCUNT  = 1 

GC  TO  140 

IF  (ICCLNT.GT.  IX/2)  GO  TO  130 

K =K-1 

GC  TC  140 

I S ET  =0 

K=K+IX-1 

I CCUNT  = 1 

CONTINUE 

GC  TO  ( 150, 15C, 50  ) , I T I ME 
I T IME= I T I ME  +1 
I C UT  =0 

I CCUNT  =0 

DC  152  J = 2 , NLAY 
LNUM ( J ) =L  NU  M J ) / 8 
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L 


L NLMF ( J)=LNLMF(J) + LNUM (J) 
152  CCMINUE 

IF  (ITIHE-GT-3)  GC  TO  160 
CC  155  J=2,NLAY 
LNLM  ( J ) =LNUN I ( J ) 

155  CCNTINUE 
GC  TO  2 

16C  CC  170  J = 2 * NL  AY 

LNUM<  J)=LNUMF(  J)/4 
1 7 C CCNTINUE 
K =N  LAY+  1 
CC  180  J=K,LL 
LMf'  ( J)  =6*LNUM<  J ) 

18C  CCMINUE 
RETURN 
END 


SUBROUTINE  KPRQB 


K PRC  E IS  DESIGNED  TO  STUDY  THE  PRCBAEILITY 
THAT  A RING  OF  ATOMS  (AT  A GIVEN  CISTANCE  FRCM 
THE  LAYER  CENTER)  WILL  SPLTTER 


CCMMCN/CGM1/RX(60C) ,RY(600 >,RZ (6CC),LCJT(60C),LL,LD 
CCMMCN/COM4/IX, IY  , IZ,  IXF,IYP,SCX,SCY,SCZ,IDEEP, 

X R EX  » RE  Z 

CCMMCN/C0M4A/  LNUM  (60G ) , LNUM I ( 6 CC ) ,LNUMF(6CC) 
CCMMCN/CCM1  1/RXI  (60C),RYl(oC0) ,RZI  (600) 

C CMMCN/ COM  2 1/ IhBPCT ( 2 J , IHTPCT ( 2 ) , SLR ( 2 ) , IH2 ( 10 ) , I CR , 

X US  UR 

CC MMCN/CCM24/  DNN ( 10G) ,KATON( IOC ) ,KSPUT ( ICO  ) , PROB ( ICO  ) 
CCMMCN/COM24A/  KSLM,KRING,RCI6T(6C0) 

CC  2 1=1,100 
KATCM(I ) =0 
K S PUT  ( I ) =0 
FRCE ( I > = 0.0 
2 CCNTINUE 

DETERMINE  THE  DISTANCE  PRCN  EACH  A T CM  TO  THE 
INITIAL  (0,0)  POINT  OF  TARGET  IMPACT  AREA 

CC  5 1=2,  LL 
XC I ST  =R  X I l I )-RBX*SCX 
2 C I ST  =R  Z I ( I ) — R B Z*  SCZ 
SCIST  = XCIST**2  +ZCIST**2 
R C I ST ( I ) = SD  1ST *=  . 5 
5 CONTINUE 
K F I RST= 1 
I L A YER= 1 
NF IRST=2 

IF  (I0R.EQ.110)  GC  TO  170 
IF  (IOR.Efl.llli  GG  TO  200 

DETERMINE  THE  FIRST  ATOM  AND  THE  LARGEST 
NUMBERED  ATOM  IN  EACH  LAYER  FOR  100  LATTICE 

1C  GC  TO  ( 20  » 3C ,40  » 5C  ,6C  ,70) , ILAYER 
2C  N LAY  =(  C L L-  1 > / ( IY/2) )/2  + 2 
NL  A Y 1 =N  L AY 
C-C  TO  8C 

2C  N LAY= ( ( LL-1 )/( IY/2 ) ) + l 
GC  TC  80 

4C  MAY  =NL A Y 1 + NLA Y -1 
C-C  TO  8 C 

5 C MAY  = NLAY1  + NL  A Y -2 
GC  TO  8 C 

60  MAY  = N LAY  1 + NL  AY  -3 
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GC  TO  ac 

7 2 NLAY  = NLAY1  + NLAY  -4 
C-C  TO  3 C 

SORT  DISTANCES  IN  A LAYER  INTO  ASCENDING  CRDER 

3C  N1=NLAY-1 

DC  100  I = NFIRST,N1 
N F 1= 1 + 1 

CC  9 0 J = NF1 .NLAY 

IF  (RCIST(jLGE.RDIST(D)  GO  TO  9Q 
SAVE=RC I ST  ( I ) 

R0 1ST ( I ) =RC  1ST ( J ) 

RC 1ST ( J )=S AVE 
9 C CONTINUE 
ICC  CONTINUE 

ESTABLISH  a RING  at  EACH  CISTANCE  (+  OR  - DR) 

K=KFIRST 

CC  12  J I = NF  IRST  , NLAY 
J = I + 1 

IF(J.GT.NLAY)  GO  TO  110 

IF  (RDISTI JJ-RDIST(I)  .GT. 0.001)  GC  TO  110 
GC  TO  120 

11C  CNN  <K ) = RG IST(  I ) 

K=K  + 1 

1 2C  CONTINUE 


REINITIALIZE  THE  DISTANCE  FROM  EACH  ATOM 
TO  THE  REFERENCE  (0,0)  POINT  OF  IMPACT  AREA 

CC  145  1=2, LL 
X C I ST  =R  X I ( I )-RBX*  SCX 
ZD  I ST  = R Z I ( I )-REZ*SCZ 
. SCIST  = XCI ST**2  + ZD  I ST  *=*2 
PCIST ( I ) = SC  1ST**. 5 
145  CCNTINUE 


DETERMINE  THE  NUMBER  OF  ATOMS  AND  SFUTS  FPCM 
ATOMS  IN  EACH  RING 

K R ING  = K -1 

CC  150  J=KFIRST ,KPING 
CC  140  I=NF IRST, NLAY 

IF(ABS(RCIST( I)-ONN(J) ) .LT. 0.001 ) GO  TO  13C 
GC  TO  140 

130  K A 7C  M J ) =K  A TOM ( J ) +1 

KSPUT< J)=KSPUT( J)+LNUM( I)' 

14C  CCNTINUE 


C 

c 


DETERMINE  THE  PROBABILITY  AN  ATOM  IN  EACH 
RING  SPUTTERING 

IF (KSUM . EG  - 0 ) KSUM=1 

IF  (IQR.EQ.100)  PROB l J) =FLO AT ( K S PLT ( J ) ) / (FLCAT(KSUM  )* 

x e .o ) 

IF  (IOR.EQ. 110)  PROS  <J)=FLOAT  (KSPL'T  (J)  ) / (FLCAT  (KSUM  ) * 
X 4.0) 

IF  (IOR.EQ. Ill)  PROB(J)=FLOAT(KSPUT(J) ) /(FLCAT(KSbM)* 
X t.C) 

15C  CCNTINUE 

ILAYER= ILAYER+1 
NF IRST=NLAY+1 
K FIRST  = K 

155  IF  ( ILAYER.GT.IY)  GO  TO  160 
IF  (IQR.EQ.100)  GC  TO  10 

IF  (IOR.EQ. 110)  GC  TO  170 

IF  (IOR.EQ. Ill)  GC  TO  200 

1 6 C RETURN 


DETERMINE  THE  FIRST  ATOM  ANC  THE  LARGEST 
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C NUMBERED  ATOM  IN  EACH  LAYER  FOR  110  LATTICE 

C 

17C  GC  TO  ( 180, 190, 195, 190, 195, 190) , ILAYER 
1 8 C N LA Y= ( ( IX+1  )/2)*(  (IZ+l)/2)+l 
NLAY1=NLAY 
GC  TO  8C 

19C  N LA Y=NLAY  + ( I X/2  ) * ( IZ/2 ) 

GC  TO  8 C 

19  5 N LAY  = NLAY+( ( IX+1 ) /2 ) * ( ( IZ+ 1 ) / 2 ) 

GC  TO  8C 

DETERMINE  THE  FIRST  ATQM  AND  THE  LARGEST 
NUMBERED  ATOM  IN  EACH  LAYER  FOR  111  LATTICE 

2CC  GC  TO  (210,220,230,2^0) , ILAYER 
2 1 C NLAY=( ( LL-2  )/( IY/2) )/2  + 2 
N L AY 1 = N L AY 
GC  TO  80 

22 C NLAY  = < ( LL-2  )/( IY/2)  ) + l 
GC  TO  3 C 

220  NLAY=NLAY1+NLAY-1 
GC  TC  8 C 

240  NLAY=NLAY1+NLAY-1 

c-c  to  ac 

ENC 
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Figure  6 - Microcrystallite  Boundaries 


Figure  9 - Expanded  Impact  Area  - (100)  Orientation 
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Figure  10  - Expanded  Impact  Area  - (110)  Orientation 


Figure  11  - Expanded  Impact  Area  - (111)  Orientation 
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Figure  17  - Angular  Pispersion  - (103),  100-eV,  Cu/I-3'*' 
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Figure  20  - Angular  Dispersion  - (100),  100-eV,  Cu/Ne+ 
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(a) Gibson-II 


(b)KSE-B 


Sputtering  Summary  - (100),  100-eV,  Cu/Ar-*- 


8 


(a) Gibson-II 
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Figure  26  - Sputtering  Summary  - (100)  , 100-aV,  Cu/Cu  + 
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Figure  27  - Ring  Probability  - (100)  , 10 O-eV,  Cu/Cu  + 
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Figure  37  - Angular  Dispersion  - (100)  , 100-?V,  C u/Au* 
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Figure  46  - Sputtering  Summary  - (100),  1-kaV,  Cu/Ne4- 
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?igure  47  - Ring  Probability  - (100),  1-keV,  Cu/Ne+ 
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Figure  48  - Angular  Dispersion  - (100),  1-keV,  Cu/Ns* 
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(a) Gibson-II 
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gure  49  - Spattering  Summary  - (100),  1-k=>V,  Ci/Ar  + 
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Figure  53  - Ring  Probability  - (100),  1-keV,  Cu/Cu-1- 
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Figure  56  - Sing  Probability  - (100),  1-keV,  Cu/Kr+ 
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Ficrure  57  - Angular  Dispersion  - (100),  1-k=V,  Cu/Kr  + 
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(a) Gibson-II 


(b) K3E-B 


Figure  58  - Sputtering  Summary  - (100),  1-ksV,  Cu/7^+ 
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Figure  63  - Angular  Dispersion  - (100),  1-k?V,  Cu/Au* 
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Figure  65  - Ring  Probability  - (100),  1-keV,  Cu/Hg+ 
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Figure  8 1 - Angular  Dispersion  - (111)  , 100-eV,  Cu/Ar  + 
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Figure  83  - Angular  Dispersion  - (111)/  100-?7,  Cu/Ar* 
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Figure  8u  - Sputtering  Summary  - (111),  100-  = '/,  Cu/Ha* 
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igure  85  - Sputtering  Summary  - (111),  100-eV,  Cu/Hg+ 
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Figure  92  - Sputtering  Ratio  vs.  Ion  Atomic  Number 
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Figure  97  - Sputtering  Summary  - (100),  200-eV,  Cu/Cu 
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Figure  100  - Sputtering  Summary  - (100) , 500-eV,  Cu/Cu 
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igure  103  - Sputtering  Summary  - (100),  1-keV,  Cu/Cu 


Figure  106  - Sputtering  Summary  - (100),  100-eV,  Cu/Hg 
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Figure  109  - Sputtering  Summary  - (100) , 200-eV,  Cu/Hg 


Figure  112  - Sputtering  Summary  - (100) , 500-eV,  Cu/Hg 
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F iaura  113  - Atom  Probability  - (100)  , 500-aV, 
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Figure  115  - Sputtering  Summary  - (100),  1-keV,  Cu/Hg 
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Figure  120  - Sputtering  Ratio  vs.  Ion/Target  Atom  Atomic  Mass  Ratio 


Figure  121  - Sputtering  Ratio  vs.  Ion/Target  Atom  Atomic  Mass  Ratio 


[ SPUTS  / RING  ] / [ ATOMS  IN  RING  x TOTAL  SPOTS  ] 


Angle  (Degrees) 

I 1 1 1 1 1—4 1 1 1 1 1— • 1—4 1 1 1 1 

0 5 10  15  20  2 5 3 0 3 5 4 0 4 5 50  5 5 6 0 6 5 70  75  80  85  90 

i 

l 

I 

► 

► 

I 

| 1 1 \ I I I < 1 1 1 1 4-  H I I f-  I 

Heavy  Isotopes 


Light  Isotopes 


Figure  128  - Angular  Dispersion  - Grid  Points 
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